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ABSTRACT 


The  purpose  of  this  work  is  to  describe  a  theoretical  formulation, 
including  a  documented  computer  program,  for  the  evaluation  of  electro¬ 
magnetic  scattering  by  perfectly  conducting  bodies  having  an  axis  of 
rotational  symmetry.  The  main  body  of  the  work  gives  the  theory,  which 
has  been  modified  considerably  from  that  given  earlier.  Appendix  I  gives 
the  analysis  and  logic  which  forms  the  basis  for  the  various  subroutines 
of  the  computer  program.  Appendix  II  gives  the  complete  FORTRAN  listings 
of  the  computer  program.  Finally,  Appendix  III  gives  the  computer  print¬ 
out  for  a  numerical  example,  scattering  by  a  conducting  sphere-cone-sphere 
obstacle  ,  as  obtained  on  the  IBM  7030  digital  computer. 
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SECTION  I 


INTRODUCTION 

GENERAL  DISCUSSION 

In  recent  years  work  has  begun  to  appear  in  the  literature  on  the 
numerical  solution  of  electromagnetic  scattering  problems  by  digital 
computer.  For  the  most  part  these  methods  have  involved  numerical 
solution  of  a  vector  surface  integral  equation.  In  any  case,  the  basic 
procedure  in  all  methods  requires  numerical  generation  of  the  elements 
of  an  N  x  N  matrix,  followed  by  subsequent  Inversion.  Because  N  increases 
roughly  linearly  with  the  size  of  the  target  (quadratically  for  bodies 
that  are  not  axially  symmetric)  there  are  practical  limitations  on  the 
sizes  that  can  be  treated  successfully.  Hence  such  computations,  exact 
in  the  sense  that  in  principle  any  desired  accuracy  may  be  attained, 
are  extremely  useful  in  the  Rayleigh  region  and  some  portion  of  the 
resonance  region but*  must  ultimately  be  supplemented  by  high  frequency 
approximate  techniques  in  order  to  obtain  the  complete  frequency  response 
of  a  given  target.  Observe  that  exact  numerical  computations  may  play  a 
useful  role  in  establishing  the  usefulness  and  accuracy  of  approximation 
techniques,  and  also  in  providing  experimental  targets  for  more  compre¬ 
hensive  range  calibration  than  is  presently  possible  using  spheres  and 
dipoles . 

An  exact  formulation  of  scattering  of  electromagnetic  waves  by 
perfectly  conducting  obstacles  was  given  in  1949  by  Maue ,  who  obtained 
a  pure  integral  equation  and,  alternatively,  an  integro-differential 


1 


equation,  either  of  which  suffices  for  determination  of  the  unknown 

surface  currents  on  the  obstacle. ^  Both  equations  have  been  discussed 

in  the  excellent  review  article  on  diffraction  by  Honl,  Maue ,  and 

(2) 

Westpfahl  (HMW) ,  and  a  derivation  of  the  pure  integral  equation  has 

(3) 

been  presented  by  Van  Bladel.  The  integro-dlfferential  equation  has 

been  programmed  and  solved  numerically  on  the  digital  computer  by 
(4) 

Andreasen,  who  considered  axially  symmetric  targets.  Similarly, 
numerical  analysis  has  been  performed  by  Oshiro  and  co-workers,  employing 
the  pure  integral  equation  for  more  general  shapes. ^ 

An  alternative  theoretical  approach,  also  leading  to  numerical 
results,  has  been  given  by  Waterman. ^  The  purpose  of  this  paper  is  to 
document  a  computer  program  for  the  implementation  of  this  method. 

Section  II  gives  the  theory,  which  has  been  modified  considerably  from 
that  given  earlier.  Appendix  I  gives  the  analysis  and  logic  which  forms 
the  basis  for  the  various  subroutines  of  the  computer  program.  Appendix  II 
gives  the  complete  FORTRAN  listing  of  the  computer  program.  Finally, 
Appendix  III  gives  the  computer  printout  for  a  numerical  example, 
scattering  by  a  conducting  sphere-cone -sphere  obstacle. 

COMPUTATIONAL  ASPECTS 

In  addition  to  their  role  in  the  present  work,  it  should  be 
noted  that  certain  of  the  subroutines  contained  in  this  report  may  be 
of  interest  for  other  applications. 
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Principal  among  these  are  those  routines  for  generating  the 
spherical  Bessel  end  Hankel  functions  by  a  combination  of  power  series 
ana  recursion  techniques,  noting  that  both  precision  checks  and  alter¬ 
native  procedures  are  included  for  those  cases  where  precision  is 
difficult  to  maintain.  The  subroutine  for  generating  associated  Legendre 
functions,  and  their  derivatives,  by  recursion  is  also  essentially  self- 
contained.  Finally,  certain  of  the  matrix  processir.g  operations,  e.g., 
orthogonalization ,  may  prove  of  use  elsewhere,  perhaps  with  modifications. 


SECTION  II 
THEORY 


MAi..j.X  FORMULATION 


Consider  an  incident  electromagnetic  wave  E*(i:)  ,  H*(r)  impinging 
on  the  closed,  perfectly  conducting  surface  0  of  Figure  1  in  otherwise 
free  space.  It  is  assumed  throughout  that  0  is  sufficiently  regular 
that  Green's  theorem  is  applicable,  and  that  0  possesses  a  continuous 

A 

single-valued  normal  n  at  each  point.  Only  simple  harmonic  time  dependence 
at  angular  frequency  U)  is  considered;  a  factor  exp  (-i'JUt)  is  suppressed 
in  all  field  quantities.  Field  behavior  is  described  by  Maxwell's 
equations  in  the  form 


V  X  V  X  e  -  k2E  «  0 , 


with  an  identical  equation  governing  H,  In  these  equations  k-w/c»2TT/\ 

is  the  free-space  propagation  constant. 

Because  the  surface  conductivity  is  infinite  on  0,  no  tangential 

components  of  electric  field  can  be  supported.  Currents  are  induced  in 

the  surface,  the  electric  field  of  which  must  precisely  cancel  the 

tangential  components  of  E*  at  each  point  on  0.  HMW  have  given  a 

representation  of  the  fields  for  this  problem  in  terms  of  surface  current, 

(2) 

After  minor  modification  their  formulas  may  be  written 


E(r)  -  El(r)  +  j  d O'V  x  V  X  i(r')go  (k|r-r'|), 


(2b) 


H(r)  -  H1  (r)  -  ik  Jda’V  x  J.(r')go  (k|r-r’|), 

where  E,  H  is  the  total  field,  gQ(kR)  ■  (4ttR)  1  exp(ikR)  is  the  (scalar) 
free  space  Green's  function  appropriate  to  outgoing  waves,  and  the  curl 
operators  are  with  respect  to  the  unprimed  (i.e.,  not  the  integration) 
variables.  The  integrals  represent  the  E  and  H  fields,  respectively,  due 
to  a  surface  distribution  of  electric  dipoles,  as  one  would  anticipate 
on  physical  grounds.  The  quantity  ^(r.)  ,  which  we  identify  with  Induced 
surface  current,  stands  for  the  jump  discontinuity  in  magnetic  field 
encountered  in  crossing  the  surface,  i.e. 

J.  ■  -  (l/ik)n  x  [l^  -  H_J  on  o  .  (2C) 

In  the  course  of  obtaining  Eqs .  (2),  the  boundary  conditions  appropriate 
to  conducting  surfaces  were  employed,  namely,  that  nxE^  *  n><E_  •  0  on  o. 

The  nature  of  jump  discontinuities  in  the  field  vectors  across  o 
can  be  shown  directly  from  Eqs.  (2),  giving 

£+  '  I.  "  (7»  •  J)  n  -  ni«)p  (3a) 

H  -  H  -  ikn  x  J.  .  (3b) 

— f  — 
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In  the  first  of  these  equations,  the  surface  divergence^  of  the 

current  may  be  defined  by  the  physical  requirement  that  it  equal  the  net 
flow  of  charge  out  of  infinitesimal  element  of  area  (per  unit  area  per 
unit  time).  The  second  equality,  involving  the  surface  charge  density  p, 
follows  from  the  continuity  equation  7S  •  J.  =  -  5p/dt. 

The  extended  boundary  condition,  requiring  that  the  total  electro¬ 
magnetic  field  vanish  identically  in  the  interior  (thus  in  particular 
E  *  0  on  ,  is  from  Eq .  (3a)  sufficient  to  guarantee  the  usual  exterior 

(S 

boundary  condition  n  x  =  0 .  Applying  the  extended  boundary  condition 
in  Eq.  (2a)  gives 

/dCT'7  x  7  x  i(r’)g0  (k|r-r'|)  =  -  El(r)  ,  (A) 

an  "extended"  integral  equation  that  is  to  hold  for  all  points  jr  in  the 
small  dashed  sphere  in  Figure  1.  By  taking  the  curl  of  both  sides  of 
this  equation,  it  follows  that  the  total  magnetic  field  H  will  also 
vanish  in  this  region,  once  Eq .  (A)  is  satisfied.  Equation  (A)  is 
equivalent  to  three  scalar  equations  for  the  two  unknown  tangential 
components  of  J[;  only  two  of  the  equations  are  independent,  however, 
in  consequence  of  the  fact  that  each  side  of  Eq .  (A)  must  have  zero 
divergence  . 

Equation  (A)  may  be  satisfied  by  expanding  both  sides  in  regular 

/Q\ 

vector  eigenfunctions'  1  M„  ,  N  of  the  vector  Helmholtz  Eq .  (1). 

°  -Cmn  -Cmn 
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To  treat  the  integral  one  writes  .  lgQ;  the  expansion  of  the 

"free  space  Green's  dyad"  lgQ  has  been  given  by  Morse  and  Feshbach.  ' 
Because  of  orthogonality  over  any  spherical  surface  about  the  origin 
shown  in  Figure  1,  corresponding  coefficients  may  be  equated  on  both 
sides  of  Eq.  (A)  to  give,  for  incident  plane  E^(r)  ■  e  e*-*1’1^ 

/  *i(£)  • 

-  -  (4TT/ika)in[n(n  +  1)]^  eQ  .  (k)  , 

/ <*ai(r)  .  >f0Bn(r) 

-  +  (An/ik3)  in[n(n  +  1)3*  eQ  .  iBj’Ck)  . 

(5) 

The  If5,  N3  are  the  outgoing  wave  functions,  and  dependence  on  the 

A 

direction  of  incidence  k  is  contained  in  the  vector  spherical  harmonics 
C  °,  B  ^ ^  These  equations  are  to  hold  for  each  triplet  of  values 

(o ,  m ,  n)  ,  with  o  -  e ,  o  (even ,  odd)  ,  m  «  0 ,  1 . n,n-l,2,  .... 

These  are  the  conditions  under  which  the  total  E,  H  field  will  vanish 
identically  in  that  volume  consisting  of  the  largest  sphere  inscribable 
within  o  about  the  coordinate  origin  employed.  As  has  been  shown 
elsewhere, ^  because  of  analytic  continuability  this  is  adequate  to 
guarantee  that  E  and  H  will  vanish  identically  throughout  the  entire 
interior  volume. 
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The  surface  current  is  next  approximaced  by  expansion  in  the 

A  A 

assumed  complete  sef  of  tangential  vector  functions  n  x  M  and  n  x  N; 
one  writes 


i(r)  -  <4/i.c)V  ,m,n,  La  .  •  .r.(r)  x  H,,  i  ,(r) 

—  (  .  u  m  n  u  m  n  — u  m  n 


~  m  n 


m  n 


(6) 


where  the  expansion  coefficients  remain  to  be  determined.  At  this  point 

one  can  expedite  the  discussion  by  introducing  a  matrix  notation.  First, 

the  triplet  of  indices  appearing  in  Eqs.  (5)  and  (6)  are  regrouped  into 

a  single  index  v  by  the  ordering  (Omn)  =  eOl,  oOl ,  ell,  oil,  e02 . 

The  vector  spherical  harmonics  may  then  be  written  as  column  matrices 

C,  B,  having  as  their  vth  elements  (i)n[n(n  +  1)]^  C  G(k)  ,  and 
—  —  inn 

O'  ^  CT  a 

(i)  ln(n  +  1) P  B  (k)  ,  respectively.  The  undetermined  expansion 
mn 

coefficients  of  Eq .  (6)  are  simply  designated  by  the  column  matrices 
a,  b. 

In  this  notation,  substitution  of  the  expansion  Eq .  (6)  into 
Eq .  (5)  yields  a  pair  of  coupled  matrix  equations 


(7) 


9 


for  the  determination  of  a  and  b.  The  matrix  elements  of  1  are  given, 
after  rewriting  the  triple  scalar  product  that  appears,  by 

*W  ’  O'*/")/  •  !f0„n< r)  x  ^,m,n,(r)  ,  (8.) 

and  the  four  matrices  I,  J,  K,  L  differ  from  each  other  only  in  the 
vector  products  appearing  in  the  integrand  of  Eq.  (8)  which  are, 
respectively,  x  x  N^i  ,  x  1^, ,  and  x  , .  By  Inspection 

of  the  integrands,  in  view  of  the  fact  that  M  ■  Re  If*  and  N  ■  Re  N3  ,  it 

is  clear  that  Rel  and  ReL  are  skewsymmetric ,  whereas  ReJ  and  ReK  are 
syranetric.  The  surface  integrals  of  Eq.  (8a)  must,  in  general,  be  done 
numerically  and  are  most  conveniently  performed  in  spherical  coordinates  8, 
0,  for  which  the  appropriate  radial  coordinates  to  employ  may  be  given  by 

the  parametric  specification  r  ■  r  (8,  0)  of  the  surface.  In  view  of 

Green's  second  vector  Identity 


/ 


don  .  [A  X  7  X  b  -  B  X  7  X  a] 


7X7XA-A.7X7XB] 


the  matrices  may  be  seen  to  be  interrelated  by 

K  -  -  J  +  l(D+)_1 

L  -  -  I  ,  (8b) 


10 


where  the  diagonal  matrix  D+  (and  D  ,  employed  below)  has  vth  elements 
defined  by 


<Vv 


£  (2n  +1)  (n  -  m) ! 
(+  l)n  -2 - - 

4n(n  +1)  (n  +  m) 


(8c) 


The  Neumann  factor  £  is  given  by  £  ■  1,  £  -2  otherwise. 

It  is  also  desired  to  compute  the  scattered  field  if  ,  if  given 
by  the  surface  integrals  in  Eq.  (2).  Specifically  for  the  electric 
field,  one  has 

E5  (r)  =  4  Y~  [f_  M3.  (r) 

—  — '  amn—  Omn  — 

0mn 


max  on  O 


F(k  r 

—  out,  in 


kr  »  1 . 


(9) 


The  vector  scattering  amplitude  F,  depending  both  on  direction  of  incidence 

k,  and  observation  k  ,  is  obtained  by  introducing  asymptotic  forms 
in  out’ 

of  the  outgoing  partial  waves  M3  ,  K3  in  the  preceding  expression  for  §f 
to  get 


F(k 


out 


,£in> 


(4/ik) 


[C'(k  ) 

-  out 


D  f  +  iB' 


(k  )D 
out  ■ 


g] 


(10) 
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where  C'  is  the  transpose  of  C  (and  hence  a  row  matrix) .  The  outgoing 
partial  wave  expansion  coefficients  f,  g  are  expressed  in  terms  of 
surface  currents  a,  b  by 


(ID 


These  formulas  have  been  obtained  by  employing  that  expansion  of  the 
free  space  Green's  dyad  valid  in  the  exterior  region  outside  the  large 
dashed  sphere  of  Figure  1. 

scat  (11) 

The  scattering  cross  section  O  is  given  by 

o3Cat  =  (16TT/ka)  (f'*D  f  +  g'*D+g)  .  <12a) 

As  a  numerical  check  on  accuracy,  one  may  also  compute  the  total  cross 
section 


otot  -  <4TT/k)  Im  [e  •  F  (£  ,  k  )  ]  ,  (12b) 

o  —  in  in 


scat  (11) 

which  must  equal  0  by  the  forward  amplitude  theorem  .  The  radar 

cross  section,  defined  as  4n  times  the  back-scattered  power  per  steradian 

divided  by  incident  power  per  unit  area,  is  given  by 


0radar  -  (64TT/ka)  |eQ  •  C'  (-kln)D.f 


+  lao  •  S'  (^in)D-g|2 


(12c) 
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If  che  return  signal  is  regarded  as  resolved  into  two  orthogonal 

linearly  polarized  modes,  then  this  equation  gives  a  measure  of  the 

power  carried  in  that  mode  having  polarization  aligned  with  the  original 

incident  wave,  whereas  the  cross-polarized  return  is  given  by  replacing 

e  by  e  '  ■  k,  *  e  in  Eq.  (12c). 
o  J  o  in  o 

EVALUATION  OF  THE  TRANSITION  MATRIX 

Instead  of  first  solving  Eq .  (7)  for  the  currents  a,  b,  then 
substituting  in  Eq .  (11)  to  obtain  the  scattered  wave  f,  g,  the  currents 
may  be  formally  eliminated  to  obtain  the  scattered  wave  directly  from 
the  incident  wave  as 


The  block  matrix, 


is  known  as  the  transition  matrix,  and  is  both  symmetric  (i.e.  T^'  ■  T^, 

* 

T2'  -  T3>  T^'  -  T^)  and  has  the  property  T  'T  ■  ReT,  i.e., 


jT1*T1  +  T2*T3  Tx*T2  +  T2*T4^ 

fT  1  T2\ 

- 

-  Re 

1 

+  t4*t3  t3*t2  +  t4*t  J 

\t3  t4/ 
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(15) 


The  property  Eq.  (lb)  is  a  consequence  of  unitarity  of  the  scattering 


matrix  S  *  1-2T,  as  may  be  verified  by  substitution  in  the  unitarity 
condition  S'*S  -  1. 


If  one  now  defines  the  matrix  Q  as 


(16) 


then  by  comparison  with  Eqs.  (7)  and  (11)  the  transition  matrix  is 
determined  by  the  matrix  equation 

QT  -  Re  Q  ,  (17) 


which  in  general  must  be  solved  numerically. 

Instead  of  working  with  the  2  by  2  block  form  of  Eq.  (17) , 
involving  in  truncation  four  N  x  N  matrices,  it  is  convenient,  for  the 
numerical  processing,  to  change  over  to  single  2N  x  2N  matrices.  Thus, 
define  the  2N  x  2N  matrix  Q  by 


* 

Q(2m-1)  (2n-l) 

■  «i>- 

\ 

a 

Q(2m-1)  (2n) 

► 

^(2m)  (2n-l) 

■  «3>- 

A 

Q(2m)  (2n) 

'  <Vmn 

i 

N  . 


(18) 
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The  matrices  T,  and  S  «  1  -  2T  are  defined  in  exact  analogy  to  this. 

At  this  point,  Eq.  (16b)  may  be  written  in  terms  of  S  as 

QS  -  -  Q*  .  (19) 

Because  of  the  behavior  of  the  radial  (Kankel)  functions  that  appear  in 
the  matrix  elements  of  Q,  the  imaginary  parts  of  the  elements  of  Q  will 
tend  to  grow  very  large  numerically  above  the  diagonal.  In  order  to 
avoid  loss  of  precision  due  to  the  finite  precision  arithmetic  employed 
by  the  digital  computer,  it  is  convenient  at  :hls  stage  to  reset  all  the 
mentioned  elements  to  zero,  by  Gaussian  elimination.  This  process  is 
straightforward,  the  net  effect  being  to  premultiply  Q  by  a  real  upper 
triangular  matrix  (all  elements  zero  below  the  main  diagonal).  Suppose 
this  conditioning  to  have  been  performed  on  Eq.  (19)  ,  which  we  continue 
to  employ  without  change  of  notation. 

To  Eq.  (19)  are  adjoined  the  constraints  of  symmetry  and  unltarlty 

A 

mentioned  above,  which  are  unaffected  by  the  S  ■*  S  transformation  and 
thus  given  by 

F  -  S’  (20) 

and 

§'*§-!  .  (21) 
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Two  extremes  of  view  with  regard  to  the  system  of  Eqs.  (19),  (20)  and 
(21)  are  as  follows:  first,  one  might  truncate  the  matrix  Eq.  (19), 
solve  numerically  by  digital  computer,  then  compare  the  resulting 
solution  with  Eqs.  (20)  and  (21',  the  latter  thus  being  employed  as 
consistency  checks.  On  the  other  hand,  one  might  attempt  to  treat  all 
three  equations  from  a  unified  point  of  view  from  the  onset,  obtaining 
a  solution  in  some  sense  of  Eq.  (19)  subject  to  the  constraints  of 
Eqs.  (2rt)  and  (21).  The  first  approach  has  been  employed  in  earlier 
work  on  the  computer  for  bodies  of  rotational  symmetry,  and  works  quite 
satisfactorily  for  a  restricted  range  of  body  shapes  and  sizes.  The 
second  approach  is  employed  in  the  present  work  in  order  to  extend  the 
range  of  bodies  that  can  be  handled,  in  view  of  the  fact  that  the 
constraints  essentially  determine  three  quarters  of  the  solution 
[i.e.,  of  the  8N^  real  parameters  appearing  in  the  2N  *  2N  (truncated) 
complex  matrix  S,  it  can  be  shown  that  only  N(2N  +  1)  are  independent, 
if  §  satisfies  Eqs.  (20)  and  (21)]. 

A 

To  develop  a  unified  analysis,  observe  first  that  if  S  could  be 
constructed  in  the  form 


§  -  U'  U 


(22) 


where  U  is  unitary,  then  both  constraints  would  be  satisfied  by  inspection. 
This  suggests  that,  rather  than  inverting  Q  directly  in  Eq.  (19),  it  be 
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made  unitary.  Inus  consider  the  upper  triangular  matrix  M  (i.e.,  all 
elements  are  zero  below  the  main  diagonal)  which  by  premultiplication 

A  A 

makes  Q  into  a  unitary  matrix  Q  ,  ,  viz. 

unit 


MQ 


Qunit 


(23) 


Premultiplying  Eq .  (19)  by  H,  one  can  write 


0  S 
\init 


-  MQ* 


MM*"1  Q* 


unit 


Upon  solving  for  S,  there  now  results 


S  -  -  Q'*  .  (m*"1)  Q*  .. 

^  unit  unit 


(24) 


Substituting  this  result  in  Eq .  (20),  the  symmetry  constraint,  it  follows 
without  difficulty  that  the  matrix  product  MM*  1  must  be  symmetric.  But 
each  of  the  matrices  appearing  in  the  product  is  upper  triangular,  and 
their  product  is  again  upper  triangular.  Consequently  the  product  must 
be  a  diagonal  matrix.  Further,  the  diagonal  elements  can  be  written 
out  explicitly,  giving 


MM* 


-1 


M  /M  * 

nir  ii 

o 

o 


M22/M22* 

0 
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If  next  one  can  arrange  to  choose  the  diagonal  elements  of  M  to  be  real, 
then 


MM*"1  *  1 


From  Eq.  (24)  the  S-matrix  is  now  given  by 


S  =  -  O'*  0* 

w  unit^  unit 


(25) 


(26) 


which  is  of  the  required  form  Eq.  (22).  Substituting  Eq.  (26),  along 

with  the  identity  O'*  .0  .  ■  1,  back  in  the  relation  S  **  1  -  2T,  the 

unit  unit 

desired  transition  matrix  is  given  by 


T 


Q'*  .  Re 

unit 


(<Wt> 


(27) 


and  the  block  form  of  T  is  readily  obtained  by  reversing  the  transformation 
of  Eq.  (18)  . 

Returning  to  M  for  a  moment,  Eq .  (25)  states  simply  that  M  is 
real.  Thus  the  process  may  be  summed  up  in  the  (formal)  theorem:  given 
the  matrix  Eq.  (19),  with  constraints,  Eq .  (20)  and  (21),  on  the  solution, 
it  follows  that  the  given  matrix  Q  cannot  be  arbitrary,  but  must  be  such 
as  to  be  factorizable  into  the  product  of  a  real  upper  triangular  matrix 
and  a  unitary  matrix,  namely 


A 

Q 


unit 


(28) 
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The  transformation  of  Q  into  a  unitary  matrix,  as  required  in 
Eq.  (23),  is  done  by  Schmidt  orthogona lization  of  the  2N  vectors  given 

A 

by  the  rows  of  Q,  beginning  with  the  bottom  row  and  working  up.  The 
pzocedure  is  straightforward,  and  details  are  described  in  a  subsequent 
section . 

APPLICATION  TO  SPECIAL  GEOMETRIES 

In  order  to  apply  the  equations  to  bodies  having  an  axis  of 
rotational  symmetry,  the  axis  of  symmetry  is  chosen  as  polar  axis  for 
our  spherical  coordinates  and,  without  loss  of  generality,  the  direction 

A  A 

of  incidence  taken  in  the  plane  of  azimuth  v  ■  0,  so  that  k^n  ■  k^(u,  0) . 
A  reduced  index  notation  may  be  employed  for  those  matrix  elements  that 
do  not  vanish  under  the  azimuthal  integration,  writing 


I  ,  E  I  ,  -  -  I 

mnn  omnemn  emnomn 


^mnn '  ^emnemn '  ^omnomn ' 


K  ,  -  -  J  ,  +  16  ,  D 

mnn  mnn  nn  mnn 


-1 


L  ,  =  I 
mnn  mnn 


(29) 
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The  independent  matrix  elements,  written  out,  are 


mnn1 


rr 

m  I  d0  te  (Pn  Pn()  (kr)2  hn(kr)  jn'  (kr) 
0 


J 


mnn' 


TT 

J  dQ  sin  0 
0 


„m 


2_m 

m  P  r  , 
n  n 

sins  0 


dpm,  ap1 


m 


n 


de  ae 


krh  (kr)  — —  ikrj  , (kr)  ] 
n  dkr  Jn 


(30) 


—  n '  (n 1 

e 

m 


« 


+  1)  J  d0  sin 
0 


ap 

_ n 

ae 


m  „m 


h  (kr) 
n 


V(kr)  ae  (kr) 


Observe  that  the  real  parts  of  all  these  matrices  are  symmetric.  Also, 
because  of  the  vanishing  of  all  matrix  elements  with  different  azimuthal 
mode  indices  (m  t  m')  ,  there  is  no  coupling  and  each  azimuthal  mode 
m  =  0,  1,  2,  ...  may  be  evaluated  separately. 

From  the  defining  Eq.  (16)  ,  the  only  non-vanishing  elements  of 
the  Q  matrix  may  now  be  written  in  reduced  index  notation  as 


«!>«.• 

^l^emnemn' 

^1^  omnomn' 

(Q2)mnn' 

^2^  omnemn 1 

^2^  emnomn' 

«3>mnn' 

^3^  emnomn' 

^3^omnemn' 

^4^  mnn' 

^4^emnemn' 

=  -  (Q.) 

4  omnomn 
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In  addition,  the  reduced  index  elements  are  related  by 


WjU'  ■  (Vmnn’ 


(Qjfflnn,  -  -  (Q,)mnn.  +  i  1 

4  mnn  l  mnn  = 


(32) 


Finally,  examination  of  Eq .  (17)  reveals  that  the  non-vanishing  elements 
in  the  four  blocks  of  the  T  matrix  are  interrelated  exactly  a9  in 
Eq .  (31),  but  not  Eq .  (32),  so  that  the  complete  solution  may  be  obtained 
by  solving  Eq .  (17)  once,  using  the  reduced  index  quantities. 

A  further  important  reduction  occurs  in  the  preceding  equations 
for  obstacles  (e.g.,  finite  cylinder)  having  a  plane  of  mirror  symmetry 
normal  to  the  axis  of  rotational  symmetry.  For  this  geometry  the  radius 
vector  r(0)  specifying  the  shape  of  the  obstacle  will  be  even  about 

0  «  tt/2,  i.e.  , 


r (6)  -  r(n  -  0)  .  (33) 

Inspection  of  the  parity  of  the  integrands  giving  rise  to  matrix  elements 
in  Eq .  (30)  readily  reveals  that  a  checkerboard  pattern  of  zeros  has 
emerged ,  i.e., 

I  ,  ■  0;  (n  +  n')  even 
mnn 

J  ,  -  0;  (n  +  n')  odd  .  (34) 

mnn 

These  elements  can  hence  be  set  to  zero  without  performing  the  numerical 


integrations . 
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Prolate  (and  oblate)  spheroids  have  a  mirror  symmetry  plane 
normal  to  their  rotational  symmetry  axis,  so  that  both  mode  and  parity 
decompositions  may  be  made,  as  discussed  above.  There  is  another 
reduction  that  occurs  here,  however,  which  from  a  theoretical  standpoint 
lays  the  Rayleigh  expansion  out  in  full  view,  and  for  numerical  purposes 
yields  extremely  well-conditioned  matrices  for  inversion. 

To  see  this,  let  us  examine  the  matrix  elements  as  given  by 
Eqs .  (30).  The  numerical  magnitude  of  these  elements  is  influenced 
mainly  by  the  radial  functions  appearing  in  the  integrand.  For  I  , , 
for  example ,  one  has 

W  ~  (kr>a  h„<kr>  V(kr)  B  <kO*  Cj„Jn,  ►  l  unin,]  . 

For  a  given  argument  x,  the  Bessel  functions  Jn(x)  decrease  rapidly  in 
magnitude,  and  the  Neumann  functions  nn(x)  increase,  roughly  as  soon  as 
the  index  n  exceeds  x.  Thus  the  real  part  of  I,  which  is  obviously 
symmetric,  will  eventually  decrease  rapidly  in  magnitude  as  one  proceeds 
along  any  row  or  column.  The  numerical  behavior  of  I  is  dominated  by  its 
imaginary  part,  for  which  elements  again  decrease  going  out  any  row,  but 
increase  going  down  any  column,  at  such  a  rate  that  diagonal  elements 
remain  relatively  constant.  These  large  numerical  values  presumably 
strongly  influence  the  truncated  matrix  inversion  procedure. 
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One  can  show,  however,  that  for  prolate  or  oblate  spheroids 

this  behavior,  specifically  the  arbitrarily  large  values  by  which 

elements  of  I  below  the  diagonal  exceed  corresponding  elements  above, 

vanishes  identically.  I  and  J  become  completely  symmetric,  and  dominant 

terms  lie  only  on  the  diagonal  once  either  row  or  column  index  exceeds 

kr  ,  where  r  is  the  radius  of  the  circumscribing  sphere, 
max  max 

(12) 

Based  on  results  given  by  Watson  one  can  show  that  the  radial 

factor  for  an  element  below  the  diagonal  in  the  imaginary  part  of  I  , 

mnn 

is  of  the  form 


nn+2  s+ 1  ^  n 


=  X2"nWl  + 


2s 


2s-2 


■-  +  1 


(35) 


where  the  equivalence  symbol  (-)  indicates  that  the  exact  coefficients 

of  inverse  powers  of  x~  have  not  been  included,  as  they  are  not  required 

in  the  present  discussion.  The  first  term  on  the  right-hand  side 

corresponds  precisely  to  the  symmetrically  placed  element  above  the 

2 

diagonal;  we  must  thus  show  that  the  inverse  powers  of  x  contribute 
nothing  to  the  integral 

TT 

Im^Im(n+2s+l)„]  '  »  (  I?  <C2s+lO  <kr>3V2s+l(kr>  Vkr>  ■  <36) 
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For  a  prolate  (oblate)  spheroid,  having  semi-axes  a,  b,  one  has 


kr  ■  ka  [cos29  +  (a/b)3  sin29]  ^ 


(37) 


which  may  be  rewritten  (identifying  x  with  kr) 


1/x2  ip  +  P, 
o  * 


(38) 


Now  the  product  of  two  Legendre  functions  may  itself  be  expanded 

in  a  series  of  Legendre  polynomials,  with  indices  ranging  from  the 

(13) 

difference  to  the  sum  of  the  original  indices,  i.e., 


P 


m 

n 


n+n' 


4  I 

p-n-n' 

(p+n+n' 


P 

P 

even) 


(39) 


where  again  explicit  numerical  coefficients  have  been  ignored. 
Substituting  Eq.  (38)  in  the  series  of  inverse  powers  of  x‘  appearing 
in  Eq.  (35)  ,  then  employing  Eq .  (39).  repeatedly,  one  can  write 


2  i 

x  n  .  1 
n+2s+l  n 


X  nn^n+2fi+l 


I 


2q 


(40) 
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This  result  may  be  put  in  Eq .  (36),  recalling  also  that  Rel  is  symmetric, 
to  get 


TT 


*m(n+2s+l)n  *mn(n+2s+l) 


m 


/  d0  §9 


s 


Y  p 


m 

n+2s+l‘n/  L-  *2q 
q«0 


n 


'  -/ 


m 


d0  sin  0  P”  ^.P 
n+2s+l 


mY? 

n  L  2q-l 

q-i 


TT 


-m 


s+n  s 

de  sin  e  Y  P2ptl  Y  P2,.l 

p-S  q-1 


J  d0  sin  0 


(41) 


where  in  the  second  step  we  have  integrated  by  parts  ,  then  employed 

Eq.  (39),  and  finally  observed  that  the  highest  Legendre  polynomial 

appearing  in  the  second  sum  is  P  . ,  while  the  first  sum  begins  at 

Z  S“  i 

P  ;  because  of  orthogonality,  all  the  resulting  integrals  vanish. 

2  s+1 

To  perform  the  analogous  calculation  for  J  n,  »  one  proceeds  by 
first  employing  Green's  identity  to  rewrite  Jmnn!  the  more  symmetric 
form 

TT 

J„nn'  •  -  r  /d6  sIn  8  Bnnn'(6>  55T  “X00  V  00  ]*-kr<e> 

m  ' 


sin  0  Cmnn,(6)[x3hn(x)jn,(x)d(l/x2)/d0]x-kr(9) 

(42) 
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valid  for  n  ^  n1  ,  with 


®  1 
mnn 


2  _m  _m 
m  P  P  , 
n  n 

sin2  9 


dpm  apm, 
n  n 

“51T  ae 


(43a) 


C  nn'  <9> 

mnn 


dPm  dpm, 

n'  (n'+i)  jgS-  Pn,  +n(n+l)  P^ 


(43b) 


It  is  convenient  this  time  to  write 


1/x2  -  const.  +  sin2 9 


Using  this  in  conjunction  with  the  inverse  polynomial  expression 


X  nn+2sJn  *  X  nnJn+2s  +  2s-l  +  2s-3  +  +  x 


Eq .  (42)  may  be  reduced  to 


^m(n+2s)n  ^mn(n+2s) 


TT 

’  <2s‘1)/€m  / d6  sin  6  Bu(n+2s)-  )  <•!"» 

0 


a 

..I 

q-l 


2q 


rn  flil 

+  1/£»  J  d6  sl"S  0  coa  0  I  <sln  0)2,1 

q-0 


(44) 
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where  in  the  first  term  the  constant  term  in  the  summation  has  been 
dropped  because  of  the  additional  orthogonality  relations 


TT 


/ 


d9  sin  9  B  , 
mnn 


0,  n  +  n' 


At  this  point,  using  the  standard  recursion  formulas  for  the 
Legendre  functions  one  can  write 


sin2 9  B 


m(n+2s)n 


=  P 


m  pm 

n+2s-l  n+1 


+  P 


m 

n+2s- 


1 


+  P 


m  m 
n+2s  n 


,m 


m 


m 


m 


+  P  P'"  f.  P”*  P" 

n+2s+l  n-1  n+2s+l  n+1 


(45a) 


sin  9  cos  9  C 


i  pm  Pm  +  Pm  Pm 

m(n+2s)n  n+2s-l  n+1  n+2s-ln-l 


,  _m  „m  m  „m 

P  P  -f*  P  P 

n+2s+l  n-1  n+2s+l  n+1 


(45b) 


The  polynomials  in  sin2 9  appearing  in  Eq .  (44)  may  be  expanded  in 
Legendre  polynomials  of  highest  index  2(s-l)  [Note  that  a  factor  sina9 
has  been  taken  out  in  the  first  case  to  employ  in  Eq ,  (45e)].  By 
examination  of  Eq.  (39)  it  may  be  seen  because  of  orthogonality  that 
only  the  first  term  on  the  right-hand  side  of  Eqs .  (45a)  and  (45b) 
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will  make  a  non-zero  contribution  to  their  respective  integrals. 
Writing  out  these  non-vanishing  terms  in  Eq .  (44)  explicitly,  one 
finally  obtains 


^m(n+2s)n  ^mn(n+2s) 


TT 


(2s-l) (n+2s+l) (n+2s+m) (n-nri-1) 


e  (2n+4s+l) (2n+l) 
m 


- 

de 


.  /  i  n\2s*"2 

Sln  D  Pn+2S-lPn+l(sln  0) 


+  same  expression 


=  0 


(46) 


Thus  I  and  J  are  symmetric,  and  one  need  only  compute  elements 

on  and  above  the  diagonal  in  Eqs.  (30).  The  matrices  are  expected  to 

be  well-conditioned  in  the  sense  that  numerical  results  will  converge 

rapidly  to  final  values  versus  truncation.  From  the  point  of  view  of 

the  Rayleigh  expansion  in  powers  of  ka ,  valid  at  low  frequencies,  observe 

that  all  matrices  may  be  expanded  in  powers  of  ka ,  e.g.,  writing 

(J  m^)  ,  ■  J  ,  one  has 

nn  mnn 


=  A  +  B(ka)  +  C(ka)2  +  ... 

=  A  [l  +  A'^J^-A)]  (47) 
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where  A  is  diagonal,  B  is  tridiagonal  (all  elements  zero  except  on, 
one  above,  and  one  below  the  diagonal)  and  so  forth.  The  inverse, 
expanded  in  powers  of  ka ,  is  readily  obtainable  by  the  binomial  theorem 
as 


(J(m))_1  **  [l  -  A"1  (J(m)-A)+  ...  ]  A_1  .  (48) 

INTERPRETATION  OF  NUMERICAL  RESULTS 

In  order  to  provide  some  insight  into  the  behavior  in  practice 
of  the  various  matrices  discussed  above,  the  numerical  printout  for  an 
example  has  been  included  (Appendix  III)  .  In  addition  to  providing  a 
test  case  for  use  with  the  computer  program,  many  features  of  matrix 
behavior  are  most  conveniently  described  by  reference  to  this  printout. 

The  obstacle  to  which  the  results  refer  consists  of  a  sphere- 
cone-sphere  ,  as  shown  in  Figure  2.  The  analytical  description  of  this 
shape  as  inputted  to  the  computer  is  detailed  in  Paragraph  9.0  of 
Appendix  I.  It  will  be  seen  that  the  printout  consists  almost  entirely 
of  matrix  quantities,  as  an  aid  in  gauging  the  numerical  effectiveness 
of  the  truncation  being  employed. 

The  first  page  lists  input  parameters.  Thus,  four  cases  (m  *-  0, 

1,  2,  3)  were  evaluated  consecutively,  with  truncation  to  6  x  6  matrices. 
The  body  is  described  analytically  in  three  sections.  The  body  shape 
"9"  indicates  that  the  body  does  not  possess  mirror  symmetry  normal  to 
the  axis  of  rotational  symmetry.  U  vector  indicates  that  46  aspect 
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Figure  2.  Geometry  of  the  Sphere-Cone- Sphere  Body  Showing  E  I !  Incident  Polarization 


angles  will  be  evaluated.  The  body  size  is  kA  =  1.0  (Figure  2),  and 

the  ratio  of  spherr  radii  is  b/a  =  (1  +  sinO.)  *  ^  0.794  (correct  to  nine 

significant  figures  in  the  computer),  with  cone  half-angle  (X  =  15  degrees. 

Numerical  integration  is  performed  by  Bode's  rule  using  64  equally  spaced 

divisions  in  each  section,  and  the  angular  end  points  in  degrees  are 

(14) 

indicated  for  each  section. 

The  basic  quantities  shown,  for  each  m  value  successively,  are 

A 

the  Q  matrix,  the  orthogonalized  Q  matrix,  the  transition  matrix  T, 
and  the  cumulative  cross  section  quantities. 

Consider  first  the  case  m  **  0.  The  blocks  and  are  zero 
for  this  case,  and  hence  are  not  shown.  The  remaining  blocks,  and  , 
are  obtained  by  numerical  integration  from  the  defining  equations 
(16,  29,  30,  31,  32).  Because  of  the  vanishing  of  the  blocks  Q2  , 
it  turns  out,  as  one  can  verify  with  some  study,  that  the  remaining 
two  blocks  actually  are  processed  with  no  interactions,  so  that  behavior 
can  be  discussed  by  examining  say,  ,  alone. 

Considering  the  imaginary  part  of  ,  which  is  the  numerically 
dominant  portion,  one  observes  that  elements  of  each  of  the  six  row 
vectors  increase  in  numerical  magnitude  as  one  moves  to  the  right.  It 
is  immediately  clear  that  row  vectors  must  be  orthogonalized  from  the 
bottom  up,  in  order  that  elements  of  the  resulting  vectors  may  settle 
down  to  constant  values  independent  of  truncation  (that  is,  if  one  began 
orthogonalizing  with  the  top  row,  then  it  can  be  seen  that  after 
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normalizing,  the  first  element  is  smaller  by  a  factor  14.1/11.1  ~  13 

than  it  would  have  been  in  5  x  5  truncation).  Observe  also  that  the 

bottom  row  is  the  best  candidate  for  a  "unit  vector  in  the  six  direction," 

in  that  the  sixth  element  is  larger  than  the  first  by  a  factor  of  about 
7  2 

10  ,  whereas  the  corresponding  factor  for  the  top  row  is  only  2  x  10  . 

The  orthogonalized  matrix  Q  unit,  Eq .  (23),  is  shown  next,  having 

dimension  12  x  12.  Odd  numbered  rows  have  come  from  the  original 

6  x  6  ,  whereas  the  even  numbered  rows  are  associated  with  .  It  is 

striking  to  observe  that  each  of  the  original  matrices  (and  hence  the 

entire  imaginary  part  of  Q  ,  )  has  become  nearly  diagonal.  This 

unit 

occurred  because  the  main  effect  of  the  sixth  row  vector  was  to  reduce 

the  last  entry  of  each  preceding  row  to  nearly  zero.  The  new  fifth 

row  vector,  in  similar  fashion,  then  served  primarily  to  reduce  the 

fifth  entry  in  each  of  the  preceding  four  rows,  and  so  on.  Thus  the 

first  row  of  ,  which  originally  increased  by  a  factor  of  about  170 

from  first  to  last  entry,  now  decreases  by  a  factor  (see  row  one  of 

-  ,  -8 

the  Imaginary  part  of  Qun^t)  of  about  7  x  10  .  The  total  relative 

reduction  is  of  order  4  x  10 

It  is  not  difficult  to  study  the  behavior  of  the  row  vectors,  or 

the  individual  matrix  elements,  versus  truncation.  For  example,  if  a 

5x5  truncation  had  been  employed,  then  the  first  entry  in  the  fifth 

row  vector  would  have  been,  after  normalization,  (-1.278  +  i0.5861)  x 

10  "V (0.5646)  -  (-2.263  +  il.038)  x  10  \  whereas  the  6x6  truncation 

(see  row  nine  of  6  .  )  actually  gives  (-2.251  +  il.036)  x  10  One 

unit 
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can  verify  that  the  sixth  row  vector  has  even  less  effect  on  rows 
earlier  than  the  fifth.  In  particular,  orthogonalizing  the  first  row 
to  the  sixth  row  will  change  the  first  element  of  the  first  row  from 
0.848  to  approximately  0.848  -  (144/0.551)  (6.99  x  10  )  ,  which 

constitutes  a  change  in  the  sixth  significant  figure. 

The  blocks  of  the  transition  matrix,  computed  from  Eq .  (27)  and 
the  reverse  transformation  of  Eq.  (18),  are  shown  next.  Again,  the 
blocks  T2  =  T^'  vanish  identically  and  are  not  shown.  Both  and 
are  seen  to  be  exactly  symmetric  to  the  number  of  digits  given,  and 
the  unitary-related  condition  of  Eq .  (15)  is  readily  verified  on  the 
desk  computer  to  within  round-off  error  in  the  last  place  shown.  Observe 
that  the  elements  of  both  and  fall  off  rapidly  in  magnitude  moving 
away  from  the  upper  left  hand  corner,  so  that  the  scattering  behavior 
would  be  efficiently  and  accurately  describable  in  this  instance  using 
only  the  first  two  rows  and  columns  of  and  T^ . 

Finally,  the  accumulated  (over  m  =>  0  only)  far  field  quantities 
are  shown  for  G  ||  (class  1)  and  EX  (class  2)  polarizations.  The  first 
column  gives  the  incident  aspect  angle  measured  from  the  axis  of 
rotational  symmetry.  For  each  aspect,  subsequent  columns  give  the 
scattering  cross  section  [Eq .  (12a)],  forward  amplitude  [the  complex 
quantity  appearing  in  Eq .  (12b)],  backscattered  amplitude  [the  complex 
quantity  appearing  in  Eq .  (12c)  before  squaring],  and  finally  the  radar 
cross  section  and  phase  of  the  back  scattered  amplitude.  Observe  that 
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both  energy  conservation  (equality  of  the  second  and  fourth  columns) 
and  reciprocity  (symmetry  of  the  third  and  fourth  columns  about  the 
aspect  angle  of  ninety  degrees)  are  satisfied  to  seven  significant 
figures . 

Turning  to  the  case  m  -  1,  the  blocks  and  are  shown, 
and  (?4  then  being  given  by  [see  Eqs  .  (8b,  16,  29)]  Q3  *  Q2  ’  ^4  =  "^1  + 

A  partial  check  on  the  precision  of  numerical  integration  is  available 
for  this  and  all  subsequent  m  values.  From  the  Wronskian  relation 
x  [j  (x)h  ,  (x)  -  j  ,(x)h  (x)]  *  i  and  the  first  of  Eqs.  (30)  it  is 
immediately  clear  that  the  imaginary  parts  of  the  first  off-diagonal 
elements  of  should  be  symmetric,  a  result  not  used  in  the  program 
(whereas  symmetry  of  Ret^  and  is  always  enforced) .  The  expected 

symmetry  is  seen  to  obtain  to  seven  significant  figures  for  the  (1,  2) 
and  (2,  1)  elements.  Precision  subsequently  deteriorates  slightly  so 
that  discrepancies  have  appeared  in  the  fifth  significant  figure 
between  the  (5,  6)  and  (6,  5)  elements. 

Numerical  behavior  of  both  the  Q  and  the  orthogona lized  Q  matrices 
appears  to  proceed  substantially  as  in  the  case  discussed  above  for 
m  ■  0,  although  the  details  are  of  course  considerably  more  complex 
because  of  the  presence  of  all  four  non-zero  blocks.  The  near-diagonal 

A 

nature  of  the  orthogonalized  Q  matrix  is  still  evident  by  inspection, 

however.  In  the  resulting  transition  matrix  the  blocks  T^  and  T^  are 

seen  to  be  symmetric,  and  the  block  T^  to  be  equal  to  the  transpose  of 

.  Comparison  of  the  far  field  results  with  those  for  m  ■  0  reveals 

significant  changes  at  all  aspect  angles. 
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For  the  two  subsequent  cases  m  =  2 ,  3a  new  effect  is  seen  over 
and  above  the  previously  discussed  features,  due  to  the  vanishing  of 
the  associated  Legendre  functions  P™  for  n  <  m.  In  consequence,  the 
first  row  and  column  of  each  block  of  the  Q  matrix  (for  m  =  3  the  first 
two  rows  and  columns)  are  identically  zero.  The  behavior  versus  truncation 
at  n  =  6,  as  judged  by  the  near-diagonal  results  after  orthogonalization , 
appears  unaffected,  however.  The  net  result  is  that  the  computation 
becomes  gradually  simpler  as  m  increases,  m  =  2  requiring  treatment  of 
5x5  blocks,  and  m  =  3  requiring  only  4x4  blocks. 

A  measure  of  error  incurred  by  stopping  at  m  *  3  may  be  obtained 
by  comparing  the  far  field  results  with  those  obtained  at  m  ■  2  (except 
for  incidence  along  the  axis  of  rotational  symmetry,  0  degrees  or  180 
degrees,  for  which  scattering  behavior  is  completely  determined  from 
the  m  =  1  results  only).  At  incidence  80  degrees  from  the  axis,  for 
example,  and  for  either  polarization,  the  scattering  cross  section  is 
seen  to  be  unchanged  to  about  six  significant  figures.  For  the  same 
cases  the  radar  cross  section,  however,  has  changed  in  about  the  third 
significant  figure.  Such  precision  is  nevertheless  quite  adequate  in 
most  practical  applications. 


35 


BLANK  PAGE 


APPENDIX  I 


ORGANIZATION  OF  THE  COMPUTER  PROGRAM 
1.0  INTRODUCTION* 

The  "IMSCAT"  Program  has  been  written  in  FORTRAN  IV  language 
for  the  IBM  7030  Computer  to  produce  solutions  to  the  electromagnetic 
scattering  problems  which  are  outlined  above  in  Section  II.  Several 
factors  were  given  consideration  in  the  design  of  the  program: 

Efficient  coding  to  reduce  computer  run  time  as  much  as  possible. 
The  routine  VECMUL  for  matrix  by  vector  multiplication  was  coded  in 
machine  language  to  take  advantage  of  specialized  coding  available 
at  that  level.  This  routine  is  also  available  in  FORTRAN  (though 
less  efficient  and  accurate)  so  that  the  program  can  be  run  on 
machines  other  than  the  7030  Computer. 

Full  single  word  accuracy  of  a  7030  register  and  where  necessary 
double  precision  accuracy  was  utilized  in  the  calculation  of  special 
functions.  Single  precision  accuracy  on  the  7030  maintains  15  digits 
of  accuracy. 

Maximum  use  of  core  storage.  The  size  of  the  solution  matrices 
(60  x  60  complex)  was  determined  so  that  secondary  storage  devices 
such  as  tapes  do  not  have  to  be  utilized  in  running  the  program. 

The  matrices  are  stored  and  manipulated  from  1  of  3  large  blocks 
of  common  storage.  Each  block  is  dimensioned  120  x  120.  However, 
through  various  equivalence  statements  in  the  proper  routines,  these 
major  blocks  are  resegmented  and  renamed  for  ease  of  programming. 

The  paragraphs  in  Appendix  I  have  been  numbered  to  facilitate  cross- 
referencing.  37 


2.0  GLOSSARY  OF  THE  SUBROUTINES 


The  program  operates  via  a  MAIN  rout^ae  and  15  auxiliary  routines 
which  are  briefly  described  and  listed  below.  Standard  I/O  and 
mathematical  routines,  e.g.,  SIN,  LOG,  etc.  are  assumed  to  be  available 
through  the  FORTRAN  operating  system. 

2.1  The  MAIN  routine  controls  overall  run  processing  and 
computes  the  I,  J,  K  and  L  matrices. 

Routines  called  are: 

RDDAIA 

GENLGP 

GENKR 

GENBSL 

PRTMTX 

PRCSSM 

2.2  Subroutine  RDDAIA  reads  the  user's  control  parameters 
and  sets  up  preliminary  output  heading  information. 

Routine  called  is:  CALENP 

2.3  Subroutine  CALENP  computes  the  sections  of  9,  the  polar 
angle,  and  the  step  size  for  numerical  integration. 

2.4  Subroutine  GENLGP  computes  the  associated  Legendre 
functions  over  the  necessary  range. 

2.5  Subroutine  GENBSL  controls  backward  recursion  of  Bessel 
functions  and  forward  recursion  of  Neumann  functions. 

2.6  Subroutine  BESSEL  computes  the  Bessel  function  for  a 
specific  argument  and  order. 

2.7  Subroutine  GENKR  computes  the  parameter  "kr"  and  its 

derivative  with  respect  to  the  polar  angle  6. 
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2.8 


Subroutine  PRTMTX  prints  the  headings  and  controls  the 


printout  of  the  I,  J,  K  and  L  complex  matrices. 

Routine  called  is:  PR1NTM 

2.9  Subroutine  PRINTM  prints  the  elements  of  a  specified 
matrix  of  specified  rank. 

2.10  Subroutine  PRCSSM  generates  the  Q  matrices  from  the 
I,  J,  K  and  L  matrices,  and  transforms  the  Q  matrix  into  the  T 
matrix. 

Routines  called  are:  NRM3MX 

CNDTNQ 

PRTRIT 

ADDPRC 

2.11  Subroutine  NRM3MX  normalizes  the  I,  J,  K  and  L  matrices 
to  obtain  the  Q  matrices. 

2.12  Subroutine  CNDTNQ  conditions  the  Q  matrix  before  trans¬ 
forming  it  into  the  T  matrix. 

2.13  Subroutine  PRTRIT  prints  headings  and  controls  printout 
of  the  T  matrix. 

Routine  used  is:  PRINTM 

2.14  Subroutine  ADDPRC  does  final  processing  of  the  T  matrix 
to  provide  the  scattering  results. 

Routines  called  are:  GENLGP 

VECMUL 
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2.15  Subroutine  VECMUL  multiplies  a  matrix  times  a  vector. 

2.16  Subroutine  DUMP  gives  a  listing  of  core  storage  when  an 
error  condition  occurs . 

Subsequent  paragraphs  detail  the  above  routines  where  necessary. 

It  should  be  noted  at  this  point  that  standard  mathematical  notation 
is  not  necessarily  followed,  e.g.,  program  notation  labels  Bessel 
functions  as  B  instead  of  j.  This  was  done  for  ease  of  relating 
program  mnemonics  to  mathematical  notation.  When  necessary,  parameters 
have  been  labeled  which  have  notation  different  from  the  earlier  text. 

3.0  THE  INPUT  ROUTINE 

Subroutine  RDDATA  reads  the  user's  control  information,  prints 
out  headings  and  obtains  information  for  numerical  integration.  The 
input  cards  and  their  formats  are  listed  below. 

3.1  Card  1  NM,  NRANK,  NSECT,  IBODY ,  NUANG 

FORMAT  (5112) 

NM  No.  of  values  of  "m" .  See  Card  3. 

NRANK  Rank  of  matrices  I,  J,  K  and  L 

NSECT  No.  of  sections  defining  body  shape  and 

integration  intervals.  See  Subroutine 
CALENP  for  fuller  description  of  body 
shapes . 

IBODY  Case  No.  or  body  shape  identifier 

* 

7  :  Spheroid 

■  8  :  Mirror  Symmetry 

9  :  General  Axisymmetric  Case 

NUANG  No.  of  aspect  angles  "u" .  See  Card  5. 
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3.2  Card  2. 


CONK 

BRXT 

ALPHA 


CONK,  BRXT,  ALPHA 
FORMAT  (3E12.7) 

ka,  scale  factor  for  r,  the  polar  radius, 
In  determining  body  shape, 
variable  parameter  to  be  used  in  computing 
body  shapes. 

(X,  or  a  variable  parameter,  to  be  used  in 
calculating  body  shapes.  For  a  fuller 
description  of  its  usage  see  Subroutines 
CALENP  and  GENKR  described  below. 


3.3  Card(s)  3.  CMI(I),  I  -  1,  NM 

FORMAT  (6E12.7) 

CMI(I)  1^  value  of  "m"  to  be  used  in  current 

solution  of  scattering  problem.  As  many 
as  30  values  of  "m",  the  azimuthal  index, 
may  be  read  in;  "m"  is  any  integer  ^  0. 


3.4  Card  4.  NDPS(I),  1*1,  NSECT 

FORMAT  (6112) 

NDPS(I)  No.  of  divisions  for  integration  in  I*"*1 

section  of  the  body  shape.  The  body  may 
be  divided  into  as  many  as  6  sections. 
These  parameters  are  used  to  calculate 
spacing  for  numerical  integration,  and 
they  must  be  a  multiple  of  4. 


3.5  Card(s)  5. 

UANG(I) 


UANG(I) ,  I  -  1,  NUANG 

FORMAT  (6112) 
til 

I  value  of  "u",  a  member  of  a  table  of 
aspect  angles  (in  degrees).  As  many  as  60 
values  of  "u"  may  be  read  in. 
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4.0 


CALCULATION  OF  END  POINTS  AND  SPACING  FOR  INTEGRATION 


Subroutine  CALENP  is  one  of  two  special  routines  that  have  to 
be  written  into  the  program  for  specific  body  shapes.  This  routine 
calculates  "NSECT"  values  of  the  polar  angle  6,  which  provide  boundaries 
for  dividing  the  body  into  sections  for  numerical  integration.  With 
each  boundary  point  a  value  of  3  is  associated.  The  spacing  for 
integration  is  then  determined  by  dividing  the  range  of  3  by  the 
correct  value  of  "NDPS".  Note  that  the  number  of  divisions  does  not 
have  to  remain  constant  from  one  section  to  the  next,  but  it  must  be 
a  multiple  of  4. 

Since  the  computations  for  each  special  version  of  CALENP  may 
vary,  the  following  parameters,  "ALPHA" ,  'feRXT",  '1)3",  "SNALPH",  and 
"CSALPH"  may  be  used  for  communicating  between  routines  special 
values  associated  with  a  particular  body  shape.  Note  the  use  of  QB 
below  for  a  variable  peculiar  to  the  sphere-  ne-sphere  shape. 

In  Appendix  II,  a  listing  of  the  routines  CALENP  and  GENKR 
are  given  for  a  sphere-cone- sphere  body  (Figure  2.). 

For  the  sphere-cone-sphere  body,  three  end  points  for  8; 

8 ^ ,  8^,  8^,  are  computed  as  follows: 


6 


1 


tan 


sinfl  cosCL  1 
L  q  -  sin^CX  J  ’ 


0  *  0  ^  105° 


NOTE:  Cl  (Figure  2)  is  an  input  parameter  stored  at  "ALPHA" 
b/a  (Figure  2)  is  an  input  parameter  stored  at  "BRXT" 
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q 


(1  -  b/a)  (1  -  sirQ.) 
2 


and  is  computed  and  stored  in  QB. 


tan'1  L 


s  ird  cos  a. 


l-q-(b/a)cos  a 


]* 


V 


0„  <  n 

2 


5.0  THE  FIRST  CONTROL  ROUTINE 

The  MAIN  Routine  controls  the  general  flow  of  the  program  and 
computes  the  real  and  imaginary  parts  of  the  partitioned  scattering 
matrices,  I,  J,  K  and  L.  After  the  user's  control  data  has  been 

r(J 

read  in,  a  numerical  integration  system  utilizing  Bode's  3  order 
(14) 

rule  is  utilized.  The  program  confutes  the  I  and  J  complex 

matrix  elements  for  the  three  cases,  axisymmetric ,  mirror -symmetric 
and  spheroidal  as  follows: 

5.1 

n 

Itj  *  m  J  (sin6)  (kr)2  B^(kr)  H^kr)  {(i  +  j)  cosQ  P®  P® 

0 

-  (i  +  m)  P®  P®_x  -  (j  +  m)  P®  P®_J  d0 


Bj (kr) 


Bessel  function  of  the  first  kind  of  order  "j"  and 
argument  "kr" 
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H^(kr)  :  Hankel  functions  which  are  defined  as 


Bi(kr)  +  /-1  N  t  (  kr  ) 

N^(kr)  :  Neumann  functions  of  order  "i"  and  argument  "kr". 

P™  (cos6) 

P"1  -  —  7T  — 

i  sinQ 


where 


P™  (cos6) 

is  the  associated  Legendre  function,  of  rank  m,  and  ord*r  i. 
r  :  Polar  radius  used  in  calculating  body  shape, 
i  :  Subscript  notation  for  iC  row  of  the  matrix, 

j  :  Subscript  notation  for  jC  column  of  the  matrix. 

5.2 

TT 

Jj,  **  ~  j  (sin0)  H^(kr)  |p™  Pm  |^kr(kr  B  ^(kr)- jB ,  (kr)  (m^+  i  jcos^O) 

J  em  J  J  J 

0 

+  i  j(j  +  1)  Bj(kr)sin0  cos6j 

-(i  +  nOP^P™  £kr  cos0(kr  Bj^krJ-j  B  ^  (kr)+(  j+l)^|p^B  ^  (kr)sinBj 
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where  : 


+(j  +  krj^kr  B^^krJ-j  Bj(kr)^)^(i  +  nOP^j-i  cosO  p™ 


m 


c 


for  m  =  0 
for  m  s  1 


Within  the  program  each  element  of  the  I  and  J  arrays  is 
used  as  an  accumulator  for  numerical  integration  under  Bode's  rule. 
Thus,  for  a  specified  value  of  0,  all  necessary  functions  are 
computed  and  added  to  the  correct  matrix  element. 

To  save  computer  time,  computations  which  would  produce  a 
null  contribution  to  the  integration  are  eliminated,  and  the  following 
symmetries  are  taken  advantage  of  in  the  direct  computation  of  the 
I  and  J  matrices. 

5.3  General  axisymmetric  bodies: 

Re(Iji>  =  Re<Iij> 

Re  ( J  j i >  =  Re<Jij> 

5.4  Mirror- symmetric  bodies:  use  paragraph  5.3  plus 

{both  odd 
both 


J^,j  =  0;  if  i  and  j  are 


f  odd, 
V,even, 


even 

even 

odd 
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5.5  Spheroids:  use  paragraphs  5.3,  5.4,  plus 


In'(I  j  i>  =  Im(Iij^ 

Im(J j i>  = 

The  K  and  L  matrices  are  then  calculated  from  the 
following  relationships  with  the  I  and  J  matrices: 


5.6  Re(Kij)  “  “Re(J  ) 


Im(K . .)  =  -Im(K  )  +  b  D, , 

ij  ij  ij 

where 


b 


{ 


1.0  General  Axisynmetric  bodies 

0.5  Mirror-Symmetric  or  Spheroidal  bodies 


0.0;  i  4  j 


ij 


[ 


e  (2i  +  1)  (i  -  m)!vl 


m 


4i( i  +1)  (i  +  m) 


3 


i  =  j 


5.7  Re(Lj.j)  =  -Red^) 

Im(Ltj)  =  -Im(I_) 

The  I,  J,  K  and  L  complex  matrices  are  then  printed  by 
Subroutine  PRTMTX  and  control  passes  to  Subroutine  PRCSSM  for  further 
processing . 


6.0  ASSOCIATED  LEGENDRE  FUNCTIONS 

Subroutine  GENLGP  generates  the  associated  Legendre  functions, 
P™(x)  for  a  given  argument  x,  a  given  value  of  the  azimuthal  index  m, 
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and  for  all  values  of  degree  i  from  0  to  "NRANK",  the  Input- 
specified  rank  of  the  matrix.  The  first  two  values  of  P  are  gener¬ 
ated  by  formula,  then  the  remaining  values  of  P  are  generated  by 
a  recursion  relationship. 

The  following  formulae  are  used  to  generate  P™(x).  Note  that 
for  this  particular  program,  the  functions  always  appear  in  the  context 


P™(cos6) 

sin8 


P™(cos0) 

sinB 


0.0 


i  <  m 


Pm(cos6) 

m 

s  in9 


(2»>;  ■n*'1(8) 


m! 


m 


P0(cos6)  =  1.0 

sin©  sinB 


i  =  m  =  0 


P^(cos0) 

sin0 


COS0 

sin@ 


Recursion  relationship: 


Pm(cos6) 

n _ 

sin0 


p”  ,  (cos0)-j 

(2n-l)cos0  l  ~  ~s7ne - J-<n  +  m  - 


P™  9(cos0) 
n- z 

sin© 


n  -  m 
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7.0 


BESSEL  FUNCTIONS 


Subroutine  BESSEL  generates  a  Bessel  function  of  the  first  kind 

Bn(x),  for  a  specified  argument  x,  and  order  n,  by  means  of  an 

infinite  series.  To  preserve  accuracy,  the  computations  are  performed 

in  double  precision  arithmetic  and  truncation  error  due  to  neglected 

-20 

terms  in  the  series  is  <  10  .  If  the  series  has  not  converged  to 

the  aforementioned  accuracy  before  the  computation  of  the  100^  term, 
an  error  indication  is  given. 

The  following  infinite  series  is  used  to  compute  a  Bessel 
function: 

00 

Bn(x)  "  1-3-5-* -(2n  +  1)  I  ai 

i-0 

where : 


2 

-x 

ai  "  2il2n  +  (21  +  1)J  ®i-l 

8.0  RECURSION  REIATIONSHIPS  F  <R  BESSEL  AND  NEUMANN  FUNCTIONS 
Subroutine  GENBSL  calls  Subroutine  BESSEL  to  obtain  two 
successive  BESSEL  functions  for  a  specified  argument,  and  then  uses 
these  first  two  values  to  recurse  backward  over  the  range  of  i  from 
NRANK  to  0.  If  the  two  computed  functions  of  order  NRANK  and  NRANK-1 
do  not  satisfy  the  accuracy  requirements  mentioned  in  paragraph  7.0, 
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the  routine  will  increase  the  order  of  the  computed  BESSEL  function 
to  4(NRANK).  If  this  fails  to  produce  a  satisfactory  pair  of 
functions,  the  run  will  abort  and  a  dump  of  core  memory  is  taken. 

The  recursion  relation  used  for  computing  BESSEL  functions  is: 

Bn_1(x)  *  (2n  +  ljx  1  Bn(x)  -  Bn+1(x)  • 

This  routine  also  computes  Neumann  functions  by  a  forward 
recursion  formula  after  the  first  two  values  are  computed  by  the 
following  formulae: 


n0OO 


-cosx 


Nl(x) 


-cosx 


sinx 

x 


The  recursion  relation  used  for  computing  the  remaining  Neumann 
functions  is: 


N+1(x)  =  (2n+l)x  1  N  (x)  -  N  .(x)  . 

n+j.  n  n-  i. 


To  test  the  accuracy  of  the  functions  over  the  range  of  computed 
Bessel  and  Neumann  functions  for  a  given  argument,  two  tests  are 
performed  in  the  MAIN  Routine  after  the  vector  of  functions  from  0 
to  NRANK  is  computed.  If  the  following  relations  are  not  satisfied 
to  an  accuracy  of  10  an  error  message  indicating  such  a  condition 
is  printed,  and  the  program  continues.  Though  the  tests  are  performed 
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In  the  MAIN  Routine  after  the  call  to  Subroutine  GENBSL,  for  convenience 
they  are  listed  here: 

Bessel  Test: 

jx2  |_B1(x)N0(x)  -  Bq^N^x)]  -  I  !  <  10"10  . 

Neumann  Test: 

i  2f  1  ,  _  1  n 

Lbnrank(x)nnrank-i(x)  "  bnrank-i(x)nnrank(x)J  "  1  I  <  10 

9.0  GENERATING  the  body  shape 

Subroutine  GENKR  is  one  of  two  custom  written  routines  which 
are  adapted  to  the  particular  body  shape  in  question.  As  noted  above 
in  Subroutine  CALENP  certain  parameters  are  available  to  the  programmer 
to  use  as  he  sees  fit  to  communicate  information  from  one  routine  to 
another.  The  basic  function  of  all  versions  of  GENKR  is  to  compute 
the  polar  radius  r  as  a  function  of  the  polar  angle  6,  to  compute 
^  and  to  scale  these  values  by  the  input  constant  ka  ■  CONK. 

To  illustrate  the  use  of  this  routine,  a  sphere-cone-sphere 
body  shape  is  used  (Figure  2).  As  a  result  of  subroutine  CALENP  the 
major  divisions  of  the  body  as  a  function  of  9  have  been  recorded. 

This  routine,  given  a  value  of  6  now  computes  (ka)r  and  ka(dr/d0)  ; 
the  scale  factor  ka  is  an  input  to  the  program. 

9.1  Section  1  0  *  9  *  0 

_  cos9  f  _  (  fl_.sin9  Y1.* 
sir CL  L  >,  sirfi  s  J 


i 


dr 

d8 


-g  sln8 
sinO. 


sin8  co>.'8 


q  stn8 
sina 


NOTE:  q  was  computed  in  Subroutine  CALENP  -md  stored  in  location  QB. 
9.2  Section  2  0  <  0  £  @2 

_  1  -  q 

r  sin(8-a) 


dr 

d8 


-(1-q)  cos(8-g) 
sin2(8-a) 


9.3  Section  3  02  <  8  £  tt 


dr 

d8 


.  1’"sina^"q]sln6  '  C  1"’sCi^h'q)28ln8  cose[(b/a)2  -  (  sin^ 


10.0  FIRST  MATRIX  PRINTOUT 

Subroutine  PRTMTX  controls  the  printout  of  the  I,  J,  K  and  L 
matrices.  Both  the  real  and  imaginary  arrays  comprising  each  of  these 
matrices  are  labeled  and  printed  out  on  the  community  output  tape. 

This  output,  which  was  originally  intended  as  an  intermediate  printout 
for  checking  the  program,  may  be  eliminated  by  removing  the  "CALL  PRTMTX" 
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statement  which  follows  Fortran  statement  860  in  the  MAIN  Routine. 
11.0  PRINTOUT  OF  AN  ARRAY 

Subroutine  PRINTM  will  print  out  a  specified  square  array  of 
given  rank. 

12.0  GENERATING  THE  Q  MATRIX  AND  THE  T  MATRIX 

Subroutine  PRCSSM  is  the  second  major  control  routine  and  it 
controls  the  transformation  of  the  I,  J,  K  and  L  matrices  to  the  "Q" 
matrices,  and  the  subsequent  solution  of  a  matrix  equation  which 
provides  the  "T"  matrix. 

Subroutine  NRMQMX  (see  below)  normalizes  the  I,  J,  K  and  L 
matrices  to  produce  the  Q  matrix. 

For  notational  convenience  we  defines 

Q  -  Re(Q)  +  i  Im(Q)  -  (  £l  ) 

^2  ^4 

where  ___ 

i  *  / -l  . 

The  method  currently  used  by  the  program  to  transform  the  Q 
matrix  into  the  T  matrix  involves  orthogonal izing  the  Q  matrices. 
After  the  complex  Q  matrix  has  been  generated  by  normalizing  the 
I,  J,  K  and  L  matrices  it  is  in  the  form  noted  in  paragraph  12.  From 
these  Q  matrices,  a  new  complex  Q  matrix  of  rank  2N  is  generate 
from  the  following  relations: 
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A 

Q(2m-1)  (2n-l)  =  (Vm  n 

^(2m-l)  (2n)  *  (Q2)m  n 

m,n  M  1,  2 . N 

A 

Q(2m)  (2n-l)  =  ®3'm  n 

A 

^(2m)  (2n)  '^4  m  n 

The  new  Q  matrix  is  next  conditioned  as  outlined  in  Subroutine 
CNDTNQ  of  paragraph  14.0  below. 

Orthogonalization  then  proceeds  as  follows. 

1)  Consider  each  row  of  Q  as  a  vector  with  ZN  components; 
e.g.  the  components  of  the  first  vector  would  be: 
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<c  i 


Q  * 


2N 


r-1 


qr  *  ^pl^ql  +  ^p2^q2  + . ^p(2N)^q(2N) 


3)  Orthognnalize  fi(2M.1)  to  $(2N)  : 


r  e-  * 


^2N-1  “  ^2N-1  ’  L^2N  ‘  ^2N-lJ& 


2N 


A)  Normalize  ,  : 

ZN-  1 


^2N- 


^2N-1 


1  ICl-SsH-lP 


5)  Orthogonal ize  ^2N_2  to  both  ^  and  ^ 


2N-2 


^2N-2  ’  L^2N-1  ‘  ^2N-2^2N-1  "  L^2N  *  ^2N-2]^2N 


6)  Normalize  fi2N_2  : 


^2N-2 


^2N-2 


pf  “ 

1&N-2  '  ^2N-2J 


7)  Continue  t.he  orthogonal izat ion  and  normalization  process 
until  2]^  has  been  orthogonalized  to  all  subsequent  rows. 
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8)  A  complex  matrix  *  Is  now  generated  from  the  complex 
matrix  3  by  the  following  relation: 


T  -  Q**Re  (Q) 


9)  The  T  matrix  is  then  decomposed  Into  the  matrices 

T  T  and  T.  by  the  reverse  of  the  procedure  In 
2*  3  4 

paragraph  12.1. 

The  complex  T  matrix  Is  printed  by  Subroutine  PRTRIT  and  then  the 
final  processing  Is  performed  by  Subroutine  ADDPRC. 

13.0  NORMALIZING  MATRICES 

The  Subroutine  NRM3MX  lormalizes  the  I,  J,  K  and  L  matrices 
to  obtain  the  Q  matrix.  The  Q  matrix  Is  blocked  as  noted  above 
In  paragraph  12.0  and  the  following  procedure  Is  used: 


-  (Z2) J'  (Z2“^)  =  Q 


-1L 


l*  /zzM 


/ Z2 ” 


q2  -  -(Z2)"%  L'  (Z2)"% 

q3  -  (Z2)'^  I'  (Z2)"* 

-  <Z2)'*  K'  <Z2)"* 
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whe  re : 


e  (2n+l)  (n-ra) !  . 

Z2n  "  L  4n(n+l)  (n+m) !  J 

The  expression  for  Z2  is  Che  same  as  that  used  in  computing  Che 

n 

K  matrix  of  paragraph  5.6.  The  prime  on  J,  etc.  denotes  matrix 
transpose  as  seen  from  the  second  half  of  the  equality  statement. 

14.0  CONDITIONING  MATRICES 

After  the  matrix  Q  of  rank  2N  has  been  formed,  the  matrix 
is  conditioned  starting  with  the  last  row  and  working  towards 

row  2^  • 

14.1 

C&2N^t  ~  [_ 1  / Im(^2N ) 2n1  5  1  “  2»  •••  2N 

The  notation  (2.2N^i  re^ers  to  the  i^  element  of 
the  (2N)t^1  (last)  row  vector.  Now  set 

^rn  ■  L  ■ 

where  the  equivalence  is  performed  for  each  of  the  2N  elements  of 
Q  ,  and  repeated  for  all  rows  m  *  1,  2,  ...,  2N-1. 

14.2  Redefine 

^2N-1  =  L1/Im<^2N-1^2N-l]^2N-l  * 
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Chen  compute 


^m  "  ^m  "  LIm^m*2N-lJ£2N-l 

for  all  rows  m  *  1,  2,  ...  2N-2. 

14.3  Continue  the  process  of  paragraphs  14.1  and  14.2  for 
all  the  remaining  rows.  The  final  step  in  the  process  is  to  generate 


$2  '  Li/i»<42)2]22  • 

2X  -  2X  -  li"<£i>2>2  • 


14.4  Set 

Im(Q^) ^  *0.0;  i  *  m  +  1,  m+2,  2N;  m  *  1,  2,  3,  ...,  2N-1  . 

15.0  PRINTING  THE  T  MATRIX 

Subroutine  PRTRIT  controls  the  printout  of  the  T  matrix, 
both  real  and  imaginary  elements,  in  the  same  manner  as  Subroutine 
PRTMTX  (paragraph  11.0  above)  controls  the  printout  of  the  I,  J,  K 
and  L  matrices.  The  community  output  tape  is  used.  Since  this 
printout  is  used  mainly  for  checkout,  it  can  be  eliminated  by  removing 
the  "CALL  PRTRIT"  statement  following  FORTRAN  statement  140  in 
Subroutine  PRCSSM. 
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16.0  FINAL  CONTROL  ROUTINE 


Subroutine  ADDPRC  is  the  third  and  last  control  routine  which 
converts  the  T  matrix  to  the  final  set  of  results.  Two  sets  of 
results  are  generated,  a  set  of  answers  for  the  current  value  of  m 
and  an  accumulated  set  of  answers  for  all  values  of  m  up  to  and 
including  the  present  value  of  m. 

To  generate  the  final  results,  the  following  procedure  is 
followed:  The  T  matrix  is  normalized 

T(l.)  -  (Z2i>^  Ttj(k)  (Z2  >'*  block 

k  indicates  1  of  4  blocks; 

T  ’  ^ 

Z2n  is  as  defined  in  Subroutine  NRhQMX  under  paragraph  13.0. 

NOTE:  For  the  reader  who  is  relating  the  mathematics  to  the  program 
listing  in  Appendix  II,  the  mapping  of  COMMON  storage  in  paragraph  19.0 
should  be  consulted. 

The  associated  Legendre  functions  of  form 

Pra(cos  u) 

n _ _ 

sin  u 


where : 
NOTE: 
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are  generated  for  each  value  of  the  aspect  angle  u,  and  for  n  ■  1 
to  NRANK.  The  derivatives  of  the  Legendre  functions  are  computed  from 


<l[p”(co»  u)]  r  *“<“■  u>  1  r  Cl<cos  u)  1 

- dJ -  ■  n  cot  u  L  ".in  u  ~  J  -  fr  +  "»  L - ,1„  u  J  • 


11  2  2 
Values  of  the  vectors  F  ,  G  and  F  ,  G 

Subroutine  VECMUL.  These  vectors  are  defined  as: 


are  generated  by 


16.1 


r(l)n  m  Pm(cos  u)/sln  u 
n 


^-(l)n+1  d[p“(cos  u)^/di 


and 

16.2 


n+^  d^P™(cos  u)J/du\ 

n  m  Pm(cos  u)/sln  u / 
n  ' 


The  final  sets  of  answers  are  generated  from  the  following 
equations : 


SCATT  1,  2 


16 

(ka)2 


NRANK 

l 

n*l 


wuuv 

I  <z2>'1Df„1,2i2 
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* 


TOTAL  1,2 


16 

<ka)2 


NRANK 

£  (Z2)~L(-t)\Fn 


1,2 


ra  Pm(cos  u) 
n _ _ 

sin  u 


n*l 


RTRAD  1,2 


8 

(ka) 


NRANK 

£  <Z2)_1(i)n[Fn1,2 
n*l 


m  Pm(cos  u) 
n 

sin  u 


The  final  results  are  divided  into  two  classes  as  noted  above 
by  the  quantities  SCATTl,2etc.  The  classes  are  two  different  incident 
polarizations.  Class  1  is  the  E-parallel  incidence;  class  2  is  the 
E-perpendicular  incidence. 


Appendix  III  contains  a  listing  of  a  sample  output  of  the 
sphere-cone-sphere-body  shape.  The  printout  titles  and  their  meanings 
are : 


ANGLE 
SCATT  1,  2 

TOTAL  1,  2 

RTRAD  1,  2 

RCS  1,  2 


u,  the  aspect  angle  (degrees) 

Scattering  cross-section  for  each 

2  -I 

class,  normalized  by  [na  J 


Complex  forward  amplitude,  normalized 
by  [TTa2]"* 

Complex  back  scattered  amplitude 

r  2--h 

normalized  by  trra  J 


Radar  cross  section  normalized  by 


[TTa2]"1 


NOTE:  RCS  =  | RTRAD |  only  is 
computed  for  the  accumulative  case . 
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PHASE  ANGLE  1  .  2 


For  the  accumulative  case,  a  phase 
angle  is  computed: 


PHANO  -  TAN'1  [  ] 

This  is  the  phase  angle  of  the  back 
scattered  amplitude  (Degrees). 

l/.o  MULTIPLYING  A  MATRIX  TIMES  A  VECTOR 

Subroutine  VECMUL  is  one  of  two  routines  coded  in  both  machine 

language  and  FORTRAN.  It  multiplies  a  matrix  times  a  vector  to  compute 

~T  I  ~ 2  2 

the  vectors  F  ,  G  and  F  ,  G  of  paragraphs  16.1  and  16.2.  The 
r.vichine  coded  version  has  the  advantages  of  higher  speed  and  accuracy. 

18.0  CORE  DUMP 

If  an  abnormal  or  uncorrectable  error  condition  occurs, 
Subroutine  DUMP  gives  a  dump  of  core  memory  as  an  aid  in  debugging 
the  error  condition.  The  Subroutine  LBPDMP  is  a  system  routine  for 
dumping  core  between  specified  limits. 

19.0  STORAGE  ARRANGEMENTS 

To  conserve  and  fully  utilize  core  storage,  three  large  matrix 
arrays  of  dimension  120  x  120  have  been  set  up  in  an  area  of  COMMON 
storage  named  "MTXCOM".  To  aid  in  programming,  various  routines  use 
EQUIVALENCE  statements  to  resegment  these  large  arrays  into  manageable 
blocks . 
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The  FORTRAN  array  names  of  the  three  major  blocks  are: 

CMTXRL  (120,  120) 

CMTXIM  (120,  120) 

SPRMTX  (120,  120) 


Within  the  MAIN  Routine  the  following  overlays  are  made: 

r 


CMTXRL 


CMTXIM 


AMXIR  (60,  60) 

RE(I) 

AMXJR  (60,  60) 

RE(J) 

AMXKR  (60,  60) 

RE(K) 

^AMXLR  (60,  60) 

RE(L) 

r 

AMXli  (60,  60) 

IM(I) 

AMXJI  (60,  60) 

IM(J) 

AMXKI  (60,  60) 

IM(K) 

w AMXLI  (60,  60) 

IM(L) 

The  SPRMTX  block  is  unused. 

Subroutine  NRMQMX  and  Subroutine  PRCSSM,  the  second  control 
routine,  allocated  storage  as  follows: 


CMTXRL 


SPRMTX 


QMTXII 

(60,  60) 

QI1 

(60, 

60) 

IMO^) 

QMTXJI 

(60,  60) 

QI2 

(60, 

60) 

im(q2) 

QMTXKI 

(60,  60) 

QI3 

(60, 

60) 

im(q3) 

^QMTXLI 

(60,  60) 

QI4 

(60, 

60) 

im(q4) 

^QMTXIR 

(60,  60) 

QE1 

(60, 

60) 

re(q1) 

QMTXJR 

(60,  60) 

QR2 

(60, 

60) 

RE(q2) 

QMTXKR 

(60,  60) 

QR3 

(60, 

60) 

re(q3) 

^QMTXLR 

(60,  60) 

QR4 

(60, 

60) 

re(q4) 
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NOTE:  After  Subroutine  NRM3MX  normalizes  and  moves  the  Q  matrix 


(complex)  into  the  SPRMTX  and  CMTXRL  areas,  the  processing  which  trans¬ 
forms  the  Q  to  the  T  mitrix  fcllowr  the  procedure  outlined  in 
paragraph  12.0.  The  storage  allocation  is  noted  above  in  the  eight 
itemized  steps. 

Subroutine  INVMBL  always  assumes  the  block  matrix  which  is  to 

be  processed  is  stored  in  the  CMTXRL  area.  The  intermediate  steps 

as  outlined  in  paragraph  14.0  are  performed  in  the  CMTXIM  area. 

The  third  and  last  control  routine,  Subroutine  ADDPRC  makes 

the  following  storage  allocations: 

r 


CMTXRL 


QMTXII  (60,  60) 
QMTXJI  (60,  60) 
QMTXKI  (60,  60) 
QMTXLI  (60,  60) 


FGVECT  (2,  120,  2) 
FGMUL  (120,  2) 
FGANS  (60,  10) 
(unused) 


CMTXIM 

SPRMTX 


TCMPLX  (2,  120,  120)  :  Real  and  imaginary 

components  of  the 
T  ratrix. 


NOTE:  That  TCMPLX  overlays  both  the  CMTXIM  and  SPRMTX  areas.  As 

~ 1  i  ~ 2  2 

noted  in  paragraph  16.0,  FGVECT  contains  the  F  ,  G  and  G  ,  F 
vectors.  The  first  subscript  refers  to  the  real  and  imaginary 
components  of  the  vectors,  the  second  subscript  refers  to  the  dimension 
of  the  vectors  which  is  2*NRANK  and  the  last  subscript  differentiates 
the  2  vectors. 

FGMUL  contains  the  vectors  which  post-multiply  the  T  matrix 
to  generate  FGVECT .  The  first  end  second  subscripts  correspond  to 
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the  second  and  third  subscripts  of  FGVECT . 


FGANS  contains  the  final  answers.  The  first  subscript  corres¬ 
ponds  with  the  value  of  aspect  angle  which  generated  it  and  the  second 
subscript  refers  to  the  answers  in  the  following  manner: 


1 

2 

3 

4 

5 

SCAT!1 

Re (TOTAL1) 

ImCTOTAL1) 

ReCRTRAD1) 

ImCRTRAD1) 

6 

7 

8 

9 

10 

SCATT2 

Re (TOTAL2) 

Im (TOTAL2) 

Re(RTRAD2) 

Im(RTRAD2) 
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APPENDIX  II 

THE  FORTRAN  IV  PROGRAM  LISTING 
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T  SU8T VP  E»  FORTRAN, LMAP ,L  STRAP 

C  SCATTERING  FROM  AXISVHHETRIC  CONDUCTORS  FOR  CASFS  7,  8  AND  9. 

COMMON  OTR.RTO.CPI 

COMMON  /CMVCOM/  NM,  CM  1 1  30  I  *CMV,  X MV ,  CM2  .  F M,  QR  m,  T WM ,  PRODM 
COMMON  /FNCCOM/  PNMll Gl 611  * 8 SSL SPI  6 II . C NF UMN I  6 1 » 

COMMON  /THTCOM/  T HE  TA , NTH ET A « DL THT A , S I  NTH ,C  (IS TH , I SMRL , I SWTCH I  7  I ,SR 
1MUL*  SMULSSI7) ,COH < 6 > , OHM, NSFC T , NDP S 1 6 ) ,FPPSI6) ,KSRCT 
COMMON  /MTXCOM/  NRANK  *  NR ANK 1 , AMX I R ( 60,60) ,  A  MX  JR  I  6  3 , 69  ) , AMXKR (SO*  6r 
1  ) , AMXLRI 60 *60  I  *  AMX 1 1 I  60,  601  ,  AMX  JM  60, 60  >  ,  AMXK  1 1  63 ,60 )  , AMXll 160,60  » 
2, SPRMTXI 1 20, 1 20) ,CMXNRM(60> 

COMMON  /BCVCOM/  CKR  ,  OCX  R  ,  CKR  2  »C  S*  R  X  ,  SNK  R  X  ,C  ONK  ,  RR  X  T  ,  Al  PH  A  ,  ( R0DV,0« 
l.SNALPH.CSALPH 
DIMENSION  CLRMTXI 402001 
EQUIVALENCE  I AMX IR.CLRMTX ) 

C  SET  PROGRAM  CONSTANTS. 

OTR  *  1.745329251 9943F- 02 

RTO  »  57.2957795131 

CPI  «  3.1415926535898 

I  SWT  CHI  1)  «  2 

ISWTCHI 2 1  '  3 

ISWTCHI 3)  ■  4 

ISWTCHI4)  «  1 

SMULSSI1)  «  32.0 

SMULSSI2)  -  12.0 

SMULSSI3)  «  32.0 

SMULSSI 4)  -  14.0 

C  call  ROUTINE  TO  KFAD  DATA  AND  PRINT  HF ADI  NOS  f OR  OUTPUT 

20  call  rdoata 

IF  I  I  BODY- 9  1 24, 22, 24 
22  BOVTCT  «  1.0 
GO  TO  26 
24  BDVFCT  ■  0.5 

C  SET  UP  A  LOOP  FOR  M  ANn  SFT  VARIAmIES  WHICH  ARH  A  FUNCTION  OF  M. 

26  00  900  IM  »  1 , NM 
CMV  «  CM  I  I  IM) 

XMV  «  CMV 
CM2  «  CMV4CMV 
PRODM  •  1.0 
I  F  I  CMV 1 40 , 43  ,  44 
40  EM  «  1.0 
GO  TO  60 
44  EM  •  2,0 

QUANM  *  CMV 
00  52  IFCT  «  1, XMV 
QUANM  >  0UANMC1.0 
PROOM  «  QUANM«PR0DM/2.0 
52  CONTINUE 
60  OEM  •  -2.0/EM 
TWM  «  CMVCCMV 

C  INITIAL IZF  ALL  MATRIX  ARF AS  TO  7FR0 

HO  80  I  »  1,28800 
CLPMTXI I  I  «  0.0 
80  CONTINUE 

C  SFT  UP  A  LOOP  FOR  ALL  VALUES  OF  ThFTA. 

THFTA  «  0.0 

C  SET  UP  GENERAL  LOOP  FOR  CORRFCT  NUMBFR  OF  INTEGRATION  SFCTIONS. 

DO  800  ISECT  •  l.NSFCT 
XSECT  »  ISECT 
NTHETA  »  NOPSI ISECT  >C1 
OLTHTA  «  COHI  ISECT! 
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n mm  ;  nirnfA/p^.s 

i  s  Mm  =  <, 

OU  7  JO  !  T  m  T  ft  =  l.MHFTft 

C  SFT  SWIIC.HIS  AM)  Mill  T  |  M|  I  rjs  FC-i.  SIMPSONS  !  HI  f  t.u  A  T I  UN  MFTHOO. 

!F(  ITMTA-  1’  ,  H  > 

1 2 J  SRMUl  -  7.C»Omm 

IM  ISFCT-  117'.i,  7T,  14? 

132  IM  ITHTA-NT'-F  TA  I  ?GO,  148, 1  4>» 

148  SRMUl  =  T.C*Oh< 

CO  TO  340 

200  ISMRl  x  I SWTCh<  ISMRt  ) 

SRMUL  -  Ml>L  SSI  I SMR  l ) *OHM 
340  ThFTA  =  THf TACOl  TMTA 
348  COST  H  =  CnS<  THF  TA ) 

SINTH  *  SIM  THFTAI 

f  GENERATE  THF  l  FC.fnORF  POL YNOMI At S. 

CAll  G^NIGP 

C  EVAUJATF  KP  AS  A  FUNCTION  OF  THFTA.  At  Srl  ITS  OFPIVATIVF. 

C At  L  GFNKR 

CSKRX  *  COSICKW )/fKB 
SNKRX  *  S  INI CKR  I /CKP 
CK  R  2  *  CKR*CKR 

C  GENERATE  PFSSFL  FUNC T ION S  ,  THE  1 R  OERIVATIVFS  ANO  NEUMANN  FUNCTIONS. 

CAU  GENPSl 

C  PERFORM  HFSSEl  TEST  ANO  NEUMANN  TF  ST 

QUANRT  *  ARSI  CKR2*(PSSlSPI2»*CNFUMNm-R5SLSPIl  )*CNFUMNI2)  1-1.0) 
OUANNT  *  APS!  CKR2*( RSSISP INRANK |)*CNEUMN( NR ANK ) -R SSL S P| NR ANK ) *CNFll 

imninranki  n-i.:i 
IF! 0UAN«T-1.0E- 1 *<>360.352,352 
352  THTPRT  *  RTO*ThFTA 

PRINT  354, THT PP T , CKR. QUANRT, OUANNT 

356  FORMAT!  1H010X.13H*****  ThFTA  *F9.4,6M,  KB  *F1C<.4,1SH,  RF5SEL  TFST 
1*E12.5,16H,  NFUMANN  TEST  *E12.S,6M  **•*•) 

GC1  TO  36  2 

360  IF  I  QUANNT  - 1.0  F- 10M62.  3  52. 352 
362  CROW  «  0.0 
CPOWM  *  CMV 
IMP  *  2 

00  600  I  ROW  =  1 . NR  ANK 
CROW  x  CROWE  1 .0 
CROWM  <=  CROWMC1.C 

C  SET  UP  A  l OOP  FOR  FACH  COLUMN  OF  THF  MATRfCFS. 

CCOl  *  0.0 
CColm  =  cmv 

GO  TO  I  364,346), |M6 
364  J«R  =  1 

|m„  i 

Cn  Ml  368 
366  J  MR  *  2 
IMR  *  1 

36R  00  40C  TC.Ol  •  l ,  NR  ANK 
CCOl  x  CCOl E 1 .0 
CPU  *  CPOWECCOl 
CRSSIJ  x  CPOw*CCnt 
CCOIM  x  CCOl  ME  1.0 

ccnu  x  ccoi  ei. c 

I F  (  |H(")V-7)  372,  369,  >T? 

369  IF  (  ICCt -IRIIW  )370.  ?72,  372 

370  GO  TO  (  390,  3<ir|  , 

372  RJltP  »  CCOl  l*ntKR*8SSl SP ( ICOLEl )*SINTH 
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c 

c 


c 


c 


c 


c 


c 


c 


8JBJ1  *  CKR*ICKR*BSSLSPI ICOL I -C COL *BSSl SP I ICOIC1I ) 

If I IBOOY-91376, 376, 376 
376  GO  TO  I376.392I.JMR 
TFST  fOR  M  «  0. 

376  IFICMVI388.388, 378 

CALCULATE  THE  TERM  FOR  THE  CURRENT  ELEMENT  IN  THE  I  MATRIX. 

378  TERMI  *  SINTH*CKR2*8SSlSP<tCOLCll*ICOSTH*PNMUG(IROwCl»*PNMiLGI  ICU 
1LC1  l*CRIJ-CROMM*PNMlLG<  ICnLtll*PNMLLGI  |ROM»-r.CniM*PNMLLGCIRnNf.n*P 
2NMI  LG  (  ICOL  1 1 

AMX  I  HI  ROW,  ICOL  »  «  AM  XIII  IRON,  ICOL  I  C'  RHUl  *CNF  IIMNI  I  ROWtl  I  *TERM| 

IFI ICOL" I  ROM  I  388, 386, 386 

386  AMXIRI  IRON,  ICOL  I  *  AMX  |R  (  IRDW,  ICOL  I  C  SR  Mill  *HSSL  SP (  IR  OMCl  I *TFRMI 
388  IFI  IBO0Y-9I390.  3<52,  390 
390  JMR  «  2 


GO  TO  60C> 

CALCULATE  TERM  FOK  CURRENT  ELFMFNT  IN  THF  J  MATRIX. 

392  PTJ1  «  PNMllG<IRPMCU*PNHlLGIICmr.l»*I8JRJl*(CM2tCRSSIJ*COSTH**2» 
1CCRSSIJ*C0STH*8J1XPI 

PTJ2  «  CRONM*CCOL*PNMll  Gl  irom i *pnmllg  < ICOL  Cl I*IC0STH*RJBJ1CRJ1XPI 
PTJ3  ■  ccolm*pnmhgi  icoi I *H JBJ1 *  1C ROMM*PNMllG < IROMI -CPOM*COSTH*PNM 
illgi  iromuii 

AMX J I C IRON, ICOL  I  «  AMX J  1 1  IRON, ICOL  I  E  SRMUl *  S I N  TH*C  NF  UMN I I  ROW Cl  I  *  I PT 
1J1-PTJ2CPTJ3I 
IFI I  COL- 1  ROM  I  39  8, 396,  396 

396  AMXJRI  IROM,  ICOL  »  *  AMX  JR  I  18  ON  ,  IC  OL  » t  SR  Mill  *9 1  NTM*BSSL  SP  (  I  ROWE  II  *  I  PT 
1 Jl-PT  J2EPTJ3 I 


398 

600 

600 

700 

800 


808 

812 

816 


820 


JMR  «  1 

CONTINUE 

CONTINUE 

CONTINUE 

CONTINUE 

SYMMFTRIZE  REAL  MATRICES  AND  IMAGINARY  SPHERICAL  MATRICES. 

00  816  IROM  ■  2.NHANK 

ie; DSY  -  IROM-1 

00  812  ICOL  «  1,1  ENOS Y 

AMXIRI IROM, ICOL I  ■  AMXIRI ICOL, IROMI 

AMXJRI IROM, ICOL I  ■  AMX JR  1 ICOL , I  ROM ) 

TEST  FOR  SPHERICAL  BODIES. 

IFI  IBOOY-7IR12,  808,  8J, 2 

AMXI II IROM, ICOL )  ■  AMX1  II  ICOL, IROMI 

AMXJII  IROM, ICOL  )  ■  AMXJII  ICOL, IROMI 

CONTINUE 

CONTINUE 

SUMMATION  FOR  ALL  MATRIX  ELEMENTS  COMPLETE.  FINISH  PROCESSING  THEM 
00  860  JROM  ■  l.NRANX 
DO  820  JCOl  ■  l.NRANK 


AMX IRI JROM, JCOL I 
AMX III JROM, JCOL I 
AMXJRI JROM, JCOL I 
AMXJII JROM, JCOl I 
COMPUTE  K  MATRIX 
AMXKRI JROM, JCOL I 
AMXKI I JROM, JCOL I 
CALCULATE  THE  l 


«  CMV* AMXIRI JROM, JCOL  I 
*  CMV* AM  XI I  I JROM, JCOL  I 

■  OEM*AMXJR I JROM, JCOL  I 
-  OEM*  AM XJ 1 1  JO OM, JCOl I 

AS  A  FUNCTION  OF  THE  J  MATRIX. 

■  -AMX JR  I JROM, JCOL  I 
«  -AMXJI I JROM, JCOL I 

MATRIX  AS  A  FUNCTION  OF  THF  I  MATRIX. 


AMXIRI JROM, JCOLI  * 
AMXL 1 1 JROM, JCOl I  * 
CONTINUE 

COMPUTE  AOOITIONAl 
CKROM  ■  JROM 
IFIK MV  1826, 82 6, 826 


-AMX IRtJROW, JCOLI 
-AMX 1 1  I JROM, JCOLI 

TERM  FOR  THE  IMAGINARY  PART  OF  THE  K  MATRIX. 
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B24  FCTKl  «  l.C 
GO  TO  840 

826  1F( JROW-KMVI 828, 830,  830 
828  CMXNRMI JPOWl  *  1.0 

Gn  TD  860 

830  IBFCT  =  JRMW-KMVS1 
IFFCT  *  J ROWtKMV 
FPROO  *  IBFCT 
FCTKl  *  1.0 

CO  832  LFCT  *  IBFCT,  IEFCT 
FCTKl  *  FCTK 1 *FPROO 
FPROO  *  FPPOOGl.O 
832  CONTINUE 

840  CMXNRMI JROW I  *  4. 0*CKR0W* (CKPOWt 1 . 0 1 *FC TK I / I F M* 1C KP OWCCKROWR 1 . 0 )1 
AMXKl(JROWfJROW)  =  AMXK  |  (  JROW,  JPOW)  £Rr>VFCT*CMXNRM(  JRtlW) 

CMXNRMI JROWI  *  SORT ( C  NX  NR  M( JROW 1  I 
860  CONTINUF 

C  PROCESS  COMPUTED  MATRICES 

CALI  PRCSSM 
<*00  CONTINUE 
GO  TO  20 
ENO 
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C  A  PROGRAM  TO  flEAO  INPUT  OATA  FOR  THF  SCATTFRING  PROGRAM. 

SUBROUTINE  R  0  OAT  A 
COMMON  DTR ,R T O, CP  I 

COMMON  /CMVCOM/  NM, CM ( ( 30 » ,CMV» KMV , CM2 , EM , OF M , TWM, PROOM 
COMMON  /FNCCOM/  PNMLL G( 6 1  I , B S SL SPt  61 I ,C NF OMNI 61* 

COMMON  /MTXCON/  NP ANK ,NRANK I , CMTXRI  I  120 , 1 20  * , CMT X  I M 1 1 20 , 1 20  * , SPR MT 

1  X (  120, 120  I .CMXNRM (601 

COMMON  / THTCOM/  THETA, NTHET A, OL THT A , S I  NTH ,COS TH , I SMRI  , I SWTCHI 7 ) , SR 
1MUL, SMULSSI7I ,COHl6l,OHM, NSEC T, NDP S I  6 ) , EPPS  1 6  * , KSEC T 
COMMON  /RDYCOM/  CKR ,OCK R , CKR 2 ,C SKR  X , SNK R X ,C ONK , PR  * T , AL PH A, ( BODY, OR 
1,SNALPH,CSALPH 

COMMON  /TOTCOM/  ACANSI6C,  10  I ,' ST  $FC  T  , R  T SF r  1 
COMMON  /UVCCOM/  UANGI60),NUANG 
OIMENS (ON  CLRTOTI600! 

EQUIVALENCE  (  AC  ANSI  1,  II  ,C.  LRTOTI 
OIMENS ION  EPOEGIIOI 
C  R E AO  NECESSARY  INPUT  DATA 

PRINT  40 

40  FORMAT! 1H1119X, 3HCVM/ //////////////////////// 1H0 39X ,*0H* ********** 
1*****************************/1H02BX,67HELECTROMAGNFTIC  SCATTFRING 

2  FROM  GENERAL  AX  I SYMMFTRt  C  CONDUC T OR S/ 1 HO 3 9 X,40H*** *••****• ******* 

3  ******•••  . . I 

RFAO  80, NM, NR ANK, NS EC T, IROOY.NUANG 
80  FORMAT  I  6 1  12  I 

NR ANK I  =  NRANKC1 
PRINT  88 

88  FORMAT! 1H129X.75H  CASES  MATRIX  RANK  SECTIONS 

1  BODY  SHAPF  U  VECTOR! 

PRINT  92, NM, N RANK, NSEC T, I BOOY.NUANG 
92  FORMAT! 1H029X.5I1SI 

READ  96, CONK, BRXT, ALPHA 
96  F0RMATI6F12.il 
RTSFCT  *  8.0/CONK 
STSFCT  *  2.0*RTSFCT/C0NK 
PRINT  100 

100  FORMAT! 1H029X ,60H«ODY  PARAMETFRS  K|A)  BFTA/RHO 

1  ALPHA! 

PRINT  104, CONK, BRXT, ALPHA 
104  FORMAT! 1H044X.3F15. 3! 

READ  96, I  CM! |  (1,1  *  1,NM| 

REAO  80, I NDPS ( I » , I  »  l.NSECTI 
FRINT  120,  (NDPS!  I  !, !  *  l.NSECTI 
120  FORMAT ! 24HC  INTEGRATIONS/SECTIONRI 12 , /( 1H023X, 81121  1 
REAO  96,(UANG! I », I  *  l.NUANG! 

C  CLEAR  AREA  WHICH  CONTAINS  RUNNING  TOTALS. 

OO  136  I  *  1,600 
CLRTOT ! I !  *  0.0 
136  CONTINUE 

C  COMPUTE  END  POINTS  FOR  THFTA. 

ALPHA  *  OTR* ALPHA 
CALL  CAL  ENP 
00  14C  I  *  1 , NSEC  T 
EPOEG!  I  I  «  RTO*EPPS! I  I 
140  CONTINUE 

PRINT  148, (EPOEG! II,!  «  l.NSECTI 
148  FORMAT! 24H0  E NO  POINTS8F12.4,/! 1H023X .RF12.4I I 

RETURN 
ENO 
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C  A  ROUTINE  TO  COMPUTE  A  BESSEL  FUNCTION  OF  SFT  ORDER  AND  ARGUMENT . 

SUBROUTINE  BESSEL (NOROER, ARGMNT, ANSMR , I FRROR! 

OOUBLF  PRECISION  ARGMNT , ANSMR , X ,CN , SUM, APR, TOPP ,C I , CNI , ACR , PROD, 

1  FACT 

(ERROR  «  0 

N  *  NORDER 

X  *  ARGMNT 

CN  =  N 

SUM  «  1.0 

APR  «  1.0 

TOPR  «  -0.500*X*X 

Cl  «  1.0 

CNI  «  2*N£  3 

on  60  i  '  l.ioo 

ACR  «  TOPR*APR/(C  I* CNI  ) 

SUM  «  SUM £ ACR 

IFIOABSI ACR/SUMI-1.0D-20) 100, IOC, 40 
40  APR  *  ACR 

Cl  *  CICl.CDO 
CNI  «  CN I £2.000 
60  CONTINUE 
I FRROR  *  1 
GO  TO  200 

C  THE  SERIES  HAS  CONVERGED. 

100  PROD  «  2*N£1 
FACT  -  1.0 
IFINI160, 160, 120 
120  00  140  IFCT  -  1 , N 
FACT  ■  FACT*X /PROD 
PROD  *  PROD-2. ODO 
140  CONTINUE 
160  ANSMR  ■  F  ACT* SUM 
200  RETURN 
ENO 
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C  A  ROUTINE  TO  GENERATE  LEGENDRE  POLYNOMIALS. 

SUBROUTINE  GENLGP 
COMMON  OTR,RTD*  CP  I 

COMMON  /CMVCOM/  NM.CM II 30  ),CMV, KMV , C M2 , EM , QE M, TWM, PRDDM 
COMMON  /FNCCOM/  PNMLL Gl 61 I,B SSL SPI  61 1 .CNFUMNI 611 

COMMON  /MTXCOM/  NRANK , NR A NK I ,CMTXR L I  120, 1 20 » ,CMTXI M II 20 , 1 20 1 , SPRMT 
1X1 120,120) , CMXNRM ( 601 

COMMON  /THTCOM/  THE TA , N THETA, OL THT A , SI  NTH ,CQSTH, I SMRL , I SWTCHI 7 ) ,SP 
1MUL,'.MU».SSI7),C0H(6),0HM,NSECT,NDPSI6)  ,  EPPS  16)  ,KSECT 
COMMON  /BOYCOM/  CKR , OCK R, CKR2  ,C SK"  X, SNKRX.CONK.BRXT, ALPHA, I  BODY, OB 
1 , SNALPH, C  S ALP  H 
OTWM  -  TWMC1.0 
IF! THETA) 16,4,16 
4  IF  I  KM V- 1 ) 6 , 1 2 , 6 
6  00  6  ILG  -  1 , NRANK I 
PNMLLGI  ILG)  -  0.0 
8  CONTINUE 
GO  TO  88 

12  PNMLLGI 1 )  >0.0 
PNMLLGI 2)  -  1.0 
PLA  «  1.0 
GO  TO  48 

16  IF ( KMV 1 20 , 20 , 40 
C  THE  SPECIAL  C  AS  F  WHFN  M  *  0. 

20  PLA  »  1.0/SINTH 
PL R  ■  COSTH*PLA 
PNMLLGI 1)  «  PLA 
PNMLLGI 21  »  PLB 
IBEG  «  3 
GO  TO  60 

C  GENERAL  CASE  FOR  M  NOT  EQUAL  TO  0. 

40  OO  44  ILG  »  l.KMV 
PNMLLGI ILG)  «  0.0 
44  CONTINUE 

PLA  »  PROOM*S INTH**|KMV-1  ) 

PNMLLGI K M VC  1  )  «  PLA 
4B  PLB  «  DTNM*CO ST H*PL A 
PNMLLGI KMVC2 )  -  PLB 
IBEG  -  KHVC3 

C  00  RECURSION  FORMULA  FOR  ALL  REMAINING  LEGENDRE  POLYNOMIALS. 

60  CNMUL  *  IBEGCIRFG-3 
CNM  «  2.0 

CNMM  *  OTWM 

00  BO  ILGR  «  IBEG, NRANK I 
PLC  «  ICNMUL*COSTH*PLB-CNMM*PLA) /CNM 
PNMLLGI ILGR)  *  PLC 
PLA  *  PLB 
PLB  «  PLC 
CNMUL  *  CNMULC2.0 
CNM  «  CNM  C  1.0 
CNMM  «  CNMMC1.0 
80  CONTINUE 
BB  RETURN 
ENO 


72 


T  SUBTYPE. FORTRAN,  LMAP.L STRAP 

C  A  ROUTINE  TO  00  PINAL  PROCESSING  ON  THE  SCATTERING  MATRIX. 

SUBROUTINE  AOOPRC 
COMMON  OTR.RTO.CP I 

COMMON  /CMVCOM/  NM.CMII  30  I , CMV, KMV ,CM2 ,FM ,QFM, TWM , PROON 
COMMON  / FNCCOM/  PNMLLGI 61  ),BSSL SP( 61) .CNFUMNI61 ) 

COMMON  /MTXCOM/  NR ANK.NRANK I.OMTXI 1 160.60) .OMTXJI 160,631 .QMTXKI 160 
1.60 I.OMTXL 1(60.60) ,PM XII 6C.  60 1.PMX2 1 60, 601 ,PMX3 (63.601 .PMX4I63.60) 
2 » OMTX |R( 60,60 ) » OMTX JR  I  60, 60 ) .QMTXKR 1 60,60 ) .OMTXLR I 6  0 ,60 ) .CMXNRMI60 
3) 

COMMON  /VCMCOM/  ISYBG, JSYBG.KSYBG.NSYMT 

DIMENSION  EGVECTI 2, 120, 2 )  , TCMPt XI 2  , 1 20, 120) .FGMULI 120,2) ,FG ANSI  60, 
1 10) 

EOUIVALENCf  (OMTX  I  l  ,  F  GV  EC  T ) ,  I  PMXi  ,  TCMPL  X )  ,  (QMTXJl  ,FGMl)l >  .IQMTXK  I  ,  F 
1GANS) 

COMMON  /THTCOM/  THE TA , NTHf T A ,01 THTA , SI  NTH ,C OS TH, I SMRl , ( SWTC  117), SR 
1MUI .SMULSSI 7) , CONI  6), OHM, NSEC T.NOP  S I  6 ) , FP PS  1 6 ) ,KSFC  T 
COMMON  /TOTCOM/  AC ANS I  60 ,  10 ( . ST SFC T ,R T SFC T 
COMMON  /UVCCOM/  UANGI 60  I  ,  NUANG 

COMMON  / B CYCOM/  CKR,OCKR,CXR2.CSKRX,SNKRX,C0NK ,BRXT  ,Al PHA.IBODY.OB 
l.SNALPH.CSALPH 

C  NORMALIZE  AND  STORE  SECTIONS  T1  AND  T3  OF  THE  COMPLEX  T  MATRIX. 

00  60  IC  «  l.NRANK 
00  20  IR  *  l.NRANK 
JR  *  IRCNRANK 

OUANNM  «  CMXNRMI IR ) /CMXNRMI IC ) 

TCMPLXI 1, IR, IC)  *  OUA  NN M* QM TX I R I  I R , I C ) 

TCMPLXI1.JP.IC)  ■  OUANNM*OMTXKR(IR,IC) 

TCMPLXI 2, IR. IC)  *  OUANNM* QMTXIIIIR.IC) 

TCMPLXI2.JR.IC)  =  OUANNM*OMTXKI  !  IR.ID 
20  CONTINUE 
60  CONTINUE 

C  NORMALIZE  AND  STORE  SECTIONS  T2  AND  T6  OF  THE  COMPLEX  T  MATRIX. 

00  BO  IC  *  1 , N< ANK 
JC  «  ICCNRANK 
00  60  IR  *  l.NRANK 
JR  *  IRCNRANK 

OUANNM  *  CMXNRMI IR | /CMXNRMI IC) 

TCMPLXI 1, IR.JC)  *  QUA  NNM*  OMTX JR  1 1 R , I C ) 

TCMPLXI 1, JR, JC)  =  QUANNM*QMTXLR( IR  » IC I 
TCMPLXI 2, IR.JC)  =  OUANNM*  OMTX J 1 1 IR , IC) 

TCMPLXI  2,  JR,  JC)  «  QUANNM*OMTXL!UR,ir) 

60  CONTINUE 
BO  CONTINUE 

C  SET  UP  A  LOOP  FOR  ALL  VALUES  OF  THE  ANGLF  U. 

00  600  III  *  1, NUANG 

C  GENERATE  LF.ENORF  POLYNOMIALS  AND  OERIVITIVES,  RESET  THE  LIST. 

I F I  UANGI  IU))96,0B,96 
Bfl  COSTH  «  1.0 
92  SINTH  =  0.0 
THETA  *  0.0 
GO  TO  112 

96  IFIUANGI IUI-100. 0)106, 103, 106 
100  COSTH  «  -1.0 
GO  TO  92 

106  THETA  *  OTR*UANG(IUI 
SINTH  «  SINITHETA) 

COSTH  *  COSITHETA) 

112  CALL  GENLGP 

00  120  IPS  =  l.NRANK 
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EGMUL ( IPS • 1 )  «  CMV*PNMLLGIIPStll 
CPS  «  IPS 

FGMUL  I  IPS  »  2)  «  CPS*COSTH*PNMLLGI  IPSC1 l-ICPSCCMVI*PNMLLGIIPSI 
JPS  «  IPS CNR ANK 
FGMUHJPS.il  «  FGMULIIPS.2I 
FGMUL I JPS  ,21  «  FGMULIlPS.il 
I  20  CONTINUE 

C  MULTIPLY  THE  T  COMPLEX  MATRIX  TIMES  THE  LFGFNORF  VECTORS. 

KMVM1  «  IKMV-ll/4 
KMVM1  *  4 *KM VM 1 
IFIKMVM1I  132.  132.  124 
124  00  128  12  -  l.KMVMl 
FGVECTI  1,  12,11  «  0.0 
FGVECT I  2,  12,11  *  0.0 
FGVECTI  1,  12,2  1  =■  1.0 
f.VFfTIZ,  12,2  I  =  1.0 
JZ  «  I2CNRANK 
FGVECTI 1,J2, 1 1  *  0.0 
FGVECTI 2, J2, II  =  0.0 
FGVECT ll.J2.2l  =  0.0 
FGVECT 1 2, J2 • 2  I  -  0.0 
128  CONTINUE 

ISVBG  «  242*IKMVM1I 
JSVRG  *  KMVM1 
KSVBG  =  2*KMWM1 
NSYMT  «  NRANK-KMVM1 
GO  TO  136 
132  ISYBG  -  0 
JSYBG  «  0 
KSYBG  *  0 
NSYMT  *  NRANK 
136  CALL  VEC.MUL 

C  A  LOOP  TO  2FR0  CURHFNT  SUMS  OF  SCA T 1 , 2 ,  THTALl,?  AND  RTRA01.2. 

DO  14C  12  *  1,10 
FG ANS I IU , I  2  I  «  0.0 
140  CONTINUE 

C  SET  UP  LOOP  FOR  CURRENT  VALUES  OF  THE  SUMS. 

IPTH  *  1 

00  200  ICMS  *  1, NRANK 
JCMS  «  I C MStNRANK 

C  COMPUTE  S CAT T 1  ANO  SCATT2  SUMS 

EGANS  I IU, II  «  FGANSI I U, 1 1  Cl FGVECT 1 1 , I  CMS . 1 1 **2CFGVFCT 12, ICMS , 1 1  **2 
ICFGVECTI  1 , JCMS, ll**2CFGVECTI2,JCMS,l l**2l /CMKNRMI ICMS 1**2 
FGANSI IU , 6 1  *  FGANSI I U, 6 1 Cl FGVECT  I  1 , JCMS ,2 1 **?CFG VFCT 1 2, JCMS ,  2 1  **2 
ICFGVECTI 1, ICMS, 2I**2CF GVECTI2, ICMS, 2I**2I ZC MX NR MI ICMS 1**2 

C  FORM  THE  REAL  AND  IMAGINARY  PARTS  OF  T0TAL1.2  ANO  RTRAD  1,2 

PFR1  «  FGVFCTI 1, ICMS, 1 l*FGMUl I ICMS, 1 1 
PFI1  «  FGVECT 12, ICMS, ll*FGMULIICMS,ll 
PER 2  *  -FGVECTI  2, JCMS, 2  1*FGMUL IICMS.il 
PF  1 2  «  FGVECTI l.JCMS, 21 *FGMUL l ICMS , 1 1 
PGR  1  «  FGV EC  Til, JCMS, II *F  GMUL I JCMS , 1 1 
PGI 1  *  FGVFCTI2, JCMS, 1 l*F GMUL I JCMS , l ) 

PGR2  «  -FGVECTI  2, ICMS, 2  l*FGMUl I JCM  S ,1 1 
PGI 2  «  FGVFCTI 1, ICMS, 2>*EGMUL1JCMS,1 I 
GO  TO  I15C, 154, 15R, 1621, IPTH 
150  SGn  -  Cl.O 
IPTH  *  2 
GO  TO  170 
154  SGN  «  -1.0 
IPTH  «  3 
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GO  TO  180 
158  SCN  =  -1.0 
IFTH  =  4 
GO  TO  170 
162  SGN  *  C1.0 
1PTH  *  1 
GO  TO  1  HO 

C  CASE  FOR  N  MOO  4  IS  1  1-1, d)  flR  3  l£l. -II 

1*0  FGANSIIU.2)  =  f G4NSI IU, 2) £SGN*< PF I  1 CPGR1 ) /CMXNRMI  ICMS) **2 
F GAMS (  I U »  31  *  FG4NSI  111,  31-SGN*!  PER  l-PGtl  ) /CMXNRMI  !C.MSI**2 
FGANSI IU. 7)  =  FGANSI  111,  7)  CSGN*<  PF  I  2-PGR2) /CMXNRMUCMS)  **2 
FGANSI  III,  8)  =  FGANSI IU, 81-SGN*! PFR2CPGI2) /CMXNRM(ICMSI**2 
FGANSI  IU»  4  )  =  FGANSI  IU,  41-SCN*!  PF|  l-PGRl  ) /CMXNRMI  KMS)**2 
FGANSI IU. 5)  =  FGANSI  I  U»  5)  CSGN* I  PFR1CPGI1)  /CMXNRMI  1C  MS  1**2 
FGANSI  IU, <»1  *  FGANSI  III, 9)-SGN*IPF|2£PGR2)/rMXNPM(  |CMS)**2 
FGANSI  III,  1C)  =  FGANSI  IU,10)£SGN*(PFR2-PGI2I/CMXNRM( |CMSI**2 
GO  TO  200 

C  CASE  FOR  N  MOO  4  IS  2  (-1,-1)  OR  4  (£1,C1) 

180  FGANSI  IU, 2)  =  FGANS I  1 1), 2 ) CSGN* I PFR 1 -PG 1 1 ) /CMXNRM I  1C  MS ) * *2 
FGANSI  IU,  3)  -  Ff.ANSI  IU,  31  £SGN*(  PF  I  1CPGR1  I  /CMXNRMI  ICMS  )**? 

FGANSI  IU, 71  =  f  GANS 1 1  U, 7) CSGN* I PFR  2CPG 1 2 ) /CMXNRM ( | C  MS  I **2 
FGANSI  IU, 8)  =  FGANSI IU,flltSGN*(PFI2-PGR2)/CMXNRM(!CMSI**2 
FGANSI  111,4)  =  FGANSI  IU,4)£SGN*(PFR1  £PGI  l) /CMXNRMI  ICMS)**2 
FGANSI  IU, 51  =  FGANSI IU, 5) CSGN* I PF I 1-PGR 1) /CMXNRMI ICMSI**2 
FGANSI IU, 9)  =  FGANSI IU, 9) CSGN*( PFR2-PGI2) /CMXNRMI ICMS)**2 
FGANSI  IU.  1C  I  =  FGANSI  IU,1G!CSGN*|PF|2CPGR2)/CMXNRM(  |G.mSI**2 
200  CONTINUE 

C  SCALF  ACCUMULATIVE  SUMS 

FGANSI  IU,1)  =  STSFCT*FGAN SI  IU, 1 ) 

FGANSI IU, 21  =  STSFCT*FGANS( IU.2) 

FGANSI  IU,  3)  *  STSFCT*FGANS<  III,*) 

FGANSI  IU. 4)  *  PT  FCT*FGANS( IU,4) 

FGANSI  IU, 8)  =  R  iSFCT*FGANS( IU.8) 

FGANSI  IU, I.)  =  STSFCT*FGANS(  IU,6) 

FGANSI  IU, 71  *  ST  SEC T*FG AN  SI  IU,7) 

FGANSI IU. 8)  *  STS FCT*FG ANSI IU,8) 

FGANSI  IU, 9)  =  RTSFCT* FGANSI  IU,9) 

FGANSI  IU,  ICt  =  RTSFCT*FGANS( IU.10) 

400  CONTINUE 

C  PRINT  PARTIAL  SUMS  AND  ACCUMULATE  TOTALS, 

01)  50C  IPR  =  1,2 
PRINT  420 , KMV , I  PR 

420  FORMAT! 1H1 37X , 3 1 h** ****** **  CURRENT  SUMS  FOR  M  =|3,11H  ***♦**•**•/ 
1 /6H0CL ASS  12, 82H  ANGLE  SCATTERING  TOTAL  I  RE  AL )  TOTAL! I  MAG 

2)  RTRADIREAl  |  R  TR  AO  I  I  MAG  )//  ) 

IRFG  =  1 C  5* I  IP'-l  1 
IFNC  *  IR6GC4 
00  460  IUP  *  1 , NUANG 
on  432  ICAL  *  I  RFC.,  IF  NO 

ACANSI  up,  ICAl  I  -  ACANSI  IUP,  ICADf.FGANSI  IUP,  ICAL  ) 

432  CONTINUE 

PRINT  440, UANGI IUP), (FGANSI IUP, LP) ,LP  *  IREG.IEND) 

440  FORMAT!  1H  F 14 .2, 1P7E 1 5. 61 
460  CONT  INUE 
500  CONTINUE 

C  PPINT  THE  ACCUMULATE  TOTALS, 

DO  600  JPR  =  1,2 
PRINT  520 , X  MV , J  PR 

520  FOKMATI 1H135X,38H**********  ACCUMULATED  SUMS  FOR  M  «I3,11H  ******* 
1 ***//6H0CL ASS  12, 112H  ANGLE  SCATTERING  TOTALIREAl)  TOTAL 
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21 (MAGI  RTRADIREAL)  RTRAOI (MAG)  RCS  PHASE  ANGLE 

1//) 

JBEG  «  It  5*1 JPR- 1 1 

JENO  -  JREGC4 

00  560  JUP  •  l.NUANG 

COMPUTE  PHASE  ANGLE  AND  RCS  1  OP  ?. 

PHANG  ■  RT0*ATAN2! ACANS(JUP.JEND) .ACANSI JUP,JFN0-1>) 

RCS  1 2  «  ACANS ( JUP  * JENO- 1 >  **2C ACANS I JUP  »  JF  NO )**2 
PRINT  440  «UANG( JUP) *  I  AC AN  SI JUP  »l  P) »l P  *  JBEG* JENO I • RCS 12* PHANG 
540  FORMAT! 1H  F 1 4. 3 , 1P7F 1 5. 6 ) 

560  CONTINUE 
600  CONTINUE 
RETURN 
FNO 
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T 

C 


SUBTVPE, FORTRAN, LMAP, L STRAP 
A  SUBROUTINE  TO  PRINT  OUT  THE  T  MATRIX. 

SUBROUTINE  PR TR I T 

COMMON  /MTXCOM/  NRANK  ,NRANK  l  ,  QMTX  I  |  ( 60 , 60  I , OM T X J I  I  60 ,60 )  .OMTXKl  (60 
1.6PI,OMTXLI(60,60I.PMX1I6C, 631 ,PMX 2 (60.601  . PM  X? 160 .601  .PMX4I60.60I 
2.QMTX  1RI60.60I.QMTXJR (60, 63  I , QM T XKR ( 60. 60 » . QMTXLR ( 60 ,60  I ,CMXNRM|60 
3  I 

EQUIVALENCE  (QMTX| I.CMTXRL I, ( PMX1,CMTXIM» , ( QMTX I R , SPRMT X » 

01  MENS  ION  CM TX»L( 120, 1201 ,CMTXIM( 120,123) .SPRMTXI 120, 120) 

PRINT  2H 

2fl  EOPMATI 1H1///1H0S2X, 16HMATRI X  Till, REAL) 

CALL  PR 1NTM( QMTX Ifl, NR ANK 1 
PRINT  128 

128  E0RM6T ( 1H1/// 1H0S2X, 16HMA TR I  X  T(2).RFALI 
CA(  L  PR (NTM( QMTXJP, NRANK I 
PRINT  2?R 

22B  FORMAT(1H1///1HOS'’X,16HMATH|X  T(3),REAL) 

CALL  PRINTMIQMTXKH, NRANK) 

PRINT  3?H 


32fl  FORMAT (  1H1 III 1H0S2X, 16HMA  TR 1 X 
CALL  PM (NT«( QMTXI R, N« ANK  ) 
PRINT  428 

428  FORMAT!  1H1///1H04RX.21HMATRIX 
CAl  L  PR (NTMIQMTX I  I ,NR ANK I 
PRINT  52B 

?28  FORMAT!  1M  1  / / / 1H049X , 2 1HMA  TR I  X 
CALL  PRINTMIOMTXJI, NRANK) 
PRINT  628 

628  FORMAT!  1H 1 // / 1H04  9X . 2 1HMA TR I  X 
CALL  PMINTMIQMTXK (.NRANKI 
PRINT  728 

728  FORMAT!  1H1///1H0<.9x,21HMATR|  X 
CALL  PRINTMIOMTXL (.NRANK) 

HF  T  URN 
FNO 


T I  4 1  .REAL  I 

Tilt , IMAGINARY) 

T( 2) , IMAGINARY) 

T(  3)  ,  IMAC.INARV) 

T(4)  .IMAC.INARV) 
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T  SUBTYPE. STRAP 

•  A  ROUTINE  TO  GfT  THE  F  l.C.l  A  NO  F2.G2  VFCTORS  FROM  T  *  Plffll  III. 
PUNREl 

PUNF PC  » L AST »  COMMON 
VFCNUL  FNT  ER  *  SVXR  S 
PUNCOC 

NTXCOM  COMflLOCK,  F  INAL 
VCMCOM  CCM8LDCK  ,  F  INAL  V 
XW, 0,0,0 

SVXRS  SX«$12,XR12  •  SAVF  I NOF  X  REGISTERS  12, 

SX,$13fX»U  *  13, 

SX,$14,XP14  *  14, 

SX,$2, XP2 
SX , (3, XP  3 
Sx.tlO.xPIO 

I  WFIUI  ,  IWEWT  y  *C-MO,nE  TOOPESs  E  \(  1  1W  toj  r,j  r>o 

0* I U ) , NP  AN  K 

SHFL , 3  R 

F-IIUT.3H 

STIU),  0,1F?A 

LX,$12,TMXWR  *RFSH  T  COMPLEX  MATRIX  CONTROL. 

VC , i 12 , ISYRGC  .32 

LCI.U2.2.0 

SX, $12 .Tmxw 

LX, $13, PMX  WR  *  RT  SF  T  MOLTIPIIFP  VECTOR  CONTROL. 

VC, $13, JSYHGCG. 32 
LC, $13.NSYMTC0.32 
SX, $13,PMXW 
LIUI.NSYmt 
ciim.nsymt 
LX,  $14,  FC.XWO 
VC,  $14,KSYRr,C0.32 
LC, $14, $1 C  0. 32 
SX,$14,F  GX  W 

FZIHU,  1,81  ,RFGGFC0.2v  ’  SF  T  SIGN  FOR  ANSWER  STORAGE,  FI  OR  G1 
FZ I RU, 1,8  1 , IFGGFC0.20 

•  00  THE  PEAL  PART  OF  FI  OR  Gl,  G2  OR  F2. 

HGANLP  LX, $2, $12 

VC,$2,G1F2ACJ.  32 

LX, $3, $13 
VC, $3, NPANKC 3, 32 
SRLF1R  KC, $14,NSYMTC0.32 
RZXF.NCCMPR 
VC, $12 .KSYRGC3.32 
SV, $12, TMXW 
VC, $  1 4 , KSYBGCO. 3  2 
LV,$2, $12 
VC, $2, GIF  2 AC 0. 32 
NOCMPR  L V I , $ 10,N IF1P 
OL  I  U I »  Z  FPO 
SMLF1R  8,0.01  $  1C  I 

N1F1R  LMPIIM  ,0.0I$12>  *N  MOO  4  «  I  OR  1.  3 

*CIN», C.OI $131 
LMPIUI ,  1.01 $2) 

•  Cl N I »  0.01 $3) 

LVI ,$10,N2F1R 
B.FMLF1P 

N2F1R  LMRI U ) , 1. 3 1  $  12  I  »N  MOO  4  »  -l  OP  I 

•NCI N ) ,0. 01  $  1 3 ) 

LNRUM  ,G.0I$2I 
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*£(NI »  0«0( *3) 

LVI.410.N3F1R 

B.EMLFIR 

N3F1R  LMRIUI ,0.0(412)  *N  MOO  4  «  -I  OR  -1.0 

♦NCINI, 0.0(4 131 
LMR(U),  1.0(42) 

♦N£ ( N I .0.0(131 
LV  I  • 4 10  »N4F1R 
B.EMLFIR 

N4F1R  LMRIUI ,1.3(412)  *N  Xdft  4  •  l  OR  -I 

*C ( N ) »  0 .0 (  SI  3  I 
LMRIUI ,0.0(12) 

*N£( N  )  ,  C.  0  (  4  3  ) 

LVI.410.N1F1R 

EMLF1R  V£I,412,MTXS7F 
VCI.42.MTXS7F 

vu.m.i.o 

VCI, 43, 1.3 
CBR • 41 3»  SMLF  IP 

RFGGF  SROINI ,0.3(4  14  I  •  S  T  OR  F  REAL  PART  OF  FI  OR  01,  02  OR  F2 

VCI  ,414,1.0 

•  no  THF  IMAGINARY  PART  OF  FI  OR  01,  02  OR  F2. 

IV, S3, 413 

VC ,  43,  NR ANKC  0. 3? 

LV,412,Tmxw 
IV, 42, 412 
V£,42,G1F2A£P.32 
LVI.410.N1F1  I 
OLCUI, ZERO 

SMLF1T  R. 0.0(4131  *N  MOO  4  »  |  OR  1.3 

N1F1I  LMRIUI, 1.0(412) 

*£(NI,0.0( 413) 

LMR(U) ,0.3(421 
•NC(N) ,0.0(431 
LVI.413.N2F1 I 
P.FMLFl I 

N2F1I  LMR(IJ), 0.0(4121  »N  MOO  4  *  -1.0  OP  I 

*£(N), 0.0(4131 
LMRIUI ,1.0(42) 

*£(N),0.0( 431 
LVI.410.N3F1  | 

R.FMLFll 

N3F1I  LMRdJI  ,  1.0(4  12  I  *NMOO  4  -  -I  OR  -1.0 

♦NCINI ,0. 3(4 131 
LMR(U1, 0.0(421 
*C ( N I , 0 .0 ( 431 
LV1.413.N4F1  I 
R.FMLFll 

N4F1I  LMRIUI, 0. 0(4121  ’N  MOO  4  ■  1.0  OR  -I 

•NCINI, 0.0(413) 

LMRIUI ,  1.0(4  2  I 
•NCINI ,0.0(43  I 
LVI.410.N1F1  I 

FMLF1I  VC1.412.MTXS7F 
VC  1 , 42 , MTX  SZ  E 
VCI, 413, 1.0 
VCI, 43, 1.0 
CRR ,41 3,  SMLF 1  I 

IFGGF  SRD(N) ,0.3(414)  • STORE  IMAGINARY  PART  OF  F1,G1,G2,F2. 

LV, 43, 413  ’SET  UP  FOR  NEXT  I TFM. 
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XR  12 

XR13 

XR  14 

XR  2 

XR  3 

XR  10 

TMXM 

PMXM 

FGXW 

TMXMR 

PMXMR 

FGXMR 

G1F2A 

TMFRTY 

ZERO 

HTXSZE 

VCTSI E 

l  AST 

COMMON 

NRANK 

NRANKI 

FGVECT 

OMIIRM 

FGMUL 

OMJIRM 

FGANS 

QMKIRM 

QMLITL 

TCMPIX 

CMXNRM 

FINAL 

isvbg 

JSYRG 


VG»$3,NRANKG0.32 

LV, $12,  TMXM 

VCI ,$12,2.0 

SV, $12. TMXM 

LV. $2. $12 

VC, $2, GIF 2 AGO. 32 

CBRG,$14,SRLF1R 

VGI,$14,MTXSZF 

LV, $12, TMXMR 

VC, $12,  ISVRGG0.32 

SV, $12, TMXM 

VG I • $  1 3, VC  TSZE 

SV,$13,PMXM 

F1IBU,  1,81 ,fl  FGGF  GO, 20 

FI! BU, 1,81 , IFGGFGO . 2C 

CR, $12,BGANLP 

LX, $2, XR2 


•PREPAWF  FOR  NFXT  ROM  Ilf  T  MATRIX. 


•RESET  XR 14  FOR  G2,F?  VECTOR 
*RF  SET  TO  FIRST  ROM  IN  T  MATRIX. 

•RFSFT  MULTIPLIER  VECTOR 

•SET  SIGN  FOP  ANSmFR  STORAGE,  F2  OR  G? 


LX, $3, XR3 
LX,$10,XR10 
LX, $12, XR 12 
LX, $13, XR 13 
LX, $14,XR14 
8,0.0! $151 

STORAGE  REQUIREMENTS  FOR  THE  VFCTOR  MULTIPLICATION  ROUTINF 
XM, 0,0,0  • SAVE  SPACE  FOR  INOFX  RFGISTER. 


XM, 0,0,0 
XM, 0,0,0 
XM, 0,0,0 
XM, 0,0,0 
XM, 0,0,0 

XM, TCMPIX, 0, TMXM 
XM, FGMUL, O.PMXM 
XM, FGVECT, O.FGXM 

xm,tcmplx,o,tmxm 

XM, FGMUL, O.PMXM 
XM, FGVECT, O.FGXM 
ORZIUl.l 

IF10I00IU). 240. 0X38 

OOINI.O 

SYN, 240.0 

SYN, 120.0 

SYN  ,  $ 

SLCRCOM 
SYN,  $ 

SLCRCOM, MTXCOM 
ORIUI,  1 
CRIU),  1 
OR  I N ) , 480 
ORINI, 3120 
OR  IN  I, 240 
ORINI, 3360 
ORINI, 600 
ORINI, 3COO 
CR  IN  I,  3600 
ORINI, 28800 
ORINI, 60 
SYN,  $ 

SLCRCOM, VCMCOM 
ORIUI,  1 
OR  I U 1 ,  1 


•INOFX  CONTROL  FOR  T  COMPLEX  MATRIX. 

• INOf  *  CONTROL  FOR  POLYNOMIAL  VECTOR. 
♦INOFX  CONTROL  FOR  VECTOR  ANSMERS. 


•STORAGE  FOR  INQEX  I NCREMFNTER  G1.F2. 
• INTFGFR  *  240 
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KSYBG  OPIUJ.l 
NSYMT  OOl'Jl.l 
MNM.V  SYS.* 

r*jr 
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T  SUBTYPF, FORTRAN, IMAP, ISTRAP 

C  A  ROUT  INF  TO  PRINT  OUT  A  MATRIX  ASSAY 

SUBROUTINE  PRINTM(P.N) 

DIMENSION  PI60.60) 

NR  «  N 

00  100  I  *  1 i NR 
IR  «  1 

20  IE  «  me? 

IFI IF-NR) 2fl«  2fl»  24 
24  IF  >  NK 
28  IFI  IB-1  136,36,60 
36  PRINT  44,  l,IP||,Jt, J  =  in, IE) 

44  FORMAT  I 5H0  R0W|3, 2X , lPflF 1 S. 6) 

GO  TO  HO 

60  PRINT  6fl,(P( I . J ) , J  »  |H, IF) 

68  FORMAT  I 1M  9X,1PHF1S.6I 

no  in  «  i ec l 

IFI  IB-NR )  20, 20,  100 
100  CONTINUF 
RETURN 
ENO 
ENO 
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T  SUBTYPE, FORTRAN, LMAP, l STRAP 

C  A  ROUTINE  FOR  BESSEL  FUNCTIONS,  DERIVATIVES  AND  NEUMANN  FUNCTIONS. 

SUBROUTINE  GFNBSL 
COMMON  OTR.RTO.CPI 

COMMON  /CMVCOM /  NM , CM  1 1  30  I , CMV, K MV , CM2  ,  EM.OFM, T WM , PROOM 
COMMON  /FNCCOM/  PNMLL G( 61 ) ,6 SSI SP( 61 1 ,C NFUMNI 61  » 

COMMON  /MTXCOM/  NRANK  ,  NRANK I , CMTXRL ( 120 , 1 201 ,CMTX I M 1 1 20 , 120  I , SPRMT 
1X1 120, 1201, CMXNRMI60I 

COMMON  /THTCOM/  THFTA,NTHFTA,OLTHTA,SINTH,COSTH,l SMRL, I SWTCHI 71 ,SR 
1MUL , SMULSSl 71, COM  161,  OHM, NSEC T,N0PS(6I , EP°S 1 6 1 , KSFCT 
COMMON  7BOYCOM/  CKR,OCKR,CKR2.CSKRX,SNKRX,CONK,BRXT, ALPHA, IBOOY.QB 
1 ,SNALPH,CSALPH 

DOUBLE  PRECISION  PCKR , ANS WR , AN SA , AN  SB , ANSC ,C ONN 
C  SFT  UP  A  LOOP  TO  GET  2  SUCCESSIVE  BESSEL  FUNCTIONS. 

NVAL  *  NR ANK- 1 
PCKR  *  CKR 
DO  40  I  »  1,4 

CALL  BESS  El < N VAL , PCKR , ANS WR , I fR R OR  I 
IF (  (ERROR  120, 20,  32 
20  ANSA  *  ANSWR 
NVAL  =  NVAltl 

CALL  BESS  El  ( N VAL , PCKR , ANS WR , I FR ROR  t 
IF ( IERROP I24,24,?H 
24  ANSB  =  ANSWR 
GO  TO  60 

28  NVAL  *  NVAl-1 
32  NVAL  >  NV  AL  CNR ANK 
40  CONT INUE 

C  PROGRAM  UNABLE  TO  GENERATE  BESSEL  FUNCTION. 

CALL  DUMP 

C  SET  UP  FOR  PROPER  RECURSION  OF  THE  BFSSFl  FUNCTIONS. 

60  IF ( NV  AL -NR  ANK I10C,1CO,64 
64  IENO  *  N V  Al —NRANK 
CONN  *  2*INVAL-1IC1 
DO  72  IP  *  1,  IFNO 
ANSC  *  CIINN*ANSA/PCKR-ANSB 
CONN  *  CONN-2. ODO 
ANSB  *  ANSA 
ANSA  *  ANSC 
72  CONTINUE 

C  PROGRAM  IS  READY  TO  RFCURSE  DOWNWARD  INTO  BFSSFl  FUNCTION  VECTOR. 

100  BSSLSPINR ANK I  I  *  ANSB 

BSSLSPINRANKl-ll  «  ANSA 
CONN  r  NR ANK6NRANK-1 
IE  =  NRANK  1-2 
JF  *  IF 

00  120  JB  *  1 , JE 
ANSC  *  CONN*ANSA/PCKR-ANSB 
BSSLSPI  IF  |  r  ANSC 
ANSB  *  ANSA 
ANSA  *  ANSC 
IF  *  IF-1 
CONN  =  CONN-2.000 
120  CONTINUE 

C  GFNFRAT E  THF  NEUMANN  FUNCTIONS. 

CMULN  *  3.0 

SNS  A  =  -CSKPX 

SNSB  =  -CSKRX/CKR-SNKRX 

CNFUMNI1I  »  SNSA 

CNFUMNI 21  =  SNSB 
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00  2A0  t  *  ?.NRANK| 

SNSC  ■  C«UlN*SNSH/f KB-SNSA 
CNFUMNt  1 >  *  SNSC 
SNSA  ■  SNSK 
SNSA  ■  SNSf 
C*nilN  -  CNULNC2.0 
2 R0  CI1NTINUF 
BFTIIRN 
END 


OUMP 

st  ART 

last 

COMMON 


SUP! VP  £ ♦ STRAP 

A  ROUT INF  TO  DUMP  CORF. 

PUNREt 

PUNFPC  ,  l  AST,  COMMON 

ENTER, START 

XW, 5,0,0 

R,*MCP 

, tAREOJ 

8,0. C(  US) 

SYN  *  % 

SLCRCOM, 

SYN,  $ 

FNO 


•RFTIIRN 
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T  SUBTVPE.F0R7RAN.LMAP.L STRAP 

C  A  ROUTINE  TO  ORTHOGAN I L 1 1 E  THE  0  MATRICES  TO  PRODUCE  T  MATRICES. 
SUBROUTINE  PRCSSM 

COMMON  SMTXCCJM/  NR  ,  NR  1 , 01 II  60, 6C I  .  01 2  (  63,601  .013(60  .601 , 01 4  I  AD ,  60  I 
l,Pl(60,60>,P2(60,60»,P3(60,60»,P4160,60) ,0R1< 60,601  , OR  2 (60, 60  I , OR  3 
2 (60. 60  I *0R4I 60. 60  I <  CMXNRM (601 
EQUIVALENCE  I  01 1, R I ! I . < OR l» RR1 » . ( TMMX , PI ) 

01  MENS  ION  RR1( 120. 120 I.R I  11120.1201  .TMMXI  120.1201 
C  NORMALIZE  AND  TRANSPOSE  THf  I.J.X.L  MATRICES  Tl)  OBTAIN  0  MATRICES. 

CALL  NRMOMX 

C  SET  UP  REAL  ANO  IMAGINARV  MATRICES  FOR  GFNFRAl  m  CASE. 

DO  6  I«1,NR 
H4-ICI 
00  4  J= 1 »  NR 
NN«JtJ 

TMMX(MM-l,NN-l»  =  Olill.JI 

TMMXIMM-l.NNI  =■  0121  I .J  I 

TMMXIMM.NN-ll  =  01 3(1 . J I 

T  J'MXI  MM  ,  NN  I  =  QI4I 1 ,J  > 

4  CONTINUE 
6  CONTINUE 
NBGR  «  NR  CNR 
00  10  I  -  l.NHGP 
00  8  J  =  1  .NBGR 
R I  1( I . J I  *  TMMXI  1 . J I 
8  CONTINUE 
10  CONTINUF 

DO  14  I  =  l, NR 
MM  *  It  I 
DO  l 2  J  -  l, NR 
NN  ■  JCJ 

TMMXI  MM-  1 ,  NN  -  1 1  =  ORKI.JI 
T  MMX  I  MM- 1 , NN I  =  0R2I  I ,J  » 

TMMX ( MM, NN-l I  =  0R3 1  I .J  I 
TMMXIMM.NNI  =  OR 4 1  I .J  I 

12  CONTINUE 
14  CONTINUE 

OO  18  I  =  l, NBGR 
00  16  J  -  l.NBGR 
RRK  1  •  J  >  *  TMMXI  I  ,J  I 
16  CONTINUF 
18  CONTINUE 

C  CONDITION  0  MATRICES  BEFORE  OR THOC ONAL I Z I NG  THEM. 

CALL  CNOTNQ 

C  NORMAl  I Z C  THE  NTH  ROW  OF  AN  N  N  MATRIX 

SUM l  =  0.0 
00  20  K  «  l.NBGR 

SUM l  *  RRKNBGR.K  !**2CR  1 1  (NBGR.K  »**2CSUMl 
20  CONTINUF 

SU«l  *  SORT! SUM  1 1 
00  28  K  *  l.NBGR 
RRKNBGR.K  )  *  RR1INBGR.KI  /SUM l 
RIIINFGR.KI  »  R  I  HNBGR.KI  /SUM  l 
28  CONTINUE 

C  SET  UP  A  LOOP  FOR  THE  N-l  REMAINING  ROMS. 

NMI  *  NBGR-1 
NROW  *  NBGR 
00  IOC  I  =  l.NMI 
NROW  •  NROW- 1 
MROW  *  NROW 
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on  36  K=  1 ,  <f,GR 

TMMX (  l.KI  *  RR1INROW.K) 

T MMX I  2.  K  I  *  Rill NROW*  K I 
36  CONT INUE 

DO  80  J  =  NROW.NMI 
SR  1  =  0.0 
S 1 1  *  0.0 
MROW  =  MR  OWC 1 
00  60  K  =  l.NAGR 

SRI  =  SP1CRR1  IMROW.K I*RR1 INROW.K  I CR 1 1 1 MROW, K I *R 1 1  (NROW,  Kl 
Sll  =  SI  1CRRJ  (MROW.KI  *H  II  INROW.K  I  —  R  I  1IMR0W.KI  *PR1  I  NROW,  K  I 
60  CONTINUE 

OH  68  K  =  l.NAGR 

TMMXll.KI  *  T  MM  X  I  1 »  K | — SR  1 *R R 1  I MROW , K I C  S I  1 *R I  1 1 MROW. K | 

TMMXI2.KI  *  TMMXI  2.K  l-SR  1  *R|  1  (MBOW.K  »- SI  1*«R1  IMROW.KI 
6fl  CONTINUF 
RO  CONTINUE 
Sum i  =  o.O 
DO  86  K  =  l.NAGR 

SUM1  *  SUM1CTMMXI 1,KI**2£ TMMXI?, Kl **2 
R6  CONTINUE 

SUM1  =  SORT  I SUM1 I 
00  88  K  *  l.NAGR 
RR1INP0W.KI  =  TMMXI l.K I /SUM1 
Rill  NROW.  K  |  =  TMMXI  2.KI  /StJMl 
88  CONTINUE 
100  CONTINUE 

C  PRINT  OUT  OR  T  HOGAN  1 L  l  /  FO  0  MATRICES 

PRINT  120 

120  FORMAT  I  1H160X.60HRFAI  SFCTION  OF  OR  THOGANl  l  I  ZEO  .3  MATRIX.  I 
CALL  PRNOOTIRRi.NAGRI 
PRINT  128 

12H  FORMAT  I 1H 1 ?7X  »65H|MAGINARV  SECTION  OF  OR TMOGA Nl l I 7F 0  Q  MATRIX.  I 
CALL  PRNOOTIRIl.NRGRI 

C  PERFORM  0  TRANSPOSE  •  RFALIOI  TO  GF  T  T  MATRIX. 

00  160  I  =  l.NAGR 

00  152  J  =  l.NAGR 

TMMXI  I, J  I  =  0.0 
162  CONTINUE 
160  CONTINUF 

00  180  I  =  l.NAGR 

00  176  J  =  l.NAGR 

00  17?  K  =  l.NAGR 

TMMXll.JI  *  TMMXI  I , J l-R | 1  IK , I  »*RR1  IK,  J» 

172  CONTINUE 
176  CONTINUE 
180  CONTINUF 

CO  196  I  =  l.NR 
MM  =  1C1 

OO  192  J  *  l.NR 
NN  =  JCJ 

01 II I.JI  *  THMXIMM-l.NN-1  I 
0I2II.J)  =  TMMXIMM-l.NNI 
0131  I.JI  =  TMMXI MW.NN-l I 
0161 I.JI  =  TMMXIMM.NNI 
192  CONTINUF 
196  CONTINUF 

00  208  I  *  l.NAGR 
00  206  J  =  l.NAGR 
TMMXI  I.JI  =  O.C 
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204 

CONTINUE 

20fl 

CONT INUC 

OD  220  I 

-  1,N«GW 

cn  216  J 

x  l.NRGR 

no  212  K 

*  1 1  NR  GR 

TMMXI 1, J) 

x  TMMXI  I,  JICPRI  I«  .1  »*R«»1  IK  ,  Jl 

212 

CONTINUE 

21ft 

CONTINUE 

2  20 

CONTINUE 

CO  236  I 

MM  x  !£! 

x  1 1  NP 

CO  232  J 
NN  *  JCJ 

x  1,N« 

CR1I It J) 

=  TMMX ( MM-l.NN- n 

Q«2I I.JI 

=  TMMX(MM-1,NN) 

QR3I I, j ) 

i  TMMX  (  MM,  NN-1  1 

QR4I I , J) 

x  TMMXI mm, NN» 

232 

CONTINUE 

2  34 

CONT INUE 

PRINT  THE 

T  MATRIX 

CAI  l  PM  T  R  IT 

C  CO  E  1NAI  PROCESSING 


CAIL  AOOPRC 

RFTlJWN 

ENO 
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T  SUBTYPE, FORTRAN, LMAP, L STRAP 

C  *  ROUTINE  TD  CONDITION  0  MATRICES  BEFORE  ORTHOGONAL I  2  I  NO  THEM. 

SUBROUT  INF  CNOTNO 

COMMON  /MTXCOM  /  NR,  NR  I,OI  II  60  ,  ftD  »  ,  0  I  2  (  *  T,  6  '  »  .(IPH  )  ,AD»  ,01  4(A~  ,A'  » 
1  *  P 1 1 6C  *  60  I  *  P  2 1 60, AC  I »  P3 ( A0»60 1  *  P4( 6C ,6 1 1  *  OR  1 ( Ao ,6 ) I  *  QR2 ( 6C  * 62  I » OR  3 
2I6C.60)  ,09M60,60I,C“XN9M(M) 

EQUIVALENCE  (01 1,»I1I,<0» l,RPl»,(TMMX,Pn 
DIMENSION  PP  1  ( 1 2o, 1 23  I «  I  1 (1 ?D, 1 20  I , TMMXC I  70 , 1 2C I 
C  SET  UP  LOOPS  POP  All  PUT  THE  FIRST  ROW, 

NBCP  -  NR  CNR 
NPOW  =  NBCP 
DO  60  KP  *  2  *  NP  GR 
C  RESCALE  THE  CURRENT  ROW. 

SCLE1  *  1  .O/PIKNPOW.NPOWI 
DO  8  1C  =  l.NRCR 

PR  1  (  NROW ,  L  C  I  =  SCIE1*RR1(NRUW,ICI 
PllINPOW.lCI  *  SC  l  F  1*R  Il(  NROW  ,(  f  I 
8  CONTINUE 

C  RESCALE  ALL  TH(  PIlwS  UP  TO  THE  f.ORRFNT  ROW. 

MRUW  =  NPOW-1 
00  20  hr  =  1 ,  MB (1W 
R  S  C  L  1  =  PIICMR.NROWI 
00  16  Mr  =  l.NP'-R 

RR 1 ( MR ,Mf  I  r  HP K MR ,MC l-M  SCL  1*RR 1< NRUW.MC » 

RI1(MP,MCI  *  PIl(MRfHCI-PSCLl*Rll(NROW,MC» 

16  CONTINUE 
20  CONTINUE 

NPOW  *  NROW-1 
60  CONTINUE 

C  SET  IMAGINARY  ELEMENTS  ABOVE  THE  MAIN  DIAGONAL  «  3. 

NROW  *  NBC.R-1 
CO  80  l  =  l.NROW 
18  *  IG1 

00  72  J  =  IH.NBGR 
RUII.JI  =0.0 
72  CONTINUE 
80  CONTINUE 
RETURN 
END 
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r  SUBTYPE, FORTRAN, ( MAP, l STRAP 

C  A  ROUTINE  TO  NORMALIZE  THE  l.J.K  ANO  L  MATRICES  TO  GET  A  0  MATRIX. 

SUBROUTINE  NR  MO MX 

COMMON  /M  IXCOM/  NR  ANY , NR  ANY  I  ,  AMx  |«  (  60  ,  6  31  ,  AMX  JR  (  6  3,  40  I  ,  AMXKR  (  60  ,  6  0 
11 , AMXt R( 60,60  I, AM  XI  I ( NO, 6 01 , AMXJI ( 60,401 , AMXKI <63,601 , AMXl I  I  NO, 63  I 
2.0MTX IRI 60,60  I .QMTX JR  I  60,  601  .OMTXKP (60,601  .OMTXLR (63,60  I ,CMXNRM(6J 
3  I 

FOU ( VALENCE  (AMXtP,OMTXIII,(AMXJR,QMTXJI!,<AMXKR,QMTXKI!,( AMXLR, OM 
1TXL  II 

HI  MENS  ION  0M1  X(  ((  60,631  .iJMTXJ  I  (  63,601  .OMTXK  I  (60,631  ,(JMTXL  I  (  63,60  I 

C  SET  UP  LOOPS  TO  PROCESS  ALL  ROWS  AND  COLUMNS  TOR  THF  REAL  MATRICES 

00  200  (P  =  1  ,  NR  A  NR 
00  100  (C  =  l.NPANY 
OUANNM  -  CMXNRMI  (P  )*CMXNRM(  (C.  ) 

CMTXIRI  IP  ,  If.)  =  A«X  ,l«  (  If  ,  IP  1  /OUANNM 
OMTXJRI  IR ,  IC I  *  -AMXLRI  1C ,|R I /OUANNM 
OM  T  XKR I  IR,  IC)  =  AMX  IR (  If,  I« 1 /OUANNM 
OMTXLRI  IR,  IC I  =  AMXYH ( IC,  IR 1 /OUANNM 

100  continue 

200  CONTINUE 

C  SET  UP  LOOPS  Of  ROMS  ANO  COLUMNS  E OR  THF  IMAGINARY  MATRICES, 

nn  40u  IR  =  1, NR  ANY 
on  300  (C  =  i.npany 
OUANNM  =  CMXNRMI IY )*CMXNRM{ IC > 

QMT  X  |  I  ( I R ,  IC 1  =  AMx. II  I  IC ,  IR I /OUANNM 
OMTXJI!  IR,  (Cl  *  -AMXL  H  IC  ,  IR  I /OUANNM 
OMTXK  I  (  (R ,  C >  -  AMX  || ( IC,  IR I/OUANNM 
QMTXLIKR,  I  *  A*»xM  (  IC,  IRI/OUANNM 

300  CONTINUE 

400  CONTINUE 

C  PRINT  OUT  NORMAL  1/(0  AND  IP  ANSPOSF  0  0  MATRICES 

PRINT  420 

420  FORMAT!  1H140X .3BHREAL  PART  OF  0 II  NORM AL ( 2F 0 , T» ANS POSF 0  I  I 
CALL  PR (NTM( QMTX (R, NR  ANY) 

PRINT  42R 

428  FORMAT! 1H140X , 38HRE Al  PARI  OF  «  ?  (  NOR M  AL  I  7 F  0  ,  T  R  ANSPOSF !3  )  I 
CALL  PRINTM( OMTXJO.NRANKI 
PRINT  436 

436  FORMAT!  1H140X , 38HMFAL  P AR  T  OF  0 3 ( NORM AL I 7 F  1 , T » ANS POSE n || 

CALL  PRINTM10MTXYR.NRANYI 
PRINT  444 

444  FnPMAT!  IH140X , 3BHPE AL  PART  OF  04(  NORMAL ! ZFO , TRANSPOSFOI  I 
CAIL  PPINTM! OMTXLR.NRANK) 

PRINT  4S2 

4S2  FORMAT!  1H13BX, 43HIMAE, INARY  PART  OE  01  (  NnR  MAL  I  /F  0  ,  T  P  ANS  0OSF  0 1  I 
CALL  PRINTM(QMTXI  I  , NR  AN Y 1 
PRINT  460 

460  FORMAT!  IHi 3BX ,4 3HIMAG INAP Y  PART  OF  02  t  NORMAL  I 7F0 , TRANSPOSFOI | 

CALL  PHINTM! OMTXJ I ,NR ANY  I 
PRINT  468 

468  FORMAT!  lHl38X,43H(MAr.|NARY  PARI  OF  03  (  NORMAL  I  7E  D  ,  TP  ANSPOSFO  II 
CALL  PR INTM( OMTXY ( ,NR ANY ) 

PRINT  476 

476  FORMAT! 1H1  38 X , 4 3H |M AG  I NAR Y  0 AR  T  OF  04 ( NORM AL I 7E0 , TR ANSPOSF 0 1  I 
CA' L  PR (NTM( QMTXI ( .NRANK | 

RF  TURN 
ENE) 


! 
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r  SURTYPF, FORTRAN, IMAP, l  STRAP 

C  A  ROUTINE  TO  PRINT  OUT  A  MATRIX  ARRAY 

SUBROUTINE  PRNOOTIP.NI 
DIMENSION  P< 1 ?0  * l 2 J  » 

NR  *  N 

Cl)  IOC  I  *  l.NH 
18  *  1 

20  It  *  I R£ 7 

I F  (  IF-NK ) 28. 28, 24 
24  IF  *  NR 
28  IF!  18-1136, 36,60 
36  PRINT  44,  |,(P(I,JI, J  =  IP, IF) 

44  FORMAT  I 5H0  ROWI3.2X.1PRF1S.6I 
GO  TO  80 

60  PRINT  68.IPII.JI.J  *  IP. IF) 

68  F0FMATI1H  9X» IP RF IS. 61 
80  18  =  I  EC  1 

I F ( I8-NR ) 20, 20, 10J 
100  CONTINUE 
RETURN 
END 
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T  SUBTYPE, ET1HTR  AN,  LMAP.l  STRAP 

f  THIS  ROUTINE  f.iir.UIATFS  r  Nil  POINT*  AND  SPACING  I  Nf  t',»  A  T  I  ( 

SURROUT INF  CALENP 
COMMON  OTR.RTO.CP I 

COMMON  /CMVCOM/  NM.CM | ( TO  I.CMV.KMV ,CM2 , FM.OF M, TWM, PMdOM 
COMMON  /F  NC.COM/  PNM|  l  C,(  61  I  ,  R  SSL  SP  (  6  1  )  ,(  NF  l|MN  t  Nil 

COMMON  /MTXCOM/  NOANK  .NRANK  I  , C «  T  XR  l  I  1 2  0 , 1  2  J  t  ,  C  MT  X  |  M  |  1  ? 0 , 1  ?(.>  I  ,  S PR MT 
1  x ( 1  20, 1201, CMxnpminO) 

common  /thtcim/  i nr ta.ntmfta.ol ihta, si nth.ccsth, i smul , i swtch'  m.sr 

1  MU  l  ,  SMUL  SSI  7  t  .COM  I  6), OHM,  NSFC  T  ,  NOP  S  I  6)  ,  FPP  SI  61  ,  K  SFC  T 
COMMON  /ROVCUM/  CKP , DC* R , CKP 7 ,C  SKR X , SNK K X , C ONK , MR  X T .ALPHA, I  BUOY, Oh 
1  .snalph.csalph 
SNAl  PH  *  S  INI  AL  PH A  t 
CS  A)  PH  -  CUSI  ALPHA  | 

OP  *  I  l,C-RR*TI*l 1. > SNAl  PHl/2. 

C  CAICUIATF  THF  FIRST  F  NO  POINT  AMI  STIP  SI /l 

TANGAM  =  SNAl  PM*CSAI  PH/ IOR-SNAI  PH* SNALPH) 

C.A“MA  =  ATANI TANGAM) 

IT IGAMMA I  20, 32, 3? 

20  CAMMA  -  f.AMMALCPI 
32  FPPSI 1 /  =  GAMMA 
COVO  *  NOPSI 1  I 
COHI  1  I  =  FPPSI 1  I/CO VO 

C  CALCUl AT  F  THF  SFCONO  FNO  POINT  ANT  STEP  SI/E. 

T  AN  PS  I  =  -RKXT*SNALPH*CSALPH/I  1.  -  JR-RRXT'-CSAI  PH*CSAl  PH) 

PSI  =  ATANI TANPSI ) 

IFIPSII  60,72,72 
60  PSI  =  PSICCPI 
72  FPPSI?)  =  PSI 
COVO  *  NOPSI 2 1 

COHI 2  I  =  IEPPSI 2)-rPPS( II  l/CDVO 
C  COMPUTE  THIRD  END  POINT  AND  STEP  Sl/F. 

FPPSI 31  *  CPI 
COVO  *  NOPSI 3 1 

C0HI3I  *  IEPPSI 31-FPPSI 2)  l/CDVO 

RETURN 

FNO 
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T  SURTYPF, FORTRAN, LMAP.l STRAP 

C  THIS  ROUTINE  COMPUTES  KR  AND  ITS  OERIVATIVF  AS  A  FUNCTION  OF  TmFTA 

SUHROUT IN E  GENKR 
COMMON  DT  R ,R  T  0, CP  I 

COMMON  /  C  MVCU  M/  NM.CM  mo  )tf.MV,KMV,CM?,EM,  JFM.TWM.PROrM 
COMMON  /FNCCOM/  PNMLLGI6U,RSSI  SPI  61  I  .CNFUMNIM  I 

COMMON  /MTXCOM/  NRANK  .NRANK  I  .CMTXRL  I  l’O.l  ??l  ♦CMTXIMI1?'!,  12;  I  .SPRMT 
1X1  12C.  120  I  ,CMXNRM(6'j) 

COMMON  /ThTCUM/  THE  TA,NTHETA,0|  THTA , SINTH.COSTH, I SMRL, ISMTCHIT) ,  S  R 
1MUI  , SMULSSI7I ,CnH(6l,nHM,NSECT,N0PSItl .FPPSI6I .KSFCT 
COMMON  /POVCUM/  CKR , OCX fl,CKR2,C SKP  X,SNKRX, CONK ,RR XT, ALPHA,  I  RODV.OP 
1,SNALPH,CS  Al  PH 
KSECT  *  KSFCT 

C  OF  T  EKM INF  SFCTION  FOP  INTEGRATION 

IF( KSECT- ? 14  3, 140, 240 
C  SECTION  1 

4 J  OU AN  1  =  SQPTI  1.0-<QB*SINTH/SNA1  ph> **21 
CKK  =  CONK  *  IQfl*C|lSTH/SNALPHC0OANl  I 

OCKR  •=  -CONK  •IOfl*SINTH/SNALPHt*  I  1.  OCOR^COSTH/I  SNALPH*OlJANll  I 
GO  TO  3'0 
C  SECTION  2 

140  QUAN2  =  THfTA-AlPHA 
SN0N2  *  SIM0UAN2I 
CKO  =  CONK*(  l.O-OM/SNON? 

nc.KR  r  -conk* 1 1. /'-cm  'r.ns (uoan'i  n  snon?*xnqn?» 

GO  TO  300 
C  SECTION  3 

243  OU An 3  *  (  l.O-ORXT-OHI /SNALPH 

CUNSO  *  SORT! RRXT*RRXT-IQUAN3*SINTH)**2I 
CKR  *  CONK*(OUNSO-OUAN3*COSTHI 

DCKR  *  CONK*(  OUAN3*SINTH-OUAN3*OUAN3*SlNTH*COSTH/(3UNSO) 

300  RETURN 
FND 


93 


BLANK  PAGE 


APPENDIX  ID 

A  NUMERICAL  EXAMPLE:  THE  SPHERE- CONE-SPHERE 
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•••••••••*•* •••••••••••••»•••••••••••••• 


INTEGRA! IONS/SECTION 
ENO  POINTS 


electromagnetic  scattering  from  general  axisymnetric  conductors 


CASES  MATRIX  RANK  SECTIONS  BODY  SHARE  <J  VECTOR 

*  A  )  4  46 

BODY  PARAMETERS  XIAI  BETA/RMO  ALPHA 

1.000  .794  IS. 000 

64  64  64 

ST. 8907  1)2. 6009  180.0000 
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ft  CAL  ft  A*  I  Oft  01(NORNALI2ED,TftANSftOSEO) 


M  *  0 


ROM  I 

-2*01  3760E-01 

-3.230234C-03 

-3.413709E-03 

1.213186E-04 

-1.278269E-09 

3.074439E-07 

ROM  2 

-3.230234E-03 

-1.98  3892E-02 

-4. 114393E-04 

-2.300333C-04 

9.430  3  34F-04 

-8, 106394E-07 

ROM  3 

-3.413T09E-03 

-4.114333E-04 

-8. 399868E-04 

- 1 . 340606E-05 

-7.489793E-06 

1.  120602E-07 

ROM  4 

1.213186E-04 

-2.  30  0333E-04 

-1.340606E-05 

-2.006002E-03 

-3.313340E-07 

-1.478181E-07 

ROM  5 

-1.278269E-03 

5. 430334E- 06 

-7. 983793F-06 

-3. 313340 E-07 

-3.098136E-07 

-4.608328F-0O 

ROM  6 

3.076439E-0T 

-8.  106394C-07 

1. 120602E-07 

-1.478181E— 07 

“4. 60 8 5 28 ft -09 

-3.3S7028E-09 

IMAGINARY  RART  OR 

01 (NORMAL IIEO. 

TRANSPOSE*)) 

ROM  I 

8.482023E-01 

-1.8413  77E-01 

-2.0B3360E-01 

4.  136367E400 

-1.U0176E401 

-1.441091E»02 

ROM  2 

•  -1.019649E-02 

6.  88031  7E-01 

8. 273333C-03 

-3.163230E-01 

2.262T33E»00 

3.443  13  7E*0  0 

ROM  3 

2.9S0026E-03 

-1.003273C-03 

6.042122E-01 

-3. 370680E-02 

-3.703436E-01 

3. 809894ft *00 

ROM  4 

-1.429732E-04 

1. 121333E-03 

-1.100414E-04 

S.7S2027E-01 

-4.1R4924E-02 

-2.319763ft-03 

ROM  3 

3.861338E-06 

-3.  128003E-03 

3. 32940 OE -04 

- 2. 4940 40 E- 04 

3 .646471F-01 

-7.883293ft-02 

ROM  6 

-6.9933T0E-08 

1.630344E-06 

-9.989378C-0S 

2.1101 16E-04 

-4.901 102E-04 

3.311401ft-01 

ftCAL  PART  Oft  04<NORNAL12E0, TRANSPOSED 


80M  1 

2.013760C-01 

3. 230234E-03 

3.413709E-03 

-1.213 1 86  E— 04 

1.278269E-03 

-S.0T6439E-0T 

ROM  2 

3.230234E-03 

1.  98  3892E-02 

4. 114333ft-04 

2.300333E-04 

-3.430334E-06 

8.106394E-07 

ROM  3 

3.41S709E-03 

4.  114333E-04 

8.399868E-04 

1 . 340606E-03 

7.983793F-04 

-1.120602E-0T 

ROM  4 

-1.213184E-04 

2. 500333E-04 

1.3406067-03 

2.09G002E-03 

3. 31 3340E-07 

1 .478 181 ft-07 

ROM  3 

1. 27826  9E-03 

-  3.  430334E- 06 

7.983793E-06 

3.  3 13340  E-07 

3 .098 136E-07 

4. 6083 2 8ft- 09 

ROM  6 

-3.076439E-07 

8.106394E-07 

-1.120602E-07 

1.478181F-07 

4.608328E-09 

3.337028E-09 

IMAGINARY  RART  OF 

04(NORMAL I2EO* 

TRANSPOSED) 

ROM  1 

1.31 7977E-01 

1.841377C-01 

2.083360E-01 

-4. 136367E400 

1 • 110 176E*01 

1.441091F402 

ROM  2 

1.019649E-02 

3. 11 9483E-01 

-8.273333E-03 

3.163230E-01 

-2.242755E400 

-3.443 137C*00 

ROM  3 

-2.980024E-03 

1.003273E-0S 

3, 937878E-01 

3.370680E-02 

3.703436F-01 

-3.  809894ft  *00 

ROM  4 

1.429732E-04 

- 1.  121S33E-03 

1.100414E-04 

4.247973E-01 

4.184924E-02 

2.319763F-03 

ROM  3 

-3.841338E-06 

3.  128003E-03 

-3. 329400E-04 

2.494044E-04 

4.333329E-01 

7.883293F-02 

ROM  4 

4.993370E-08 

- 1. 430344C-06 

9.989378E-06 

-2. 1 101 16E-04 

4.901 I02E-04 

4.488399E-01 
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M4TMII  Till, MU 


MM 

1 

I.IMIIK-02 

2,8114862-01 

1.0  125862  -  01 

-1.4801882-05 

4.0287782-06 

-1.6178217-07 

KM 

2 

MllMtf-01 

8.  6888882-04 

6.5402812-05 

8.1825152-06 

-4.7454702-08 

2.8517852  -  08 

Km 

3 

1.0UMU-0) 

6.2402812-05 

2.1124822-05 

-4.4  2  8  0  70  2-07 

8.8052222-08 

-2.4716682-08 

KM 

4 

-MHIIH-OI 

8.  18  25122-06 

-4.4280702-07 

1 .7180572-07 

-5.5848682-08 

5.8886287-10 

DOM 

3 

4.0207  7  82  -  06 

-4.  T4  54T02-08 

8.8052222-08 

-5.5848682-08 

5.4118152-10 

-2.5180517-11 

MOM 

* 

-i.»imie-oT 

2,8517852-08 

-2.6716682-08 

5.8886287-10 

-2.5180112-11 

2.1858812-12 

N4TMIX 

Tlll.lM*ClN6«r 

40M 

1 

2.  2840  7  72-01 

1.0227407-02 

4.2162522-01 

-1.6207157-04 

1.6787782-05 

-7.2208227-07 

ROM 

2 

1.0287402-02 

2  .  82  1  4  8  77  -  02 

6.22161  87  -  04 

1.7887157-04 

-8.2067472-06 

1.2451112-06 

MOM 

3 

4,2162822*01 

6.  2216182-04 

1.1741222-01 

2.8012702-05 

1.1218682-05 

-1.7782182-07 

KM 

* 

-1.6207152-04 

1.  7887157-04 

2.8012702-05 

1.4115187-05 

6. 4158082-07 

2.55C2467-07 

MOM 

3 

1.6  7  8  7  7  82  -  05 

-8  .  20  6  7  4  77  -  06 

1.1218682-05 

6.4158002-07 

5.4118182-07 

8.8668257-08 

MOM 

6 

-7.2205222-07 

1.2451112-06 

-1, TT821 82-07 

2. 550246f -07 

8.8668257-08 

5.8841177-08 

MiTMli  TUI, MU 


MOM 

1 

6.1010652-01 

-1.1876462-02 

8.8074012-01 

-1.8128612-04 

1.1175847-05 

-1,1 185812-06 

MOM 

2 

-1.18764  62-02 

4.1277212-01 

-8.2725042-05 

5.4212502-05 

-1.8871507-06 

1.6788647-07 

KM 

1 

8.8074012-01 

-8.2725042-05 

1.1118117-04 

-4.5758512-06 

4.6101562-07 

-1.1611847-08 

MOM 

4 

-1.8128612-04 

5.4212502-05 

-4.5758512-06 

8.0725522-07 

-1.0861287-08 

2.4618857-08 

MOM 

5 

1.1175842-05 

-1.5  8  71  502  -  04 

4.6101562-07 

-1.0861282-08 

2.2085812-08 

-8.8241617-11 

mom 

4 

-1.1185812-06 

1.6788647-07 

-1.1611842-08 

2.4618852-08 

-8.824J617-11 

7.7212817-12 

M6TM1X 

TI4I.IM6C1N6MY 

MOM 

1 

-4.8  241282-  01 

7.88104  87  -  01 

-6.8128652-01 

2.8664762-04 

-2.1451782-08 

8.1068887-07 

MOM 

2 

7. 88104  82-01 

-6.2421782-02 

-1.0505512-01 

-7.5182882-04 

1.6675802-05 

-2.1615617-06 

MOM 

1 

-4.8128652-01 

-  1.0805512-01 

-2.0720872-01 

-1.6858122-05 

-1.8806012-05 

2.  8TP 1  777-07 

Km 

4 

2.8644747-04 

-7.5182882-04 

-1.6858122-05 

-4.8084102-05 

-4.5114C87-07 

-1.8108402-07 

KM 

5 

-2.1451782  -  05 

1.4675802-05 

-1.8806012-05 

-6.5114082-07 

-7.2621202-07 

-8.1148117-08 

KM 

6 

8.1068882-07 

•2.  1615617-06 

2.8701772-07 

-1.5108402-07 

-8.1148112-08 

-7.6481517-08 
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HA  T  ft  I  X  Illl.fttU 


tow  1 

6.  8448  798  -02 

1.933149F-05 

8.4442998-04 

-2. 1761068-09 

2.909714F-06 

-1.079394F-07 

tow  2 

1.9331498-03 

1.4T2004F-03 

-6.1793098-06 

1. 1471198-09 

-2.36681tF-07 

3.0964408-08 

tow  \ 

8.644249F-04 

-6.  1T9309F-06 

1. 7201778-04 

-2.4476048-07 

9.999387F-0t 

- 1. 24400OF-09 

tow  4 

-2.  1  761068-09 

i.  147U9F-09 

-2.4476048-07 

1.3097628-07 

-3.376299F-09 

4.042027F-  10 

tow  5 

2.9097148-06 

- 2«  366tl tF-0 7 

9.9883878-08 

-3.3762998-09 

3.467906F-10 

-1.427183F-11 

tow  6 

-1.0793998-07 

3«C96540F-0t 

-1.244000F-09 

4.0420278-10 

-1.427183F-11 

1.378734F-12 

HATftlX 

TIlt.lHAGlNAftY 

tow  1 

2.909806F-01 

3.  79T793F-03 

. 3.29*8«28-03 

-9.7728068-09 

1.102171F-09 

-4.297141F-07 

tow  2 

3.74774  38-03 

2. 43C924F-02 

8.704882F-04 

3.1936948-04 

-8.9788627*06 

1.019428F-06 

tow  \ 

3.2968828-03 

t#  704892F -  04 

1. 2491188-03 

3.1914098-09 

1.190297F-08 

-1.46061  8F-07 

tow  4 

-9.  7  72  *068  -  09 

3. 19  3A94F- 04 

3.1914048-04 

3. 1 78  «8 1 8-04 

6.811860F-07 

2.  364842F-07 

tow  S 

1.1021718-04 

-9.978962F-06 

1. 190297F-08 

6.R11440F-07 

8.106828F-07 

9.  16282 7F*  09 

tow  6 

-4.2971418-07 

1 • 01 942 tF*  04 

-1.4606188-07 

2. 364*428-07 

9. 1 6?  827F-Q9 

8.6709468-09 

HATftlX  rm.ftFAl 


tow 

1 

8.  2201  338  -0  3 

-  2.  «6  21 2«8-03 

3.3929028-04 

-2.916401F-04 

2. 1264128-06 

-9.9«Ol«3F-0ft 

tow 

2 

1 .3412  2  08  -  02 

7.02  86948-04 

2.06396*8-04 

-4.9430948-06 

7.460290F-07 

-3.000  796F-0« 

tow 

3 

-9.0*99398-04 

-6.0476698-09 

-1.1 320698-09 

9.  1093B2F-07 

-4 . 1028*  98- Oft 

2.2337848-09 

tow 

4 

9.8149108-09 

3.  7090628-06 

7.3002628-07 

-2.9790138-0* 

2.4786978-09 

-1.2102618-10 

tow 

8 

-2.9734188-06 

-  1.4927098-07 

-2.7989788-0* 

1 , *8*1 60F-09 

-1.9734098-10 

7.3124038-  12 

tow 

6 

9,4329428-08 

7.4412398-09 

8.0931988-10 

-4.9167478-11 

2.341240F-12 

-1.7936*78-13 

HATftlX 

T 1 2 1 . 1 HAGINAftV 

tow 

1 

-2  .  7111478  -  02 

-1. 74  *  3  6  98  -  02 

3.291«llF-04 

- 1.292999F-04 

4.7813788-06 

-3.8  *69  9  28  -  07 

tow 

2 

-  2.  398436F-02 

2.  7891918-04 

-7.0673768-04 

2.4!l«208-04 

-7,817  4O0 F- 06 

1.4792998-07 

tow 

3 

1.3196718-03 

-8,  3*672  98-04 

1 , 7«9304F-09 

-1.9940918-04 

4,191 471F— 07 

-6,4677078-0« 

tow 

4 

- 1 . 090897F-04 

2.99  43928-09 

-1.42949*8-09 

3.9494938-07 

-1.89467*8-07 

4.1191048-09 

tow 

8 

4.1886218-06 

-3.  37887*8-06 

3.66261*8-07 

1.866«968-07 

3 , 3 16S96F-09 

-1.7479248-09 

tow 

6 

-1.9266748-07 

1.2974078-07 

-4.43614*8-0* 

4.0978988-09 

-1.6971228-09 

2.9464498-11 

106 


MATRIX  T(3I.R0Al 


SOW 

1 

S. 2201336-03 

1.  3312200-02 

-4.0*49340-04 

3. *139 lOt-05 

-2.37341 **-0* 

9.A3234K-0I 

sow 

2 

-2.8621286-03 

T.  0206941-0* 

-6. 0376*9* -03 

3. 7040620-06 

-1.432704t-07 

7.94l239f-09 

sow 

3 

3. 3929026-04 

2. 06  356*0-04 

-1.132063* -03 

7. 300262 f-OT 

-2. 73*37**-0* 

4. 03315*0-10 

sow 

4 

-2.9 1640 16-03 

-6.9*30530-06 

3. 1033*20-07 

-2.9T9013f-0* 

1 . *0* 160F-09 

-4.91674  77-11 

sow 

5 

2. 1264126-06 

7.4602900-07 

-4.102**30-0* 

2.47*6931-09 

-1.073409F-10 

2.331240f-12 

ROW 

6 

-9.9  8  918  36  -  08 

-3  .  00  0  7960  -  0* 

2. 2337*50-09 

-1.210261 *-10 

7 .312303*- 12 

-1.74*6*70-13 

MATRIX 

7(3). IMAGINARY 

ROW 

1 

-2.711147F-02 

-2.  39  04  36* -02 

1.3196710-03 

- 1.090*37f-04 

4. 180621F-O6 

-1.926*740-07 

ROW 

2 

-1. 7463696-02 

2.7*91410-04 

-*. 3*6723* -OA 

2.993392F-05 

-3. 77**70* —06 

1. 23760 TF-07 

ROW 

3 

3.23181  16-04 

-7.0633760-04 

1.7*33040-03 

- 1 .42949SE-03 

3.662616*— 07 

-5.4361460-08 

ROW 

4 

-1.2329336-04 

2.53102OF-O5 

-1.3940410-03 

3.349493E-0T 

-1. *66*66*- 07 

4. 037*9**- 09 

ROW 

•s 

3.  7813786-06 

- 3« 81 7400F-06 

4. 131671*-07 

— 1 . *9467**— 07 

3.316*96*— 09 

— 1.69712  2F -09 

ROW 

6 

-3.8869926-07 

1. 4742900-07 

-6.46270TT-0* 

4. 1 19103F-09 

-1.747323*- 09 

2.356344F-11 

MATRIX  TI41.R0AI 


ROM 

1 

2.  6*733  4* -01 

1. 64 2962F-02 

4. 209066C-03 

-1,1 19041*— 04 

1.4793000-0* 

-*.336*940-07 

ROM 

2 

1.642962F-02 

4.13R103F-O3 

2.7973720-04 

3.11*6410-05 

-4.3122200-0* 

*.  *3404*0-0* 

ROM 

) 

4.  20906 6F -03 

2.  79  7372F-04 

7.04*0*0*-03 

- 1 . 36063*0—06 

2 .6073310—07 

—7. 799»3Q*-04 

ROM 

4 

-1.1  1904 1*-04 

3.11*6410-05 

-1.36063*0 -06 

3.1119170-07 

-1.701116t-0* 

1.64*72*0-09 

ROM 

9 

1.4T9300*-03 

-6.312220F-0* 

2.607331*— 07 

-1.7013160-0* 

1.43*4740-09 

-7.022**4*-ll 

ROM 

6 

-3. 156594F-0T 

*•  854O90F-O* 

-7. 799*300-09 

1.69*72**- 09 

-7.022444*-U 

3.  910  1  930-12 

MATRIX 

7(4). IMAGINARY 

ROM 

1 

-4. 40  34  3  7*- 01 

-2.3102740-02 

-6. *6112*0 -03 

2.2731240-04 

-2.5*23070-09 

1.0100  340-06 

ROM 

2 

-2.31*274*-02 

-3.4931610-02 

-1.002 *39* -03 

-6. 51*9*3*— 04 

1 .434334*— 03 

-2.061**00-06 

ROM 

) 

-6.*6U2**-03 

-1.00  2*39*- 03 

-1.93*4900-03 

-3.*»3159t-03 

-1. 79**060-03 

2.2*34500-07 

ROM 

4 

2.  27  31240  -04 

-6.  31*9*3*-04 

-3.  **5139* -03 

-4.  3  237*1  *-03 

-7.673277F-0T 

-3. 21*0021-07 

WM 

9 

-2.3123O70-O3 

1. 4543341-05 

-1.7933660-03 

-7.6T327T*-0T 

-6.7324190-07 

-1.0443330-0* 

ROM 

4 

1.01Q034F-06 

-2. 04 1**0C- 04 

2. 2*3930*  -07 

-3.215002I-07 

-1.0943530-0* 

-7.1*44370-09 

107 


••*••*••••  i  -  m  »oj  sons  O3i»vmno3» 


>  9  9  < 

*  •  I*  « 


inr  n 

f*\  n  *■ 

9  9  9 


NOPOK'if  4  cNN*«co«<HOin04Pini»ifv«^«iroMnirr 
ostroeg  —  egir^9tf'9a39^'©99*'9f,*K99,o^*.099#'9 

Nc^c(ir9HiAKr9^f'<^H«cirHK(hc^f^4>M^H« 

eg  ®>  s  9  9  ciriNKxirK  Nr 

J5555SS5JS5555}J5557^SS5^7iI5J5?J'S^Ic 


r  —  * 

—  fs, 

•  •  f» 


OhOC 
K  »«* 
r-  —  it  it 
9  *  bc  9 
r  t  c  m 
9  9  ^  K 


?  S  r 
**  —  * 
«r  *r  « 

9  ir  tr 


I  I  I  i  I  I  I  I  I  •  I  I  i  I  I  I  I  I  t  i  I  I  I 


l  t  I  i  I  i  i  t  I  I  I  t  I  I  I 


oo< 

O  o  < 


>  o  o 

>  o  o 


o  o  < 
o  o  < 


iS 

•  • 

fg  eg 


J  Ubi  U>i  £  UJ  tti  k 

>  9  S  f-  9  9  4 

*  —  9  —  9 

—  9  eg  ^ 

-  n  r  r  N 


Ui  Uu  UJ  h 

9  9—* 


—  9 

?s 


7#? 

•  •  • 


*  * 
—  9* 
*  0* 


C  « 
9  9 
9  9 


O  —  i 
m  eg  i 
IT  K  i 

*r  ®  * 

fN-  *  i 


)  *e  9 
)  N  ^ 
>  <C  eg 
*  tr  9 

lit 


'  o  o  o  < 

>  o  o  o  < 

•  ♦  ♦  ♦  . 

»**! 
9  l*  * 

—  tr  9 

9  fsj  ^ 

«  —  9 

eg  »*  eg 


>  O  O 

>  o  o 


eg  « 

9  a 
9  a 
tr  r- 
—  9 


eg  O 
**  9 
-•  9  i 


a  K  9  ►  N  it 


•  -C  fg 

■  ® 
l  9 

•  9  9 


—  -#99 

—  K  f»  9 


O  N 
9  IT 
tr  — 


^  yj  ^  uj  yj  yj  ^  yj  yj 

r-  a^orNNa-fi 
9«rrg«90l*>K0 

ir  i*a<(«<oK  n 

g)inv«rg9r  «««• 


K9*»9MT- 


9  *•> 
eg  * 
tr  * 
«r  f* 

9 


e*  9  < 

9  9aiTir99#^^#<NN«gNMMVNfg 


eg  eg  eg  eg  f-g  eg 


94r¥>tr99®f-Ke-f»^e»9 


a  eg  9  9N  i 
—  f*-  Q  9 
#  fg  r  O  1 
tf>  •  9  —  < 

9  9  «r  e-  i 

9  9  9  9 


>  O  < 

3  iL  V 


«  o  o  o  o  o  < 


9  9  9  9 

)  fg  9  «  O 

i  h  h  rg  i 

;3£^g 

0^99 


9  O  V  K 


C  C  C  C  I 

I  I  I  I 

0  0  0  0  4 

?$f  ?' 

Hi  UlU  UM 
9  9  fg  9  4 
9  9  tr  O  ( 
9  —  9  *  I 
0^9*. 
e  O  i 


eg  —  fg  r*-  eg  9  , 
r* 

9 
9 
*» 
r> 

— 

1 


oooocoooooooooooooooo-  —  —  —  —  —  —  —  —  —  —  — 

OOOOOOOOOOOOOOOOOOCOOOOOOOOOOOOOO 
♦  ♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦♦♦♦  1  1  1  1  1  1  1  1  1  1  1  1 

m  u  ui  u.  in  ui  y  u.  y  u.  it/  y  ui  ui  y  i^i  y  y  Hi  r-  Hi  H1  H1  HUl1  tl 


>  m  9 

.  eg  « 

*  g  9 


o  « 

9  K 

9  4T 


9  — 

eg  9 
9  9 


>000 
>000 
♦  ♦  ♦ 
H.  Hi  Hi 

—  eg  — 

O  eg  eg 

*  fg  * 
c  ce  * 
9  9  9 


MT9NO^>,OKa-9  99e#-9e9 

SC  9BeCONf«Kx**ireCeN#e« 
— ■  e-  e-  K  990  — Oeeegir— 9^ert90'< 

0N09tf'9*'«>N  —  909egp*99C«o9« 


0000 
0000 
♦  ♦  ♦  ♦ 
Hi  Hi  Hi  HI 

nr  9Nir 
9  r-  —  »e 
i  fg  ®  e-  « 


I  I  I  I  I  I  I 


99ee9eN**09N-9f  ». 

9  « 


I  I  I 


I  I  I 


i  I 


lilt 


I  I  I 


o  o 
o  c 


'  tr  9 
V  m  e- 

t  9  * 

J  9  8 


oooooooooooooooooooooooooooooo 
ooocooooooooooocoooooooopoooop 
♦  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦* 
HiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHiHi 

99®  —  nnomt  9  —  oe(eoNeM9e«--c«N^4eN9 
fgege909£)9ir9  9#e«--Ofg99ef-Oe90K'^0 
egK094T«990Ke-e-  999CN-09N99-iMMf|if0» 
eM'g^-K9ae999tgeeNrgMeiroee-a9  9<^r»>9 
ei-NnN-9e9-Kf  O9NC9O#e9  99-9<0K99 
00—  —  —  —  —  00009  99®«KKK994rtr»ir»*99999 


I  I  I  I  I  I  I  I  I 


I  I  I  I 


9 

« 


3 

5 


8888! 
♦  ♦  ♦  ♦  * 


8888! 


JOOOOOOOOOOOOQOOOOOOOOOOOOOOOOp 

>0©0©000000000000c0^c0000000000 
♦  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦•  >♦♦♦♦♦♦♦♦♦♦ 
IHiHiHiHiHiHiHiHiHiHiHiHi^HiHiHiHiHiU.  ..  UittiUiltllfciltiWgUJUlWg 

'9e#9eaee«e9  9e9ee-eN«e#ir90#990e 

•  9trK994r994r94re.9tr— —  K9egfM99999fgOC  —  — 

)9^»eB9ef49ce9e909ee9eee9einr-f40f 

•  egffK9<fg--N99eaeg-e9»<«e-e  90Neif  B  9  M 

- -  -  - - 5«T9  —  oe»*^e*9  —  e  9  99egK9NN9aiT9 

- - - - *  ~ii  —  09SP*Kf» 


I  I  I  I  I  I  I  I 


Hi  Hi  Hi  Hi  h 
(>>099, 
—  —  9  O  r 

p  f*y  —  « 

fg  9  •  4T  f 

jtpss; 


0  0  0  0  4 
OOOOI 
♦  ♦  ♦  ♦  < 
HI  H  Hi  Hi 

O  9  4T  9 
9  9  9  9 

IT  f*  9  IT 
eg  O  9  f- 
eg  9  e-  eg 

O  —  eg  9 


O  O 
O  O 
♦  ♦ 


«r  BC  eg  4T  —  t*  e> 
9  eg  eg  9  e-  —  — 


1  9  f  ft  - 


9  O 
e>  — 
9  f* 
9  9 


I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 

OOOOOOOOOOOOOOOOOOCOOpOOOOOOOP 
OOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
♦  ♦♦♦*  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 
yj  yy  yj  y^  yj  yj  y  yj  yj  yj  yj  yj  y  yj  y  y^  yj  yj  y  yj  y^  y^  y  yg  y  yj  y  y 

ffeo9e>ieeBef9e9e^-eegBf^ftBOfMOf' 
irtrMr9e  9  f  e»f»e  9>9>*4*-e-  9egrg9f  99^NO<f-  — 

fMMTBf1'ffBir9e»o9eMteef9MN'H  mor 

eg  9  e>  9  9NpigN99MMVi»K  9  —  9^  -  *»90egf*tif>«9*f 
fteO-^Bi  e9-OMM‘9-^f9»Nf*9NNf#f9 
a  KrgN»"^ft  «tfgpg  —  9B099  e-?ffcr9M-09#f*e-N 


—  —  —  —  —  egpgegfg 


(*>9999999  *99m*'l**>*>®>fgfgpgfgegfgfg  —  —  —  —  — 


108 


i 

1 


»CCU«m»TEO  SUMS  F0« 


a 


u.  rgrs.rs.r'  *\  fs.  r  r,  r,  rv.  r\.  r\.  a  tv  r  iv  rv  r>  N^N^NNKNNNNNfvjNNNNNMNNNNNNNNN 

w  ercoforcc  c'onoorococ  r'ccooocor'oorrrrcr'cooccooooo 

o  ♦+♦♦•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

2  U.UJUjU<UjUjU.U.UjUJU.U.U-U.U 

<  o*-**4cw49.o*‘o*-r  «i?  >'ir^e‘«(irKN<oo940ir9^9f<hifiao>c«ff«^MMOM 

h-^irco^<<^NM«M<<;^K(rNV'MrtMV>V'<o«c9M4air«^v>«N09irir^nttN^9 

w"  4t>4*9C*fsn4c9*4C'tr*f*'<NC449'»«C4f90C4ir*99uOC49*449>C9*«»‘49j49‘4 

«  mnf -Mir<hf-c90*<fjr  •4*r  449-<^44C**9j94**4*r49-9-*4o*"9jfrii*44tr»rtrtr 

z  44444444*44^U'»rtru'tf,'irir4*r*ru'444444444<C449*9*9“9‘*9*9‘9-9-9“9-9* 

a  . . . . . . . 

I  I  »  »  i  i  i  *  i  »  i  i  i  I  i  •  i  I  I  t  i  i  i  i  *  i  i  I  I  l  l  i  i  i  i  i  i  i  i  i  i  •  i  i  l  i 


%A 

O 

a 


>oooooooooo< 


+  4  +  ■*  +  +  +  +  +  +  '4  4'4'4'4  +  tr\r*ie\ff\tr<r*t*r*t*r4tsir*f*r*tsi 


O 

4 

o t 


»ooocoaoooooo 


-  *  o 

'  4  r* 
k  O  4 


uj  u. 

9*  4 
C  fs< 

•-  o 

j  r 
*  «  — 
)  4  f* 


I 

uj  *. 
9*  »r 

O'  <*■ 

4  O 
O  4 
f-  — 


9*  O  9-  « 
*N«( 
4lT("l 
9U  4  IT  . 


K  r  CCK  z  c 


o  o 

11 
4  — 
9-  *« 
4  4 
9*  4 
4  4 
«  9- 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
I  I  I  II  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 
iLlilU'UIUtuwUUUUUUiWUlLliltLtLlk  UlUIUUlUJUwIbWw 

tr<oo*'oo*9tf'049*9*o4*-tro9**9©--9*tro4 
*"  4  9j  k  44* ^r»M*»pMpMN449^^9M«^cr»F«CM<NCN4 
n  4  N.4N49  fmfMNm^^mjNfMfmN^  4  4  O  N  O  4  9*  9*  44^ 
NN<ON9«mpM9  944^49  9INNBK  N^>  4  h  C  4  -«  *9 

N  N9  4f40^  ^944^4  4inNNmK4m<ONNlMMi«ltt4 
K^^r'H^h.4«^^C4fl404UCir^(rK4'(»'M-00 


ICC  B  • K  h  444CC4  44 

I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I  I 


till 


I  I  I  I  I  I  I  I 


OOOOOOCOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPOOO 
OCOOOOOOOOOOOOOOOOOOpOOOOOQOOOOOOOOOOOOOOOOOOO 
♦  ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦  + 
^  ^  yj  ^  yj  yj  yj  y^  yj  yj  m  y  ^  y  yj  y  y^  yj  yj  yj  y^  y^  yj  y  yj  yj  ^  yj  ^  yj  yj  yj  yj  yj  ^y  yj  yj  yj  ^  yj  yj  yj 

B("he<<B-o4iPi«O9h4*'Mf'(F|-4B4Niri*l49t^0®f*N0  0O4«-9N9B4 
4N9«D9r9i"4M('9N*4h-O-4ir4KM.f»44^irNl"4404W®Wi-h4N04e0 
4®NeN0444e4H»0^4Kf4B^494ir90-®Mf'-B®ff-^K-9lNi*i^9Ha 
04f'^4N4«»*B^F-irONN®»B4®HeeC4NN9'l*l4BO«CCF<4ir9‘4449»®4 
®B4®F-ir949BF-C4HMr«  mBNBC90-09^4-*K«94940B-K  4  **  4  9*  4  4 
■m  —  pu*Nf^^r'4*r»r4)^9-«<rO‘aoooOp--*M-ooOo<r«r®«r^h.9-444\*rtr4444 


1*111111 


<I<IIIII<II 


o 

; 


5} 

4  4 


9*  4 
IT  9* 
«  O  1 

o  o 

*r  4  1 


.  UJ  u.  t 

:zz\ 

'ire. 

'  o  o . 

■  4  -  * 

■  9*  4  I 


SO  1 
O  ( 
♦  ♦ 

UJ  X  1 

o  > 
o  *r  > 
*v»  rr 
O  *v» 


>  O  O  O  O  < 

>  O  O  O  O  < 


8000000000000 
000000000000 


9-  4 
*-  o 
»r  INP 

4  4 


*  4 
4  9* 
*T  4 
4  4 
9-  O 
O  — 


4yjyjyjyjyju.yjyjyjyju.yju.i 


4  *9  < 
9*  W  • 

S«  * 
IT  t 


m  in 
4  <r 
4  4 
4  4 


:444000— —  m* 


0O9N4-9-4C49W0 
44N^O«4K409-IPO 
«r*r«r4»-^h-4*rfs»tr4>M 
H-44N9IO*  444940 
004  BMTKIO^  *  4  tf\ 

NN--UU.M-000  4  4 


OOO 

OOO 


4  4  • 

4  •  *9 
M9N 
4  4  0 
—  •  4 
4  4  4 


sss 

♦  ♦  ♦ 

UJ  UJ  UJ 

—  —  4 

4  94  O 

—  40  «T 
4  0  0 
9I«4C 
4  4  N 


88888 


WUUWII 
9496* 
-99C- 
90  9M‘  ‘ 

«*o  O  1 

949H 


»  4  *0 

:h 


opopoooooooooopoooooooopooooooooooooooopooooop 

0000000000000000000000000000000000000000000000 

n  ^  yj  yj  y,  yj  ^  yj  yj  yj  yj  yi  y  yj  yj  ^  yj  y  yj  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y  y 

("N04«904  4990B**OBNir0OBN®4N®OOirN9'O9*ir099  4»o5440N9 
00  00«04p.fB**4irirt440  40ir4i44irO®04  4®9B49B4-40  40  00  0 
9  N  4  <  NC  -  MM«9®9NW9.9C'tr44'9«ir4»ir9‘l»l9'l9ir9-9»9N949'CMJ*JlP944N9 
S4<n090N9B4B®4lnB04B9N04C®4QN9944B949B494NQ90®99 
BB49990«N®4B9'®4444444BBIf'Bi44444449'ir49Nrt0f99®ir9 
NPNNNWW9(nf"(Fi^999J'J'999(n999Ji9JM»'9999999999999NNNNNN 


ooopooooooooooooQOooooooopopooooooooooopopOOOO 

OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOPOOOOOOOOOOOOO 

lillilUIUIUIUIUUIUUUUUpWUIU4UlWWU<IU  y  y  y  m  m  u*  ui  m  m  ui  iai  «1.  u  1  uj  uj  y  ui  m  U|  UI  y  m  ill 

i«094i'1f--C490i"*»4449*‘4N9  009n9n94®9  9®599®--49^P0iH 
h  4  IT  9>  hON4p^C  <0*\N04**  <909N44n9094900MP04®94NO^N9®h 
9  49eirircFMNBf*MF.irNir9M"^4iriririr4N«94irNC9NMrN4Circo9fir 
N90090040CBON444409N44f«f«44N904444NOCC04009009N 
JT4»T49-4  —  f^44  —  IT4  —  49-0*949*44  004  «9l9909fM-#fM-«C9lM-f9*4f4^ 
F*9-9-r-F-9*Cei)«D4C9C7‘C>00*-*-*-*-*-*-9ur<J*-*-*-*U*U— OOOe4444449-9»^9*9*9* 

-.  —  ,Z^_^-.J^INlNINru<N*N*N<NfMINfN<NIN<NiNfN<N(NIN-U  —  —  —  —  — 


OOOOOOOOOOOOOOO 
OOOOOOOOOOOOO  OO 
•  •••••••••••••a 

9J40449J404«<U4044 
4  4«r»T  444 


ooooooooooo< 

•  •  at  aaaaaaaaaaa 
N404<DN4044N4044 
9*9*44444000^^949494 


i  4  0  4  4* 
'*444* 


888808 
•  ••••• 

84  •  94  4  O 
4  4  9*  9*  4 


in 


*9 

*s> 


u 


109 


*1 


w  •  i 


REAL  PART  OR  01 INORNAllIED.TRANSPnSEOI 
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ROW 
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IMAGINARY  MART  OF 
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1.169331E-01 

-3.969337E»00 

2.20S317E*01 

ROW 

3 

-.OOOOOOE 

00 

6.430992F-03 

3.837S32E-CI 

-5.333773E-02 

4.11 3904E-0 1 

-1.241397E*00 

ROW 

4 

.OOOOOOE 

00 

5.295720E-04 

6.48768  7E-04 

3. 721704E-01 

-3.886486E-02 

3. 694933E- 01 

ROW 

5 

-.OOOOOOE 

00 

-  1 • 950544E -05 

2.9372«.2E-04 

-1.693886E-04 

5.623191E-01 

-6.428394E-02 

ROW 

6 

.OOOOOOE 

00 

1. 129103E-06 

-1.121  7-6E-05 

1.972807E-04 

-3.162321E-04 

3.481905F-01 

REAL  PART  OT  Q2I NORMAL  I  T  E  O.  TR ANSPOSEOI 


RP‘< 

1 

-.OOOOOOE 

00 

-.OOOOOOE  00 

-.OOOOOOE  00 

.OOOOOOE  00 

-.OOOOOOE  00 

.OOOOOOE  00 

R';m 

2 

-.OOOOOOE 

00 

1.  100635E-03 

3.630192E-04 

—6, 016391 E-06 

1.630004E-06 

-6.199754E-0B 

ROW 

3 

-.OOOOOOE 

00 

5.65C192E-04 

1 . 063367E-05 

8.89607BE-06 

-1.S17073E-07 

3.633249E-0B 

ROW 

4 

•  OOOOOOE 

00 

-6.016391E-06 

8. 8960 7 BE -06 

2 . 33B92 1  E-OB 

1. 301362E-07 

-1. 96790 3E-09 

ROW 

5 

-.OOOOOOE 

00 

1. 6*00041-06 

-1.31  7073E-07 

1.301362E-07 

-5.4B3412E-10 

1.2B4772E-09 

ROW 

6 

.OOOOOOE 

00 

-6.199T34E-08 

3.633249E-0S 

- 1 . 9 6790 3 E “09 

1.2B4772E-09 

-7.9903B3E- 1 2 

IMAGINARY  PART  OF 

02 1  NORMAL (2E0, 

TRANSPOSED  1 

ROW  1 

.OOOOOOE 

00 

•  OOOOOOE  00 

-.OOOOOOE  00 

.OOOOOOE  00 

-•OOOOOOE  00 

.OOOO00E  00 

ROW  2 

.OOOOOOE 

00 

-  1.074392E-02 

-4. 99B032E -03 

“1.741 323F-01 

-2.407710  E-01 

;. 70031 if»o; 

ROW  3 

.OOOOOOE 

00 

-4.99  8052E-03 

8.  377618E-04 

-1.6B7382E-03 

-1 •045580E-01 

-6.703992!— 01 

ROW  4 

-.OOOOOOE 

00 

2. 4549558-04 

-1.687381E-03 

7 . 143379E-04 

-6.836783E-04 

-7.37««.19E-02 

ROW  5 

.OOOOOOE 

00 

-  1 • 59  7296E-05 

7.47804  3E -03 

-6.836761 E-06 

3.91271BE-04 

-3.09M82E-04 

ROW  6 

-.OOOOOOE 

00 

6.6081  59E-07 

-4.177438E-06 

2.334003E-03 

-3.0980I3E-04 

2.027106E-04 
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MATRIX  Tilt, REAL 


ROW  1 

.00000  0* 

00 

-.000000*  00 

.000000*  00 

-.000000*  00 

.000000*  00 

-.000000*  00 

tow  2 

-.000000* 

00 

4.  5061642-04 

1.I54021E-05 

4.  174120E-04 

-1.548481E-0* 

1. 147400**00 

KOW  3 

.000000* 

00 

1.154021**05 

2. 241014* -06 

1.444441  f -07 

4.747446E-04 

2. 545754*- 10 

ROW  4 

-.000000* 

00 

4.  174120E-06 

1 • 464441 E— 07 

4.2256* >*-08 

-2.1026042-10 

1.174540*- 10 

ftOW  5 

.0000  00* 

00 

-1.54f4R}*-0« 

4.  747644* —04 

-2. 302404**10 

4.202441E-11 

-2. 162042*-12 

ROW  4 

-.000000* 

00 

I.  14740  0F-O0 

2. 545754F-10 

1. 176540**10 

-2.I42042E-12 

)«  7050402-15 

MATRIX 

TUI  .IMAGINARY 

ROW  I 

• OOOOOOF 

00 

-.000000*  00 

.000000*  00 

-.000000 *  00 

.OOOOOOf  00 

-.000000*  00 

ROW  2 

-.OOOOOOF 

00 

2.  1071812-02 

7.  441645*  -04 

1.485817F-04 

-2.181067F-06 

5.572444*-07 

ROW  3 

.00000  Of 

00 

7.  49  5655F-04 

4.  7R441 1*-04 

1.648742F-05 

8.440075**06 

-5.514081**08 

ROW  4 

-.OOOOOOF 

00 

1.95  58  37F-04 

1.648742* -05 

2.5186 70F-05 

7. 412458* -07 

1.  8255082-07 

ROW  5 

•OOOOOOf 

00 

-2.  18 106  7F-06 

8.440075E-06 

7.412458F-07 

4. 188411*— 07 

1.077144E-08 

ROW  6 

-.OOOOOOF 

00 

5. 572444E-07 

-5. 516O01F-O8 

1.825508F-O7 

1 .077164**08 

4.740585F-04 

MATRIX  TI2I.REAI 


ROW  1 

.OOOOOOF 

00 

.000000*  00 

-.000000*  00 

.OOOOOOF  00 

.OOOOOOF  00 

.OOOOOOF  00 

ROW  2 

-. OOOOOOF 

00 

2.8670442-05 

-0.  741 1 84f-06 

8.425469F-07 

-4.R51214F-0B 

4.542408E-04 

ROW  3 

•  OOOOOOF 

00 

2.6584142-05 

1.147160* -06 

2.051412E-O7 

-5.510068**09 

8. BB9B04E- 10 

ROW  4 

-.OOOOOOF 

00 

-1.146474**06 

-1.481712**07 

-7. 0811 41  *-04 

-1.0 11842* -10 

-4.294581E-12 

ROW  9 

.OOOOOOF 

00 

1.  185140*-07 

7.  651 1 81* -09 

1.077570 *-05 

-9.715518F-12 

2. 812226C- 12 

ROW  6 

-.OOOOOOF 

00 

-5.566215F-04 

-5.177404E-10 

”4,061 104*— 1 1 

— 1.574511F-11 

— 7.419104F— 14 

MATRIX  T 1 2  1 1 1  MAG  INARY 


ROW 

1 

•  OOOOOOf 

00 

•  OOOOOOf  00 

-.000000*  00 

.000000*  00 

.OOOOOOF  00 

•OOOOOOf  00 

ROW 

2 

-.OOOOOOF 

00 

-2.  31  7323F-03 

-6.611714* -04 

1.150744F-06 

— 1. 178168F-06 

1. 136766F-07 

ROW 

3 

•  OOOOOOf 

00 

-0.76V226F-O4 

-2.681184* -05 

-1.864  765*  -05 

2.467211f-07 

-7.0936180-00 

ROW 

4 

-•OOOOOOf 

00 

2.  22  3034F-09 

-1. 821 041C-05 

5.544575E-04 

-2.711 752E-07 

4. 127000F-O9 

ROW 

9 

•  OOOOOOf 

00 

-3.  6410060-06 

1.110814* -07 

-2.654008*-07 

1 . 447578**04 

-2.6904900-04 

ROW 

6 

-.OOOOOOf 

00 

1.3429730-07 

-7. 2 1441 2* -08 

4.24*4152-04 

-2.55*5l0*-04 

1.7044370-H 
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MATRIX  TDItllH 


ROM 

1 

•  00000 oc 

00 

-.000000*  00 

.000000*  00 

-.000000*  00 

.000000*  00 

-.090000*  00 

ROM 

2 

•eoooooc 

00 

2. 047O446-O9 

2.490419t-O9 

-1.194474*— 04 

1.10919OE-O7 

- 0. 044219*- 09 

ROM 

5 

-.000000 f 

00 

-R.TAURAE-O* 

1.147 140E— 04 

-1.441712C-07 

7 .491101*— 09 

—9. 1799040-10 

ROM 

4 

•  OOOOOOE 

00 

0.4294696-07 

2.0919126-07 

-7.001191 *-09 

1.O77970E-O9 

-4. 04110 90-11 

ROM 

5 

.000000* 

00 

-4.09121 AE-Ot 

-1.910040* -09 

-1.O11092E— 10 

-9.719910*- 12 

-1.979911E-11 

ROM 

4 

.000000* 

00 

4.9424016-04 

0.0040O4E-1O 

-4.294901E-12 

2.012224E-12 

-7.4191040-14 

MATRIX 

Till .imaginary 

ROM 

1 

.000000* 

00 

-.000000*  00 

.000000*  00 

-.000000*  00 

.000000*  00 

-.000000*  00 

ROM 

2 

.000000* 

00 

-2.11  T121E-01 

-0.T44224E-O4 

2.221019E-O9 

-1.4410O6C-O4 

1.1429710-07 

ROM 

9 

-.000000* 

00 

-6.6117196-04 

-2.401104* -09 

-1.021O41E-O9 

1.11O014E-O7 

-7.2144120-00 

ROM 

4 

.000000* 

00 

1. 1907496-06 

-1.444749E-09 

9. 9449 79 *-09 

-2.499009E-07 

4.2404  1  96  -  09 

ROM 

9 

1 

• 

o 

o 

s 

o 

o 

m 

00 

- 1.  1701  60E-O6 

2.96721 1C-07 

-2.711T92E-07 

1.497970E-09 

-2.9949106-09 

ROM 

4 

-•OOOOOOE 

00 

1.1 16766E-0  7 

-7.  09141 0E-O0 

4.1 27000*— 09 

-2.6904906-09 

1.7049176-11 

MATRIX  TUI, REAL 


ROM 

1 

.000000* 

00 

.0000006  00 

-.0000006  00 

.0000006  00 

.0000006  00 

.0000006  00 

ROM 

2 

.  000000* 

00 

1.1U91TE-01 

4.9991106-09 

1.1406216-09 

-1.1219696-07 

1.2942216-00 

ROM 

1 

-.0000006 

00 

4.  9991106-09 

9.  979140C-O4 

7.1109996-07 

1.2109276-00 

1.6966196-04 

ROM 

4 

.0000006 

00 

1.  14R421E-09 

7. 110999* -07 

1.2144466-07 

-0.0000116-10 

1.4700046-10 

ROM 

9 

.000000* 

00 

-  1 . 12 1949*-07 

1.21092TE-OR 

-0.0000116-10 

1.0606276-10 

-9.7211776-12 

ROM 

4 

.0000006 

00 

1.2942216-00 

1.4964196-09 

1.4700046-10 

-9.7211776-12 

1.041H26-12 

MATRIX  TIAI t IMAGINARY 


ROM 

1 

.000000* 

00 

.0000006  00 

-.000000*  00 

.0000006  00 

.0000006  00 

.0000006  00 

ROM 

2 

.000000* 

00 

-1.1174796-02 

-1.0907416-01 

-1.4217996-04 

4.1001*46-06 

-1.0000996-04 

ROM 

9 
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